Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,951,964	T→C	intergenic (+1407/‑237)	ilvX → / → ilvM	uncharacterized protein IlvX/acetolactate synthase II subunit IlvM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,951,964	0	T	C	98.6%	250.4 / ‑6.2	69	intergenic (+1407/‑237)	ilvX/ilvM	uncharacterized protein IlvX/acetolactate synthase II subunit IlvM
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (30/38);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGCGCCAGTAATTCAGAAATGTTGGAGAAATTATCATGATGCA  >  NC_000913/3951728‑3952208                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     ccGGAAAATTTCAGCCCCTCCATGGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGGTTGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCGCCGTAAAACGCAAGCA‑GTTACC‑GTTTAAAATCGTCTTACTCGAT‑AATCAACGGGTAGGGATGGTTCGATTT‑TGGCAGCAACTGTTTTTCCAGTAACtatatag                                                                                                                                                                                                                                        <  2:429795/251‑4 (MQ=255)  cGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGc                                                                                                                                                                                                                                       <  1:512591/251‑1 (MQ=255)   ggAAAATTTCATCACCTCCAGCGGCTTAGGTACCATTGGTTGTGGTTTACCGGCGTCGGTTGGCGCGCAAGGCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGTACCGTAAAACGCAAGCA‑TTTACC‑GTTTAAAATCTTCTTACTCGAT‑TACCAACGGTTAGGGAGGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCg                                                                                                                                                                                                                                      <  1:134587/251‑1 (MQ=255)     aaaaTTTCATCACCTCCAGCGGCTGAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGaa                                                                                                                                                                                                                                    <  2:351471/251‑1 (MQ=255)        aTTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGGTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGGGACGGCCCTTTCATGATGAATGTGCA‑AGAGCTGGGGCCCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AAACAACGGATAGGGATGGTTCGACAA‑TGGCAGCAAACGTTTTTTCAGGaaagaaatagcggaacc                                                                                                                                                                                                                                 >  1:223598/1‑236 (MQ=255)                aCCTCCCGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGCCGCGCGACCAAACGATACCGGAGCCTGTAACTCCGGAGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCCA‑CACGGTACC‑GGTGAAAATCGTCTTACTCGATCAACAAACG‑TTAGGGAAGGTTCG‑CAACTGGCAGCAACTGTTTTTCCAggcgagataagatcggcagagcatcgg                                                                                                                                                                                                                         >  2:478765/1‑226 (MQ=255)                aCCTCCAGCGGCTTAGGCACCATGGGTTTTGGTGTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATa                                                                                                                                                                                                                                           <  1:478765/233‑1 (MQ=255)                    ccAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGGCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGAt                                                                                                                                                                                                                     <  2:307672/251‑1 (MQ=255)                       ggggcttAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAAc                                                                                                                                                                                                                  <  2:223543/246‑1 (MQ=255)                         ggcttaagcaCCATGGGTTTTGGTTTACCGGCGGCGGTTGGGGCGGAAGACCCGCGTCCAAACGTTACCGTGGGCGGTATCTCCGGGGCCCGATCTTTCATGTTGGATGTGGACAG‑GCTGGGCACCGGAAAAACCAAGCA‑GTTACCAGTGGAAAA‑CGGCTTACCCGAT‑AACCCACGGTTAGGGGTGGTGTGACAA‑TGGCCACAACTGTTTTTCCACGGACGAGATAGCGGAACCACCCTTACCGATacacc                                                                                                                                                                                                                >  2:477413/10‑247 (MQ=255)                          gcttAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGGCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGGGCA‑AGAGCTGGGCACCGTAAAAAGCAAGCA‑GTTACC‑GTTGGAAAACGGCTTACTCGAT‑AACCAACGGTTAGGGGTGGGTTGAAAA‑TGGCAACAAATGGTTTTCCAGGGACGATATAGCGAAAACACCCTTACCGATAAcccc                                                                                                                                                                                                               >  2:153931/3‑251 (MQ=255)                                     aTGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCAt                                                                                                                                                                                                    <  1:188442/251‑1 (MQ=255)                                            ttGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGc                                                                                                                                                                                              <  1:421657/250‑1 (MQ=255)                                                   aCCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGcc                                                                                                                                                                                       <  1:476189/250‑1 (MQ=255)                                                      ggcggcggTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCg                                                                                                                                                                                   <  2:178976/251‑1 (MQ=255)                                                       gcggcggTTGGCGCGCAAGTCGCGCGACCAAACGATTCCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCTCCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCgg                                                                                                                                                                                  <  1:293882/251‑1 (MQ=255)                                                            ggTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATcc                                                                                                                                                                             >  1:128251/1‑251 (MQ=255)                                                                ggCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTgg                                                                                                                                                                         <  1:425910/251‑1 (MQ=255)                                                                     gcaAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGAGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGCGATGGTTCGACAA‑TGGCAGCCACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCTTGTTAGCCAGCGCCTTCGGCATCCCAGGCCAac                                                                                                                                                                    >  2:385785/2‑251 (MQ=255)                                                                             cgcgACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCAcc                                                                                                                                                            >  2:43731/1‑251 (MQ=255)                                                                             cgcgACCAAACGATACCGTCGTCTGTATCGCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGTCCAACACATCAcc                                                                                                                                                            <  1:263510/251‑1 (MQ=255)                                                                                   caaaCGATACCGTCGCCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCgtaa                                                                                                                                                       >  1:458041/3‑247 (MQ=255)                                                                                        gATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGAcc                                                                                                                                                  >  2:542309/1‑250 (MQ=255)                                                                                         aTACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCCACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCCCCCGTAACGAcc                                                                                                                                                  >  2:347633/1‑249 (MQ=255)                                                                                              gtcgTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGa                                                                                                                                           <  1:361016/248‑1 (MQ=255)                                                                                                 gTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCNAACACATCACCCGTAAAGACCAGGTTGAAGc                                                                                                                                        >  2:270090/1‑251 (MQ=255)                                                                                                      tATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCAc                                                                                                                                   >  1:129408/1‑251 (MQ=255)                                                                                                                gACGGATCTTCCATGATGAATGTGCA‑AGAGCTGGCCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATAACCCGTACAGACCAGGTTGAAGCGGCACTCGACACCAt                                                                                                                         >  2:159228/1‑251 (MQ=255)                                                                                                                   ggatttcgcATGATGGATGTGAA‑AGAGCTGGGCACTGTAAAACGCAAGCA‑GTGACC‑GGTGAAAACCGTCTGCCTCGAT‑AACCAGAGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACTGATACCCCCGATTTCATCATGTGAGCCAGCGCCGTCGGCCTCCTTGGCCAACACATCCTCCGTAAAGACCAGGTTTCAGCGGCACTCGACACCATGCt                                                                                                                      <  2:170828/243‑1 (MQ=255)                                                                                                                            atgatgAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTg                                                                                                              >  1:512417/1‑250 (MQ=255)                                                                                                                            atgatgAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTg                                                                                                              >  1:175640/1‑250 (MQ=255)                                                                                                                                 gaATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGc                                                                                                        >  1:39964/1‑251 (MQ=255)                                                                                                                                     gtgCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCAt                                                                                                     <  1:176203/250‑1 (MQ=255)                                                                                                                                                tGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCAt                                                                                          >  1:225458/1‑251 (MQ=255)                                                                                                                                                   gCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGtc                                                                                       >  1:154595/1‑251 (MQ=255)                                                                                                                                                       cGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCaa                                                                                   <  1:416683/251‑1 (MQ=255)                                                                                                                                                       cGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCa                                                                                    >  1:411778/1‑250 (MQ=255)                                                                                                                                                        gTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCCCAAt                                                                                  >  1:47677/1‑251 (MQ=255)                                                                                                                                                          aaaaCGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATc                                                                                 <  1:211567/250‑1 (MQ=255)                                                                                                                                                            aaCGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTNCTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCTACTGGTTTGCCAGGAACGATATAGCGAAACCACCCTTGCCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTCCATGTCCCAATcgac                                                                              <  2:234146/251‑1 (MQ=255)                                                                                                                                                             aCGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATcgacg                                                                             <  1:136572/251‑1 (MQ=255)                                                                                                                                                                cAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAAc                                                                          <  2:90575/251‑1 (MQ=255)                                                                                                                                                                 aaGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACt                                                                         >  2:133702/1‑251 (MQ=255)                                                                                                                                                                    cA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTg                                                                       >  1:13321/1‑250 (MQ=255)                                                                                                                                                                    cA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGAAAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCCTCCCTGGCCAACACATCACCCGTAAAGACCAGGATGAAGCGGCCCTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCCATCGACGAACTTg                                                                       >  2:288810/1‑250 (MQ=255)                                                                                                                                                                    cA‑GTTACC‑GTTGAAAATAGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATATTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGACACTCGACACCATTCTGACCAGTGATGGGCCATACCTGCTTCAGGTCTCAATCGACGCACTTg                                                                       >  2:299726/1‑250 (MQ=255)                                                                                                                                                                            c‑GTTGAACATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGTAACGATATAGCGAAACCCCCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACTCATTACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTc                                                               <  1:275858/251‑1 (MQ=255)                                                                                                                                                                            c‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTc                                                               <  1:351651/251‑1 (MQ=255)                                                                                                                                                                                    aaTCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGCTTTCCTCATGTTAGCCAGCGCCTTAGGCATCCCTGGCCAACACATCCCCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGCGCACGTCTGGCCGc                                                        >  2:375745/1‑251 (MQ=255)                                                                                                                                                                                       cGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTgg                                                     <  1:288810/251‑1 (MQ=255)                                                                                                                                                                                        gTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGt                                                    >  2:166784/1‑251 (MQ=255)                                                                                                                                                                                        gTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGt                                                    >  1:427930/1‑251 (MQ=255)                                                                                                                                                                                            tACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGCTTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTgccg                                                >  1:293832/1‑251 (MQ=255)                                                                                                                                                                                            tACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTgcc                                                 <  1:153931/250‑1 (MQ=255)                                                                                                                                                                                                  aT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTggg                                          >  1:335705/1‑250 (MQ=255)                                                                                                                                                                                                  aT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTggg                                          >  1:142694/1‑250 (MQ=255)                                                                                                                                                                                                         accGGTTAGGGATGGTTCGACAA‑TGGCAGCAACGGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGt                                    <  2:405724/249‑1 (MQ=255)                                                                                                                                                                                                          aCGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGATCGATATAGCGAAACCTCCCTTACCGATAACCCCGTTTTCCTCATGTTAGCCAGCGCCTTCGCCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGCCACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGATCTTTAGAACGTCTTGCCGCTTGTGCCGCCTGGTGCCAGTa                                   <  1:125082/251‑1 (MQ=255)                                                                                                                                                                                                             gTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAAtt                                <  1:170384/251‑1 (MQ=255)                                                                                                                                                                                                                 gggATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGGGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGa                            <  2:494076/251‑1 (MQ=255)                                                                                                                                                                                                                  ggATGGTTCGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGa                            >  1:85787/1‑250 (MQ=255)                                                                                                                                                                                                                         gcGACAA‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCAGACCTGCTTCATGTCTCAAGCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTgg                    <  2:489178/250‑1 (MQ=255)                                                                                                                                                                                                                              aa‑TGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGaaa               <  1:542309/251‑1 (MQ=255)                                                                                                                                                                                                                                    cagcaACTTGTTTTCCAGGAACGATATAGCGAAACCCCCCTTACCGATACCCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCCTCCCTGGCCAACACATCACCCGTAAAGACCAGGTTTAAGCGGCACTCGACACCATGCTGAACAGTGATGGTCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGGAATTCAGAAATGTTGGAGAAAGTATc          <  2:350562/251‑1 (MQ=255)                                                                                                                                                                                                                                     agcaACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCa         <  1:181873/251‑1 (MQ=255)                                                                                                                                                                                                                                       caACTGTTTTTCCAGTAACGAGATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGGCCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGTGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTTGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCatg       <  1:72534/251‑1 (MQ=255)                                                                                                                                                                                                                                       caACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCatg       <  1:292776/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCatgat     <  1:375745/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACAGCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGGCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCatgat     <  2:227192/251‑1 (MQ=255)                                                                                                                                                                                                                                             tttttcCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCAGCTTCATGTCTCAATCTACGAACTTGAGAACGTCTGGCCGATGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCATGATGCa  >  2:79916/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     CCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCA‑GTTACC‑GTTGAAAATCGTCTTACTCGAT‑AACCAACGGTTAGGGATGGTTCGACAA‑TGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGCGCCAGTAATTCAGAAATGTTGGAGAAATTATCATGATGCA  >  NC_000913/3951728‑3952208

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