Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,579,943	G→A	G103D (GGC→GAC)	tktB →	transketolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,943	0	G	A	97.3%	228.8 / NA	74	G103D (GGC→GAC)	tktB	transketolase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (41/31);  minor base C (1/1);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGGTTTTTACGATC  >  NC_000913/2579698‑2580182                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                cAACCCTGGTCATCCCGGCGCTCCGATGGGCATGGCAGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACAAAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACtata                                                                                                                                                                                                                                            <  1:376745/251‑1 (MQ=255)         gTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGc                                                                                                                                                                                                                                    <  1:449562/251‑1 (MQ=255)                   cgcgcCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAaccc                                                                                                                                                                                                                          >  1:288405/1‑250 (MQ=255)                      cccgattgtcATGTCTGATATTTCCGAACTGCTGTGGACCGATTTTCTGAATCATAACCCTACCGCCCCACCCCTGTTTGATCGCTCCCGCTTTATTCTTCCCACCGGTCACGCGTCGATGCTCCTCTACTGTTGGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCTCCCGGAGATTGACTATACGCCAGGCGTTGAAaccacca                                                                                                                                                                                                                       <  2:636417/242‑1 (MQ=255)                       ccGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGCAAACCCCAGGCCACCCGGAGATTGACTCTACGCCAGGCGTTGAAaccaccac                                                                                                                                                                                                                      >  2:3424/1‑251 (MQ=255)                               cATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCa                                                                                                                                                                                                                                       >  1:25840/1‑226 (MQ=255)                                  ggcggATATTGCCGTATTGCTGTGGAACGATTTTCTTAAACATACCCCTACTGACCCATCCTGGTATTATCGCGACCGCTTTCTTCTTTCCAACCGTCACGCGTCGATGCTGCTCTACAGTTTGCTCCTTCTGTCCGGTTACGACCTGTCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGTAGATTGCCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTg                                                                                                                                                                                                           <  2:550011/247‑1 (MQ=255)                                    cTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTGAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTCCAGTTTGCTACATCTGACCGGTTCCGCCCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCTTTCGAAAACCCGAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCAACACCGGCCCGCTTgg                                                                                                                                                                                                          >  2:265563/1‑250 (MQ=255)                                               gAAGCGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATTCTGACCGCTTTACTCCATCCACCGGTCCCGCGTCGCTGCTGCTCTACCCTTTGCTTCATTTGACCGGTTACGACCTGCCGCTGGACGAAATGGAGAATTTCCGTCAGCTACATTCGAAAACCCCAGGCCACCAGGAGCTTGACTATACGCCAGGCGTTGGACCCAGCACCGGCCCGCTTGGACCAGGCTTGGc                                                                                                                                                                                              >  1:653837/1‑251 (MQ=255)                                                aaGTGCTGTGGAACGATTTTCTTAAACATAACCCGACCGACCCAACCTGGTAGGATCGCGAACGCTTTATTCTTTCCAACGATCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGGTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCg                                                                                                                                                                                             <  2:361263/251‑1 (MQ=255)                                                               aTTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTg                                                                                                                                                                              >  1:457929/1‑251 (MQ=255)                                                                ttttCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTgg                                                                                                                                                                             >  1:559976/1‑251 (MQ=255)                                                                                ccTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTac                                                                                                                                                             >  1:650844/1‑251 (MQ=255)                                                                                   aCCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTCCAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCCCCACCGGccc                                                                                                                                                                                                                >  2:404194/1‑197 (MQ=255)                                                                                   aCCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGccc                                                                                                                                                                                                                <  1:404194/197‑1 (MQ=255)                                                                                         ccAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAGACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTggcggc                                                                                                                                                    <  1:564743/251‑1 (MQ=255)                                                                                         ccAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTggcggc                                                                                                                                                    >  1:464077/1‑251 (MQ=255)                                                                                          cAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACCACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGCACGTAAACTggcgg                                                                                                                                                     >  2:58479/1‑249 (MQ=255)                                                                                             ccTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACAGGCCCGCTTGGACAAAGTTTGGCGCACGCCGTCGGGCTGGCGATCGCAGAGCGTACACTGGCGgcgccg                                                                                                                                                >  2:16875/1‑249 (MQ=255)                                                                                                ggttagttCGCGACCGCTTTATTCAGTCCAACGGGCCCGCGTCGCTGCTTCTCTACAGTTTGCTACATCTTACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGGTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATTCGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGTCGAACGTCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGttt                                                                                                                                             <  2:83129/244‑1 (MQ=255)                                                                                                ggTATGATCGCGACCGCTTTATTCTGTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGtt                                                                                                                                              <  1:16875/250‑1 (MQ=255)                                                                                                   aTGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAc                                                                                                                                          >  1:488522/1‑251 (MQ=255)                                                                                                        cgcgACCGGTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGGTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAccagcc                                                                                                                                     <  2:328283/251‑1 (MQ=255)                                                                                                        cgcgACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAccagcc                                                                                                                                     <  1:505918/251‑1 (MQ=255)                                                                                                           gACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAacc                                                                                                                                                                                                                           >  2:235596/1‑162 (MQ=255)                                                                                                           gACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAacc                                                                                                                                                                                                                           <  1:235596/162‑1 (MQ=255)                                                                                                            aCCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGAc                                                                                                                                 >  1:676981/1‑251 (MQ=255)                                                                                                            aCCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGAc                                                                                                                                 <  1:651168/251‑1 (MQ=255)                                                                                                               gCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGAGTTGAAACCACCACCGGCCCGCTTGGACAAGGGTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTAAACTGGCGGCGCAGGTTAACCAGCCAGACCAt                                                                                                                              >  1:446996/1‑251 (MQ=255)                                                                                                                   aggctcttccgagcttCACGCGTCGATGCGGCTCTACAGTTTGCAACAACTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATGCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATgaga                                                                                                                          <  2:378523/236‑1 (MQ=255)                                                                                                                     ttCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATg                                                                                                                             >  1:122705/1‑246 (MQ=255)                                                                                                                           ccAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGAt                                                                                                                  <  1:554380/251‑1 (MQ=255)                                                                                                                           ccAACGGTCACGCGTCGATGCGGCCCTACAGTTCGTTACAACTGACCGGTTACGACCTTCCGCTGGAAGAACGGAAGAACTTCTGTCAGTTGCATTCGAAAACCCCAGTCCTCCCGGAGATTGACTATACGCCAGGCGTGGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGa                                                                                                                   <  1:106177/250‑1 (MQ=255)                                                                                                                              aCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGAGCAATTTAACCCGCCCGAcccggggatcctcggccac                                                                                                               >  2:282765/1‑234 (MQ=255)                                                                                                                                ggTCACGCGTCGATGCTGCTCTACAGTTTGCTACAGCTGACCGGTGACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACtt                                                                                                             <  1:403973/251‑1 (MQ=255)                                                                                                                                  tCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGCCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATgaga                                                                                                                          >  1:378523/1‑236 (MQ=255)                                                                                                                                  tCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGATCGTACACTGGCGGCGCAGGTTAACCAGCCAGACCATGAGATCGTCGATCCCTTc                                                                                                            >  2:169782/1‑250 (MQ=255)                                                                                                                                          cGATGCTGCTCTACAGTTTGCTACCTCTGACCGTTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGGCGATCACTTCACCTAtgt                                                                                                    <  2:676981/250‑1 (MQ=255)                                                                                                                                           gATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGCCCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGCCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGAAGTTTAACCAGCCAGACCAGGGGACCGTCGATACCTTCACCTATgtgt                                                                                                  >  2:22990/1‑251 (MQ=255)                                                                                                                                              gctgctCTTCAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGttt                                                                                                <  1:696785/250‑1 (MQ=255)                                                                                                                                                    ctACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCg                                                                                         >  1:500091/1‑251 (MQ=255)                                                                                                                                                      aCAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCg                                                                                         <  1:169782/249‑1 (MQ=255)                                                                                                                                                           ttGCTACATCTGACCGGTTACGCCCTGCCGCTGGAAGAACTGAAGAACTTCCATCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTAGACGCCAGGCGTTGAAACCCCCACCGGACCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCACCCAAGCCATGAGATCGTCGATCACTTCGCCTATTTGGTTCTTGGCAACGGCTg                                                                                  >  2:286429/1‑251 (MQ=255)                                                                                                                                                                acaccTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAAATTCCGTCAGGTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTg                                                                              <  2:488522/246‑1 (MQ=255)                                                                                                                                                                aaatcggacgggTTACGACCGGCCGCTGGAAGAACTGAAGAACGTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATCCGCCAGGCGTGGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGATGTCTGCCTg                                                                              <  2:296355/240‑1 (MQ=255)                                                                                                                                                                aCATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTg                                                                              >  1:266864/1‑250 (MQ=255)                                                                                                                                                                  aTCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACATACGTCAGTTGCATTCGAAAACCCCAGGCAACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAACAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGGCGGCTGCCTGATg                                                                           >  1:590372/1‑251 (MQ=255)                                                                                                                                                                      gACCGGGTACGCCCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGTCTGCCTGATGGAAg                                                                       <  2:457929/251‑1 (MQ=255)                                                                                                                                                                                    tGCCGCTGGAAGAACTGCAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTaa                                                                                                                                           >  1:225448/1‑169 (MQ=255)                                                                                                                                                                                    tGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTaa                                                                                                                                           <  2:225448/169‑1 (MQ=255)                                                                                                                                                                                        gCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCt                                                     >  1:591510/1‑251 (MQ=255)                                                                                                                                                                                           ggAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAg                                                                       >  1:670464/1‑230 (MQ=255)                                                                                                                                                                                           ggAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGtt                                                  >  1:323031/1‑251 (MQ=255)                                                                                                                                                                                           ggAAGAACTGAAGAACGTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACGATACGCCAGGCGTTGAAACCACCACCGGCCCGCGTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAg                                                                       <  2:670464/230‑1 (MQ=255)                                                                                                                                                                                                   tgaagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggc                                           >  1:607806/1‑250 (MQ=255)                                                                                                                                                                                                    gaagaactTCCATCATTTGCAGTCGAAACCCCCAGGCCACCCGGATTGTGACTATACGCCAGGCGGTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTTCACTGGCGGCGCAGTGTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCGGTTCGCTggc                                           <  2:50532/249‑1 (MQ=255)                                                                                                                                                                                                       gaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAAGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATATGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGACTGGTCGCTggcaggc                                       >  1:321543/1‑250 (MQ=255)                                                                                                                                                                                                           ttCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCAAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGGTTAACCAGCCCAACCATGAGATCGTCGATCACTTCACCTATGTGGTTATGGGCGACGGCTGCCTGCTGGGAGGGATTTCCCCCCAAGTATGTTCGCTGGCAGGCACgat                                  >  2:678080/1‑249 (MQ=255)                                                                                                                                                                                                            tCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATga                                                                                                                            <  1:35271/160‑1 (MQ=255)                                                                                                                                                                                                            tCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGGTTAACCAGCCAGACCATga                                                                                                                            >  2:35271/1‑160 (MQ=255)                                                                                                                                                                                                               gTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTggg                               >  2:554524/1‑250 (MQ=255)                                                                                                                                                                                                               gTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGGTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCACTGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTggg                               >  2:508267/1‑250 (MQ=255)                                                                                                                                                                                                                    ttGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCCGTTTAACCAGCCAGACAATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGGAGGTATTTCCCACGAAGTCTGTACGCTGGCAGGCACGCTGGGACTggg                         >  2:101029/1‑251 (MQ=255)                                                                                                                                                                                                                       cATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTTACCTATGTGTTTATTGGCGACGGCTGCCTGATGGACTGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGAATGGGCa                       >  1:21217/1‑250 (MQ=255)                                                                                                                                                                                                                        aTTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGt                                                                     <  1:408781/203‑1 (MQ=255)                                                                                                                                                                                                                        aTTCGAAAACCCCAGGCCACCCGGAGATTGAATATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGTGACTGCTGCCTGGTGGAAggg                                                                     >  2:408781/1‑202 (MQ=255)                                                                                                                                                                                                                           cGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCt                   <  1:634995/250‑1 (MQ=255)                                                                                                                                                                                                                                   ccAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGGttt          >  1:155347/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGGTTTTTAc      <  1:297462/251‑1 (MQ=255)                                                                                                                                                                                                                                       gCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGCATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGGTTTTTAc      >  1:92288/1‑251 (MQ=255)                                                                                                                                                                                                                                           cccGGAGATTGACTATACGCCAGGCGTTGAAACCACCCCCGGCCCGCTTGGACAAGGTTTGGCGATCGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGTGGTGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGGCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGGTTTTTACGATc  <  1:451210/251‑1 (MQ=255)                                                                                                                                                                                                                                            ccGGTGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGGTTTTTACGATc  >  1:225972/1‑250 (MQ=255)                                                                                                                                                                                                                                                  aTTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTcc                                                               <  2:545545/183‑1 (MQ=255)                                                                                                                                                                                                                                                  aTTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTcc                                                               >  1:545545/1‑183 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                CAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGGTTTTTACGATC  >  NC_000913/2579698‑2580182

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