Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,323,477	C→T	G198G (GGG→GGA)	glmM ←	phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,323,477	0	C	T	97.1%	223.7 / NA	69	G198G (GGG→GGA)	glmM	phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (35/32);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACCACGCAGCTGGCCCTGACGAAGACCTTCACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGT  >  NC_000913/3323234‑3323721                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     gcCACCACGCAGCTGGCCCTGACGAAGACCTTCACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGtt                                                                                                                                                                                                                                               >  1:15586/1‑251 (MQ=255)           agatGGCCCTGACGAAGACCTCCCCGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAACGCAATGCCGATATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCTTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCTGTTCGCGCAGCAc                                                                                                                                                                                                                                     <  2:398110/248‑1 (MQ=255)           aGCTGGCCCTGACGAAGACCTTCACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCAc                                                                                                                                                                                                                                     >  1:133762/1‑251 (MQ=255)               ggCCCTGACGAAGACCTTCACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTc                                                                                                                                                                                                                                 >  1:144934/1‑251 (MQ=255)                  ccTGACGAAGACCTTCACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCa                                                                                                                                                                                                                                                  >  2:7076/1‑231 (MQ=255)                  ccTGACGAAGACCTTCACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCa                                                                                                                                                                                                                                                  <  1:7076/231‑1 (MQ=255)                            aCCTTCACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGAt                                                                                                                                                                                                                    >  1:414036/1‑251 (MQ=255)                                 cACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTc                                                                                                                                                                                                               <  2:15586/251‑1 (MQ=255)                                 cACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGGCAACCATAATCACGCGATCTCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTc                                                                                                                                                                                                               <  1:171224/251‑1 (MQ=255)                                   cgcgcgATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGc                                                                                                                                                                                                             <  1:13848/251‑1 (MQ=255)                                   cgcgcgATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGc                                                                                                                                                                                                             <  1:202490/251‑1 (MQ=255)                                     cgcgATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCAc                                                                                                                                                                                                           >  1:233283/1‑251 (MQ=255)                                      gcgATGATATCCAAGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGAGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCAcc                                                                                                                                                                                                          >  2:27727/1‑251 (MQ=255)                                                  aTGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATc                                                                                                                                                                                              >  2:271154/1‑251 (MQ=255)                                                         gatcGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATccac                                                                                                                                                                                           <  1:394687/247‑1 (MQ=255)                                                         gatcGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCAGCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGTGATGTGATAAGTCGCCCCGTTTGCTCAATCccccac                                                                                                                                                                                        <  2:234195/250‑1 (MQ=255)                                                          atcGCCCTAGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCTCCATCGCCGTCGAAGGCAATAACGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTAGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGGGATAAGTCGCACCGTTTGCACAATCCACcaca                                                                                                                                                                                       >  2:238320/1‑250 (MQ=255)                                                              ccATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGCAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCtt                                                                                                                                                                                  >  2:68746/1‑251 (MQ=255)                                                               cATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTc                                                                                                                                                                                 <  1:157518/251‑1 (MQ=255)                                                               cATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTc                                                                                                                                                                                 <  1:343342/251‑1 (MQ=255)                                                                aTCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTc                                                                                                                                                                                 >  2:107340/1‑250 (MQ=255)                                                                aTCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTc                                                                                                                                                                                 >  2:79495/1‑250 (MQ=255)                                                                       ttATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACt                                                                                                                                                                         >  1:188150/1‑251 (MQ=255)                                                                       ttATTGCCTGCAGGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGATATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCGCCACAATCTTCAGTTCa                                                                                                                                                                           <  1:195282/249‑1 (MQ=255)                                                                             ccTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCCCTGAGGCt                                                                                                                                                                   >  1:237090/1‑251 (MQ=255)                                                                                tCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTa                                                                                                                                                                  <  2:23302/249‑1 (MQ=255)                                                                                    ggTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATCCCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttc                                                                                                                                                             <  2:408450/250‑1 (MQ=255)                                                                                     gccAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttcgt                                                                                                                                                           <  2:287933/249‑1 (MQ=255)                                                                                     gTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttcgt                                                                                                                                                           >  2:101434/1‑251 (MQ=255)                                                                                       cAACCATAATCACGCGATCGCCATCGCCGGCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAgttcgttc                                                                                                                                                         <  2:215351/251‑1 (MQ=255)                                                                                         aCCATAATCACGCGCTCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGAGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTNCgg                                                                                                                                                       >  2:376546/1‑251 (MQ=255)                                                                                          ccATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCggg                                                                                                                                                      >  1:98022/1‑251 (MQ=255)                                                                                              aaTCACGCGATCGCCATCGCCGGCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCCCACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAAc                                                                                                                                                   <  2:144934/250‑1 (MQ=255)                                                                                                   cgcgATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCt                                                                                                                                             >  1:345748/1‑251 (MQ=255)                                                                                                   cgcgATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCt                                                                                                                                             <  2:123201/251‑1 (MQ=255)                                                                                                    gcgNTCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCtt                                                                                                                                            <  2:378726/251‑1 (MQ=255)                                                                                                      gATCGCCATCGCCGTGGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTg                                                                                                                                          <  2:345748/251‑1 (MQ=255)                                                                                                        tcgccatagccGTCGAACGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGCGCGCGCGAACGTCGGTAGCCCCAACTTCGAAATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGCTAACGTTCGCTCCCAGTTCGCGCAGCAGGTTCGGCGCGATGTGATAAGTCGCACCGTTTACACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTAGGAAACGTGGCTTTGc                                                                                                                                         >  2:263302/1‑250 (MQ=255)                                                                                                             atcgccGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGAGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACt                                                                                                                                   >  1:255127/1‑251 (MQ=255)                                                                                                             accgccGTCGAAGGCAATACCGAGATCCGCTGTTTCAGCCAGCACCCGAGCCTGGTGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATATTTACACCGTTTGGCTCACATCCGATAGCGATATCGTTCGCTCCCAGTTCGCGCAGCGCGTTCTGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCGCCACAATCTTCAGTTCATTGAGGCTAATTTCTTTCGGGAACTTGGCTTTGCAAAACt                                                                                                                                   <  2:18685/249‑1 (MQ=255)                                                                                                                               agcGAGATCCGCATTTTCCGCCAGAACACGAGCCTGGAGCGCCCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAAAGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCCCAATCTTCAGTTCACTGCGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCCATATAGAGACCCGCgg                                                                                                                 >  2:271099/3‑251 (MQ=255)                                                                                                                                 cGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGc                                                                                                                >  1:18722/1‑250 (MQ=255)                                                                                                                                      tCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAAc                                                                                                          >  1:303109/1‑251 (MQ=255)                                                                                                                                               tCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCTGGAACTTGGCTTTGCAAAACTCGATATAGCGACTCGCGGCATCAACGATACggc                                                                                                  <  2:399368/250‑1 (MQ=255)                                                                                                                                               tCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggct                                                                                                 >  1:249375/1‑251 (MQ=255)                                                                                                                                                cagccagcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGAGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACggctg                                                                                                >  2:151874/1‑251 (MQ=255)                                                                                                                                                     agcACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCttt                                                                                           <  1:107340/251‑1 (MQ=255)                                                                                                                                                       cacaCGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTAc                                                                                         >  1:160683/1‑251 (MQ=255)                                                                                                                                                          acGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTAccc                                                                                       >  1:34402/1‑250 (MQ=255)                                                                                                                                                             ggCCTGTAGCGCGCGAACGTCTGTAGCCCCCACTTCGGCATTGATGTTTTCACCGTTTTGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCATGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCTCCACAATCTTCAGTTCACTGAGGCTAACTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACTATACGGCTGGCTTTACCCAGtt                                                                                   <  2:255127/250‑1 (MQ=255)                                                                                                                                                                   gAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCg                                                                             <  1:7334/251‑1 (MQ=255)                                                                                                                                                                      cgcgcgAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTGCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAAt                                                                          <  2:69969/251‑1 (MQ=255)                                                                                                                                                                        cgcgAACGTCGGTAGCCCCCACTTCGGCATTGATTTTTACCCCGGTTTGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAAGCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCa                                                                        <  2:251002/251‑1 (MQ=255)                                                                                                                                                                                       cccccACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTcc                                                         >  1:275099/1‑251 (MQ=255)                                                                                                                                                                                           cACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCtttt                                                     >  1:411189/1‑251 (MQ=255)                                                                                                                                                                                            aCTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTc                                                    <  2:133762/251‑1 (MQ=255)                                                                                                                                                                                                     ttGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGc                                           >  2:355270/1‑251 (MQ=255)                                                                                                                                                                                                      tGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCt                                          >  2:362395/1‑251 (MQ=255)                                                                                                                                                                                                        aTGTTTACTCCGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTc                                        <  2:403941/251‑1 (MQ=255)                                                                                                                                                                                                                  cGTTTGGCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCctct                              >  1:392970/1‑251 (MQ=255)                                                                                                                                                                                                                   gTTTGGATCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCtctt                             >  1:318738/1‑251 (MQ=255)                                                                                                                                                                                                                        gCTCACAACCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTAcc                        <  1:101434/251‑1 (MQ=255)                                                                                                                                                                                                                             caACCGATAGCGCTAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGCTATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCCTTTCCGCTTCGATGGCCTCTTCTACCGCATc                   >  2:233322/1‑251 (MQ=255)                                                                                                                                                                                                                               aCCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCg                 >  1:49992/1‑251 (MQ=255)                                                                                                                                                                                                                               aCCGATAGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCg                 <  1:68746/251‑1 (MQ=255)                                                                                                                                                                                                                                     aGCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGt            <  2:8849/250‑1 (MQ=255)                                                                                                                                                                                                                                      gCGATAACGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGttt          <  2:18722/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCg   <  1:334585/251‑1 (MQ=255)                                                                                                                                                                                                                                              gTTCGCTCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGAACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGt  >  2:231873/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     GCCACCACGCAGCTGGCCCTGACGAAGACCTTCACGCGCGATGATATACATGATCTGATCGCCATCGACTTTATTGCCTTCATGGTCAACCATAATCACGCGATCGCCATCGCCGTCGAAGGCAATACCGAGATCCGCTTTTTCAGCCAGCACACGAGCCTGGAGCGCGCGAACGTCGGTAGCCCCCACTTCGGCATTGATGTTTACACCGTTTGGCTCACAACCGATAGCGATAACGTTCGCCCCCAGTTCGCGCAGCACGTTCGGCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGT  >  NC_000913/3323234‑3323721

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