Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,925	T→A	64.7%	intergenic (‑1008/‑160)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,925	0	T	A	64.7%	89.5 / 81.4	91	intergenic (‑1008/‑160)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base A (29/29);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 8.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAT  >  CP000730/1929788‑1930034                                                                                                                                                 |                                                                                                              tGTATATTTGCCGATGATATAATAGGAGCGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GC‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                                            <  1:41819/141‑1 (MQ=255)  gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGACTGTGCTGTTGCGtattata                                                                                                          >  2:25995/12‑138 (MQ=255)   tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                                                                          <  1:396339/141‑1 (MQ=255)   tataTTCACAGATGGTATAATAATAGAGCCGCCTATCTCTCAGGCGTCAATT‑‑GACCAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                                         <  1:159110/141‑5 (MQ=255)     tattcactgatggtATAATAATAGATTCGCCTATCTCTGAGGCGTCAA‑‑TGGACGAAGAGAGGAGGTGCGT‑TAC‑T‑T‑T‑GC‑TAAGCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                                                                       <  1:178478/133‑7 (MQ=255)      attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                      >  2:327679/8‑141 (MQ=255)      attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑T‑T‑TACTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                      >  1:379275/8‑141 (MQ=255)        ttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑CGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAt                                                                                                     <  1:60499/141‑1 (MQ=255)         cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGCGCTGTTGCGTATTATACGTATTg                                                                                                   >  1:210517/5‑141 (MQ=255)         cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATCCTTATTg                                                                                                   >  1:134399/5‑141 (MQ=255)         cacagatggtATAATAATAGAGTCGCCTATCGCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                   >  1:331669/5‑141 (MQ=255)         cacagacggaaaaataatAGATTCGCCTATCATTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                   <  1:138513/129‑1 (MQ=255)          acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTGTTGCGTATTATACTTATTgg                                                                                                  <  1:85048/138‑1 (MQ=255)           cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTGC‑T‑T‑T‑GC‑TAATCATTTTCGTTCAGATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTCTTATACTTATTggt                                                                                                 >  2:223025/3‑140 (MQ=255)           cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                 >  1:154708/3‑141 (MQ=255)              atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                                              >  1:317410/5‑139 (MQ=255)               tggtataATCATAGGGCCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                             <  2:173884/138‑4 (MQ=255)               tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCACTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                             >  1:53428/4‑138 (MQ=255)               tggtagaataatAGAGTCGCTTATCTCTCAGGCGTCAATT‑‑GACGGAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TACTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                             <  2:69275/135‑4 (MQ=255)               tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGGCTaa                                                                                              >  2:63992/1‑141 (MQ=255)                ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGcctagt                                                                                            <  2:240761/139‑6 (MQ=255)                ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                            <  1:198195/139‑1 (MQ=255)                ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                            <  1:117363/139‑1 (MQ=255)                ggtataATAATAGAGTCGCCTATCGCTCAGTCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                            <  2:276286/139‑1 (MQ=255)                gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                                             >  2:323567/1‑141 (MQ=255)                gatATAATAGTAGAGTCGCCTATCTCTCAGGCGACAATTTAGACGCAGAGAGGAGGGG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCATCTTGGCTAAg                                                                                             >  2:211056/1‑141 (MQ=255)                 gtataATAATAGAGTCGCCTATCCCTCAGGCGCCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                           <  1:755/140‑1 (MQ=255)                 atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACAGCGTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                            >  2:217750/1‑141 (MQ=255)                  tataATAGTAGAGTCGCCTATCTCTTAGGCGTAAATTTAGACGCAGAGAGGAGGGG‑TG‑TAA‑G‑G‑G‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGtt                                                                                           >  2:273363/1‑140 (MQ=255)                  tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCCCATCATAGCACCAGTCATCAGTAGCTGTGCTGTTGCGTATGATACTTATTGGCTGAGTaa                                                                                          >  2:236414/1‑140 (MQ=255)                  tataATAATAGAGTCGCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGCGGAGGTG‑CA‑TTA‑G‑T‑TAGT‑TACTCATTTTGGCTTATATCATAGAACCAGTTCTCAGTGGCTGTGTGTTTGCTGTTTATACTTATTGGCTTAGTaa                                                                                          >  2:329658/1‑140 (MQ=255)                   ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACTAAGAGAGGAGGTG‑CA‑TTTCT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                         <  1:178075/141‑3 (MQ=255)                   ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGCAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACGAGTCATCAGTGGCTGTGCTGTTGCGGATTATACTTATTGGCTTAGTaaa                                                                                         >  2:86108/1‑139 (MQ=255)                   ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                         >  1:27226/1‑139 (MQ=255)                   ataataataGAGTCGCCTATCGCTCAGGCGTCAATT‑‑GACGAAGAGAGGGGGTG‑CA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                         <  1:99730/141‑3 (MQ=255)                    taatgatagAGTCGCCTATCTCTCAGGCGGCA‑‑GTGGACGAGGAGAGGAGGGGCAA‑TAC‑T‑T‑G‑GG‑TAATCATTTTCGTGCACATCATAGCACCAGTCAGCAGTGGCCGTGCGTTTGCGTATTCTAACTATTGGCttggtaacc                                                                                        >  2:355167/7‑133 (MQ=255)                      ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAA‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                      >  1:32389/5‑141 (MQ=255)                      ataataGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                      <  1:107586/141‑1 (MQ=255)                      ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAGTGGCTTAGTAAACGc                                                                                      <  1:294771/141‑1 (MQ=255)                      ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                      <  1:120251/141‑1 (MQ=255)                      ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAATCACTAAGTTGAcc                                                                                       >  1:148590/1‑140 (MQ=255)                      agaagagcggcGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTTTGCCGTTGCTTATTATACTTATTGGCTTAGTAAACGc                                                                                      <  2:143877/131‑1 (MQ=255)                       tagtagAGTCGCCGATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                      <  1:147474/141‑1 (MQ=255)                       atgtagAGTCGTCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                      <  1:267220/139‑1 (MQ=255)                       aaatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                     <  2:203161/138‑1 (MQ=255)                        aatcgagTCGCCTCTCTCTCAGGCGTCAATT‑‑GACGAAGCGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                    <  1:99205/137‑1 (MQ=255)                        aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                    >  1:69935/3‑141 (MQ=255)                        aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                    <  1:59507/139‑1 (MQ=255)                        aatagAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCCCCCATTAGTAGTGGCTGGTTTCTTTCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                    <  1:92973/139‑1 (MQ=255)                         atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                   <  1:61452/140‑2 (MQ=255)                           agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                 >  1:231743/1‑138 (MQ=255)                           agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                 <  2:360076/141‑4 (MQ=255)                           agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACAGCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTAGTATACTTCTTGGATTAGTAAACGCaacaa                                                                                 >  2:169649/1‑138 (MQ=255)                           agagTCCCCTATTTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                 <  2:10224/141‑4 (MQ=255)                            gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCGCCGGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                >  1:389272/1‑137 (MQ=255)                             agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                                                                                <  1:370019/141‑1 (MQ=255)                             agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                               <  2:175061/141‑6 (MQ=255)                             agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                               <  1:228455/141‑6 (MQ=255)                                cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GT‑TAATCATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTCAACGCaacaaataat                                                                            <  1:294894/141‑9 (MQ=255)                                cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                            >  1:32827/1‑133 (MQ=255)                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAGGCaacaaataatt                                                                           >  2:223362/1‑132 (MQ=255)                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGccacaaataatt                                                                           >  2:128686/1‑130 (MQ=255)                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                           >  1:353055/1‑132 (MQ=255)                                  ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGGAGCAATACCAAATa                                                                           >  2:350998/1‑141 (MQ=255)                                  ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                                           <  2:224689/141‑1 (MQ=255)                                  ccTATCTCTCAGGCATCAAATTAGAAGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                                           >  1:216671/1‑141 (MQ=255)                                   cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                         <  2:367696/141‑12 (MQ=255)                                    tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                        <  2:315908/141‑13 (MQ=255)                                    tATCTCTCAGGCGTCAACT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                        <  1:359212/141‑13 (MQ=255)                                     aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                                                        <  1:326067/141‑1 (MQ=255)                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCCTATTAGCTTAGTAAACGaaacaaataattaata                                                                       >  1:304464/1‑128 (MQ=255)                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                       >  2:349620/1‑128 (MQ=255)                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                       >  1:149819/1‑128 (MQ=255)                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                       >  2:134852/1‑128 (MQ=255)                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaaa                                                                       >  1:385276/1‑128 (MQ=255)                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattcatc                                                                       >  2:293934/1‑128 (MQ=255)                                     aTCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTG‑CA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                       >  2:121391/1‑128 (MQ=255)                                     aTCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTG‑CA‑TTA‑C‑T‑T‑TACTAATCATTTTCGCTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGTaacaaataattaata                                                                       >  1:240248/1‑128 (MQ=255)                                      tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                      <  1:323452/141‑15 (MQ=255)                                           cAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                                                  <  1:273782/141‑1 (MQ=255)                                            aGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                                                 <  2:262361/141‑1 (MQ=255)                                             ggCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACAtata                                                                >  1:159331/1‑141 (MQ=255)                                              gCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAg                                                               <  1:335541/141‑1 (MQ=255)                                                 ccAATGTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCg                                                            <  2:170324/140‑1 (MQ=255)                                                    aTTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTTAGTAGACGCAATATAAAATAGGTGACATATAGCCGCAc                                                         <  1:254576/141‑1 (MQ=255)                                                    aTTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAc                                                         >  1:307362/1‑141 (MQ=255)                                                          aCGCAGAGGTGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaa                                                   <  2:220766/141‑1 (MQ=255)                                                          aCGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaa                                                   <  2:287627/141‑1 (MQ=255)                                                                 gaggagGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATccc                                            <  1:86747/141‑1 (MQ=255)                                                                                             cTTATTTTCGTTCATATCATAGCACCAGTCATCAGTGGCTGTGCCATCGCGTTTTTTTATTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGTACCCATAAAAATCCCCTCACTACCGCAAATAGTGAg                       <  2:101414/141‑1 (MQ=255)                                                                                              ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACACAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                      <  2:121359/141‑1 (MQ=255)                                                                                                aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATCGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                    <  2:203355/141‑1 (MQ=255)                                                                                                            atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa        >  1:198931/1‑141 (MQ=255)                                                                                                            atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa        >  1:349501/1‑141 (MQ=255)                                                                                                            atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa        >  1:259325/1‑141 (MQ=255)                                                                                                                  gCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATACAAAATAGGTGACAAATAGCCGCGCCAATTAAGATCTCCTCACTACAGCAAATAGTGACGGGATGGGGGTGAAAGTAAAt  >  1:128711/1‑141 (MQ=255)                                                                                                                                                 |                                                                                                              TGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAT  >  CP000730/1929788‑1930034

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑