Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,941	A→T	59.1%	intergenic (‑1024/‑144)	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,941	0	A	T	59.1%	40.6 / 85.8	77	intergenic (‑1024/‑144)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/15);  new base T (20/24);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

TGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTT  >  CP000730/1929802‑1930074                                                                                                                                                 |                                                                                                                                      tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGGCTaa                                                                                                                                      >  2:63992/1‑141 (MQ=255)  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                    <  1:198195/139‑1 (MQ=255)  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                    <  1:117363/139‑1 (MQ=255)  ggtataATAATAGAGTCGCCTATCGCTCAGTCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                    <  2:276286/139‑1 (MQ=255)  gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                                                                                     >  2:323567/1‑141 (MQ=255)  gatATAATAGTAGAGTCGCCTATCTCTCAGGCGACAATTTAGACGCAGAGAGGAGGGGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCATCTTGGCTAAg                                                                                                                                     >  2:211056/1‑141 (MQ=255)   gtataATAATAGAGTCGCCTATCCCTCAGGCGCCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                                   <  1:755/140‑1 (MQ=255)   atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACAGCGTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                                    >  2:217750/1‑141 (MQ=255)    tataATAGTAGAGTCGCCTATCTCTTAGGCGTAAATTTAGACGCAGAGAGGAGGGGTG‑TAA‑GG‑G‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGtt                                                                                                                                   >  2:273363/1‑140 (MQ=255)    tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCCCATCATAGCACCAGTCATCAGTAGCTGTGCTGTTGCGTATGATACTTATTGGCTGAGTaa                                                                                                                                  >  2:236414/1‑140 (MQ=255)    tataATAATAGAGTCGCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGCGGAGGTGCA‑TTA‑GT‑TAGT‑TACTCATTTTGGCTTATATCATAGAACCAGTTCTCAGTGGCTGTGTGTTTGCTGTTTATACTTATTGGCTTAGTaa                                                                                                                                  >  2:329658/1‑140 (MQ=255)     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACTAAGAGAGGAGGTGCA‑TTTCTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                 <  1:178075/141‑3 (MQ=255)     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGCAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACGAGTCATCAGTGGCTGTGCTGTTGCGGATTATACTTATTGGCTTAGTaaa                                                                                                                                 >  2:86108/1‑139 (MQ=255)     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                 >  1:27226/1‑139 (MQ=255)     ataataataGAGTCGCCTATCGCTCAGGCGTCAATT‑‑GACGAAGAGAGGGGGTGCA‑TTA‑CTTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                 <  1:99730/141‑3 (MQ=255)        ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAA‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                              >  1:32389/5‑141 (MQ=255)        ataataGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                              <  1:107586/141‑1 (MQ=255)        ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAGTGGCTTAGTAAACGc                                                                                                                              <  1:294771/141‑1 (MQ=255)        ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                              <  1:120251/141‑1 (MQ=255)        ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAATCACTAAGTTGAcc                                                                                                                               >  1:148590/1‑140 (MQ=255)        agaagagcggcGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTTTGCCGTTGCTTATTATACTTATTGGCTTAGTAAACGc                                                                                                                              <  2:143877/131‑1 (MQ=255)         tagtagAGTCGCCGATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                              <  1:147474/141‑1 (MQ=255)         atgtagAGTCGTCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                              <  1:267220/139‑1 (MQ=255)         aaatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                             <  2:203161/138‑1 (MQ=255)          aatcgagTCGCCTCTCTCTCAGGCGTCAATT‑‑GACGAAGCGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                            <  1:99205/137‑1 (MQ=255)          aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                            >  1:69935/3‑141 (MQ=255)          aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                            <  1:59507/139‑1 (MQ=255)          aatagAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCCCCCATTAGTAGTGGCTGGTTTCTTTCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                            <  1:92973/139‑1 (MQ=255)           atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                           <  1:61452/140‑2 (MQ=255)             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                         >  1:231743/1‑138 (MQ=255)             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                         <  2:360076/141‑4 (MQ=255)             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACAGCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTAGTATACTTCTTGGATTAGTAAACGCaacaa                                                                                                                         >  2:169649/1‑138 (MQ=255)             agagTCCCCTATTTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                         <  2:10224/141‑4 (MQ=255)              gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCGCCGGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                        >  1:389272/1‑137 (MQ=255)               agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                                                                                                                        <  1:370019/141‑1 (MQ=255)               agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                       <  2:175061/141‑6 (MQ=255)               agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                       <  1:228455/141‑6 (MQ=255)                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GT‑TAATCATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTCAACGCaacaaataat                                                                                                                    <  1:294894/141‑9 (MQ=255)                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                                    >  1:32827/1‑133 (MQ=255)                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAGGCaacaaataatt                                                                                                                   >  2:223362/1‑132 (MQ=255)                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGccacaaataatt                                                                                                                   >  2:128686/1‑130 (MQ=255)                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                   >  1:353055/1‑132 (MQ=255)                    ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGGAGCAATACCAAATa                                                                                                                   >  2:350998/1‑141 (MQ=255)                    ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                                                                                   <  2:224689/141‑1 (MQ=255)                    ccTATCTCTCAGGCATCAAATTAGAAGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                                                                                   >  1:216671/1‑141 (MQ=255)                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                 <  2:367696/141‑12 (MQ=255)                      tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                <  2:315908/141‑13 (MQ=255)                      tATCTCTCAGGCGTCAACT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                <  1:359212/141‑13 (MQ=255)                       aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                                                                                                <  1:326067/141‑1 (MQ=255)                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCCTATTAGCTTAGTAAACGaaacaaataattaata                                                                                                               >  1:304464/1‑128 (MQ=255)                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                               >  2:349620/1‑128 (MQ=255)                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                               >  2:134852/1‑128 (MQ=255)                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                               >  1:149819/1‑128 (MQ=255)                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaaa                                                                                                               >  1:385276/1‑128 (MQ=255)                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattcatc                                                                                                               >  2:293934/1‑128 (MQ=255)                       aTCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTGCA‑TTACTT‑T‑GC‑TAATCATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                               >  2:121391/1‑128 (MQ=255)                       aTCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTGCA‑TTA‑CT‑T‑TACTAATCATTTTCGCTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGTaacaaataattaata                                                                                                               >  1:240248/1‑128 (MQ=255)                        tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                              <  1:323452/141‑15 (MQ=255)                             cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                                                                                          <  1:273782/141‑1 (MQ=255)                              aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                                                                                         <  2:262361/141‑1 (MQ=255)                               ggCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACAtata                                                                                                        >  1:159331/1‑141 (MQ=255)                                gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAg                                                                                                       <  1:335541/141‑1 (MQ=255)                                   ccAATGTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCg                                                                                                    <  2:170324/140‑1 (MQ=255)                                      aTTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTTAGTAGACGCAATATAAAATAGGTGACATATAGCCGCAc                                                                                                 <  1:254576/141‑1 (MQ=255)                                      aTTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAc                                                                                                 >  1:307362/1‑141 (MQ=255)                                            aCGCAGAGGTGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaa                                                                                           <  2:220766/141‑1 (MQ=255)                                            aCGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaa                                                                                           <  2:287627/141‑1 (MQ=255)                                                   gaggagGTGTA‑TAA‑GG‑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:349501/1‑141 (MQ=255)                                                                                            atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa                                                >  1:259325/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)                                                                                                                                          tttttcTAAGTAGACGCCATACAAAATAGGTGACATATAGCCGCACCGATAAAAATCCCCTCACTCCCGCATATATTGTGGGGTTTGGTGTATAAGTAAAAACTTATTTTCGTTGTCTTAATTATACTGCTAAttttttct  >  2:301139/6‑138 (MQ=255)                                                                                                                                                 |                                                                                                                                      TGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTT  >  CP000730/1929802‑1930074

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