Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,930	T→A	65.7%	intergenic (‑1013/‑155)	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,930	0	T	A	65.7%	98.2 / 83.1	91	intergenic (‑1013/‑155)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/17);  new base A (29/29);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTT  >  CP000730/1929790‑1930070                                                                                                                                                    |                                                                                                                                             tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                                                                                                              <  1:396339/141‑1 (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                                                                                                           >  1:149819/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                                                                                                           >  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:259325/1‑141 (MQ=255)                                                                                                          atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa                                            >  1:349501/1‑141 (MQ=255)                                                                                                          atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa                                            >  1:198931/1‑141 (MQ=255)                                                                                                                gCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATACAAAATAGGTGACAAATAGCCGCGCCAATTAAGATCTCCTCACTACAGCAAATAGTGACGGGATGGGGGTGAAAGTAAAt                                      >  1:128711/1‑141 (MQ=255)                                                                                                                                           gttggttttttttttttGATAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAACAATCAAAATGCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTATATACTTATGTTCGTTGTCTTAATtatatt           >  2:299475/6‑139 (MQ=255)                                                                                                                                                    tCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtttt  >  1:166392/1‑141 (MQ=255)                                                                                                                                                    |                                                                                                                                             TATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTT  >  CP000730/1929790‑1930070

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: ‑