Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,580,388	G→A	P251P (CCG→CCA)	tktB →	transketolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,580,388	0	G	A	95.6%	140.1 / NA	45	P251P (CCG→CCA)	tktB	transketolase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/26);  minor base C (2/0);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GGCACGCTGGGACTGGGCAAGCTGATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCGAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACGGAT  >  NC_000913/2580143‑2580618                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       ggtaCGCTGGGACTGGGCAAGCTTATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCCAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATaa                                                                                                                                                                                                                                   <  2:171113/248‑1 (MQ=255)      gCTGGGACTGGGCAAGCTGATTGGTTTTTACGATCACAACGGTTTTGCCATCGACGGTGAAACAGAAGGCTGGTTTACCGTCGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATTCGCGGGCGGTGAAGGAAGCTATTCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAg                                                                                                                                                                                                                              <  1:357961/251‑1 (MQ=255)              tGGGCAAGCTGATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATAGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATATGCCGTAAGGGGATAGGATTTGGGTCGCCAAATAAAGCAGGTAAGGa                                                                                                                                                                                                                       >  1:271016/1‑250 (MQ=255)                 gCAAGCCGATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAag                                                                                                                                                                                                                     >  1:118976/1‑249 (MQ=255)                   aaGCTGATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTGTGGTTCGCCAAATAAAGCAGGTAAGGAAGAggc                                                                                                                                                                                                                 <  1:23459/251‑1 (MQ=255)                                   aCGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTg                                                                                                                                                                                                 >  1:349079/1‑251 (MQ=255)                                   aCGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTGCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTg                                                                                                                                                                                                 <  2:204095/251‑1 (MQ=255)                                      aTCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTgggg                                                                                                                                                                                              <  1:376089/251‑1 (MQ=255)                                          caACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGgaag                                                                                                                                                                                          <  1:71031/251‑1 (MQ=255)                                            aCGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGAAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCAAGATCCGCAGGCGGTGGAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGGGATTGGCTTTGGGTCGCCAAATAAAGCAGGTAAGGCAGAGGCGCACGGCGCACCACTGGGGgacgc                                                                                                                                                                                         >  1:217043/1‑247 (MQ=255)                                              ggaaTTTCCATCGACGGGGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGAGCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGgaagaaga                                                                                                                                                                                      <  2:153928/248‑1 (MQ=255)                                               gTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAACGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGGTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCCCGGAGCACCACTGGGGgaagcaga                                                                                                                                                                                      >  1:68297/1‑250 (MQ=255)                                                       aTCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCt                                                                                                                                                                             <  1:198296/251‑1 (MQ=255)                                                                gAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGTTTATCTGCCGTACGGTGATTGGCTTTGGGTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGGGGCGCTGGCACGGc                                                                                                                                                                     >  1:138862/1‑250 (MQ=255)                                                                gAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGc                                                                                                                                                                     >  1:371104/1‑250 (MQ=255)                                                                 aaaCAGAAGGCTGGGTTACCGAAGATACGTCAAAACGTTTTGAAGCCTATAACTGGCATGTGATCCATGAAATCGACGGTCACGATCAGCAGGCGGTGAAGGACAAGATCCGTGAAGCGCTAAGCGTGAAAGATACGACGTCGCTGATTTTCTGACGTACGGGTATTGCCTTTTGGTAGCACAAAAAAGAAGGTAAGGAAGAGGCGCACGGCGCACCCCTGGGGGAAAAAGAAGGGGCGCTgggcaggca                                                                                                                                                                    >  2:305542/1‑243 (MQ=255)                                                                       aaGGCTGGTTTACCGACGATCCGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATACGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAAATAAGCCGTCGCTGATTAACTGCCGTACGGTGATTGGCCTTTGTTAGCCAAATAAAGCAGGTAAGGAAGAGGCGCCCGGCGCACCCCTGGGGGCAGCCGACGTGGCGCTGGCACgggccaaac                                                                                                                                                               >  2:308262/1‑242 (MQ=255)                                                                            tGGTTTACCGACGATACGACAAAACGATTTGAAGCCAACCACTGGCATGTGCTCCATGAAATCGACGGTCACGATCCGCAGGCGGTGCAGGAAGCCATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGGCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGACACATAAAGAAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGGAGAAGAAGTGGCGCTGGCAAGGCAAAAACAGGGCTg                                                                                                                                                        >  2:383796/1‑251 (MQ=255)                                                                             ggTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTgg                                                                                                                                                       <  1:404316/251‑1 (MQ=255)                                                                             ggTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTg                                                                                                                                                        <  2:118976/250‑1 (MQ=255)                                                                                   ccGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCAt                                                                                                                                                 <  1:341884/251‑1 (MQ=255)                                                                                    cgacgaTACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATc                                                                                                                                                <  1:184370/251‑1 (MQ=255)                                                                                     gacgaTTCGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGGTCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATcc                                                                                                                                               >  1:121331/1‑251 (MQ=255)                                                                                       cgaTACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCGAAAACTGGGCTGGCACCATCCGc                                                                                                                                             >  1:304745/1‑251 (MQ=255)                                                                                                  aaCGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCGGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATc                                                                                                                                  <  2:362594/251‑1 (MQ=255)                                                                                                    cGTTTTGAATCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCTGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATccc                                                                                                                                <  1:79210/251‑1 (MQ=255)                                                                                                     gTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAACGCGTGAAAGATAAGCCGTCGCTGGATATCTGCCGTACGGTGATTGGCTTTGGTTCGCCACATAAAGCATGTAAGGAAGGAGCGCAAGGCGCACCCCTGGGGGGAGAAGAAGGGGCGCTGGCACGGCAAAAACATGGCGGGCACCGTCCACCATGTGAGATcccc                                                                                                                               >  2:5656/1‑250 (MQ=255)                                                                                                          gAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAag                                                                                                                           <  2:371104/250‑1 (MQ=255)                                                                                                           atgCCTAGCACGGGCAGGTGATCCATGAAATCGACGTTCACGATCCGCAGGCGGTGAAGGAAGCGATCCGTGAAGCTCAAAGCGTGAAAGAGAAGCCGTCGCTGATTATTTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCTCGGCGCACCACGGCGGGAAGAAGAAGTGGCGCTGTCACGGCAAAAACTGGTCTGGCACCGTCCGCCATTTGAGATTCCTAAagag                                                                                                                         <  2:280138/249‑1 (MQ=255)                                                                                                             gCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTAGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAag                                                                                                                           >  1:430126/1‑247 (MQ=255)                                                                                                                                   aTGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGt                                                                                                 <  1:88827/251‑1 (MQ=255)                                                                                                                                    tGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGGTTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGGGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTg                                                                                                >  1:95098/1‑251 (MQ=255)                                                                                                                                                           aGGCGGTGCAGGAAGCGATACTTGAAGCGCAAAGCGTCAAAGATAAGCCATCGCTGATTACCTGCCGTACGGTGATTGTCTTTGGTTCGCCCAATAAAGAAGGTCAGGAAGAGGCGCACGGCGCACCACTGGGCGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCTTTTGAGATCCCTCAAGAGATTTAACACGCCTGGGAGGCCCGTGAAAAAGGCGAAAAAgagccgc                                                                           >  2:294363/1‑243 (MQ=255)                                                                                                                                                                  gaaggaagCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTgg                                                                   <  1:403852/250‑1 (MQ=255)                                                                                                                                                                  gaaggaagCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTgg                                                                   <  1:34992/250‑1 (MQ=255)                                                                                                                                                                     ggacgctaTCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCTCTGATTATCTGCCGTACGATGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATGTGAGATCCCTAAAGAGATTTATCACGCCTGGTATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATg                                                               <  2:185065/244‑1 (MQ=255)                                                                                                                                                                       aagCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATgag                                                             <  2:121331/251‑1 (MQ=255)                                                                                                                                                                              ccTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGttt                                                       <  1:323148/250‑1 (MQ=255)                                                                                                                                                                                   aaGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTgccgcc                                                 >  1:399239/1‑251 (MQ=255)                                                                                                                                                                                                    gATAAGCCGTCGCTGATTATCTGCCTTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATcc                                <  1:197456/251‑1 (MQ=255)                                                                                                                                                                                                        aGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCTTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCaa                            <  1:383796/251‑1 (MQ=255)                                                                                                                                                                                                           cGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCAAATAAAGCAGGCAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTg                         >  1:430680/1‑251 (MQ=255)                                                                                                                                                                                                                 tGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCACATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAACAGGCTCATCCGCCACTGGCaga                    >  1:250405/1‑250 (MQ=255)                                                                                                                                                                                                                               gcGGTGATGGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGAGCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCGGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACg     <  2:52422/250‑1 (MQ=255)                                                                                                                                                                                                                                  gTGATTGGCTTTGGTTCGCCAAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACGGAt  <  2:349079/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       GGCACGCTGGGACTGGGCAAGCTGATTGGTTTTTACGATCACAACGGTATTTCCATCGACGGTGAAACAGAAGGCTGGTTTACCGACGATACGGCAAAACGTTTTGAAGCCTATCACTGGCATGTGATCCATGAAATCGACGGTCACGATCCGCAGGCGGTGAAGGAAGCGATCCTTGAAGCGCAAAGCGTGAAAGATAAGCCGTCGCTGATTATCTGCCGTACGGTGATTGGCTTTGGTTCGCCGAATAAAGCAGGTAAGGAAGAGGCGCACGGCGCACCACTGGGGGAAGAAGAAGTGGCGCTGGCACGGCAAAAACTGGGCTGGCACCATCCGCCATTTGAGATCCCTAAAGAGATTTATCACGCCTGGGATGCCCGTGAAAAAGGCGAAAAAGCGCAGCAGAGCTGGAATGAGAAGTTTGCCGCCTATAAAAAGGCTCATCCGCAACTGGCAGAAGAGTTTACCCGACGGAT  >  NC_000913/2580143‑2580618

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