Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,208,132	0	G	A	64.8%	24.5 / 51.1	54	intergenic (‑600/‑15)	purH/rrsE	bifunctional AICAR transformylase/IMP cyclohydrolase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (10/9);  new base A (18/17);  total (28/26)
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.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGC  >  NC_000913/4207901‑4208360                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      cggttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTATGCATCAAACTTTTAAAtt                                                                                                                                                                                                                    >  2:466513/1‑251 (MQ=255)       ggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAagag                                                                                                                                                                                                              <  1:340496/251‑1 (MQ=38)              gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACTCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATc                                                                                                                                                                                                       >  1:116836/1‑251 (MQ=37)                 aaaaaTCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATg                                                                                                                                                                                                    >  1:408523/1‑251 (MQ=25)                  aaaaTCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGTGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATgg                                                                                                                                                                                                   <  2:532655/251‑1 (MQ=25)                                            tGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGcc                                                                                                                                                                         >  2:139603/1‑251 (MQ=25)                                                   ggAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAacac                                                                                                                                                                   >  1:311945/1‑250 (MQ=25)                                                                 aCGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAATATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACgg                                                                                                                                                   >  1:429897/1‑251 (MQ=25)                                                                     cACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAAc                                                                                                                                               <  1:139603/251‑1 (MQ=25)                                                                               cAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAgc                                                                                                                                     <  1:463114/251‑1 (MQ=37)                                                                                      cTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGAACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACCCGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGCACGCTGGTGGCAGGCCTAACACA‑TGCACGTCGAACGGTAACAGGCAGAAGCTTGcttc                                                                                                                              >  1:24542/1‑251 (MQ=11)                                                                                           gcgcGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGGCGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGc                                                                                                                         <  1:389774/251‑1 (MQ=14)                                                                                                  ggcaCTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAg                                                                                                                  <  1:169616/248‑1 (MQ=11)                                                                                                   gcaCTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGt                                                                                                                 <  1:190373/249‑1 (MQ=11)                                                                                                     aCTGCTCTTTCCCAATTTATCAGACAAACTGTGTGGGCACTCGAAGGTACGAATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGGCTCAAGAGTGAGCACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTGAATTGAAGCGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAAC‑CAGTGCAAGGCGGACGGGAACGGGAAGAAGCTTGCTTCTTTGcgtacgagggg                                                                                                               >  2:53848/1‑241 (MQ=11)                                                                                                       tGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGCTACGGATTCTTCACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGCTTGAACGCAGGCGGACGGACTAACCCA‑TGGAAGTCGACCGGTAACAGGCAGAAGcttccttcctggcttcagggcgcgg                                                                                                             >  2:240406/1‑234 (MQ=11)                                                                                                        gCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCgg                                                                                                            <  1:324159/251‑1 (MQ=31)                                                                                                        gCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTTAAGAGTTTGATCATGGTTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAATCTTGCTGCTTTGCTGACGAGTGGCg                                                                                                             <  2:290473/250‑1 (MQ=9)                                                                                                        gCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCgg                                                                                                            >  1:494543/1‑251 (MQ=14)                                                                                                           ctTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGa                                                                                                           <  2:164198/249‑1 (MQ=31)                                                                                                             ttAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACg                                                                                                         <  1:486743/249‑1 (MQ=31)                                                                                                             ttAACAATTTATCAGACAAACTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGATTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACTAGTGGCGGACg                                                                                                         <  2:542932/249‑1 (MQ=9)                                                                                                                   aTTTATCATACAATATGTATGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAATATGGATACCAAGTCTAAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAATCTTTTACATTGAGGAGTTTGATCATGGCTCAGATTGACCGCTGGCGGCAGGCCTCACCCA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCGGCCGAGTGGCGGaatggagggt                                                                                                  >  2:464659/1‑241 (MQ=11)                                                                                                                             cAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGAAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGGGTAATgggtgggg                                                                                       >  2:417368/1‑244 (MQ=9)                                                                                                                              aaTCTGTGTGTGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGATGGAACACGTAATCCCTTACGAAGTTTAATTCTTTGAGCATCATACTTTTAAATTGAAGAGATTGATCATGGCTATGATTGAACGCTGGCGGAAGGCCTAACACA‑TGCAAGTCGAACGGAAACAGGAAGCAGCTTGCTTCATTGCTGACGGGCGTCGGAAGGGAGAGTAATGTCTggggc                                                                                      >  2:79280/1‑249 (MQ=11)                                                                                                                                   gtgtgGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTg                                                                                  <  1:300332/250‑1 (MQ=14)                                                                                                                                       gggCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGGACGGTAACTGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGGGGAATGTCTGGGGAACTGCCTgg                                                                             >  2:236156/1‑250 (MQ=255)                                                                                                                                           aCTCGAAGATACGGATTCGTAACGTCGCAAGACGAAAAATGAATTCCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGa                                                                         <  2:303495/251‑1 (MQ=14)                                                                                                                                            cTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGGGGCGGACGGGTGGGGAATGGCTGGGAAACTGCCTggagggg                                                                        >  2:106503/1‑245 (MQ=255)                                                                                                                                            accgatcttACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGa                                                                         <  2:304178/242‑1 (MQ=14)                                                                                                                                                aaGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAggggg                                                                    <  2:494543/251‑1 (MQ=14)                                                                                                                                                   aTACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATa                                                                 <  1:167909/251‑1 (MQ=14)                                                                                                                                                    tACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGa                                                                         >  1:304178/1‑242 (MQ=14)                                                                                                                                                     ccGGAGTCTGAACGTCGGAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACATGTTATTCATGACGAAGGTTAATGCTTTAAGCACCAAACTTTCAAAGTGAAGAGTGTGATCATGGCCCAGATTGAACGCGTGCGGCAGGACTCCCTCA‑TGCAAGGCGAACGGTACCAGGAGGAAGCTTGCTTCTTCGCTGACGAGTGGCGGACGGCCGAGTAATGTCGCGGAAACCGCCTTATGGATGGTGATacc                                                               <  2:262857/250‑3 (MQ=2)                                                                                                                                                      cGGATGCTTAACGTCGCAAGACGAAAAATGAATTCCAAGTCTCAAGAGTGAACACGTAATTCAATACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACTGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCGGATGGAGTGGGATAacg                                                              <  2:343812/251‑2 (MQ=255)                                                                                                                                                       ggATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGTAGGGGGATAact                                                              >  1:143529/1‑250 (MQ=14)                                                                                                                                                             ttAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCTAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGa                                                        >  1:162759/1‑250 (MQ=14)                                                                                                                                                                 cGTCGCAATACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACgg                                                   <  2:429897/251‑1 (MQ=14)                                                                                                                                                                  gTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACTGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGt                                                  >  1:485653/1‑251 (MQ=14)                                                                                                                                                                     gCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCTTGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGCAAGAAGCTTGCTTCTTTGATGACGAGAGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTTCTGGATACGGGAGc                                               >  1:95198/1‑251 (MQ=11)                                                                                                                                                                           cGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGGGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATAc                                         >  1:284493/1‑251 (MQ=14)                                                                                                                                                                            gAAAAATGAATACCAAGTCTCAAGAGTGAACACGAAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGTGGGGGGTAACTACTGGAAACGGGAGCTACTAcc                                        >  1:164413/1‑251 (MQ=255)                                                                                                                                                                             aaaaaTGAGTACCAAGTCTCAGGAGTGAAAACGTAATTCATTACGAAGTTAAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGGGTAATGTATGGGAAACTGCCTGGTGGAGGGGGAGTACTAATGGGAACGGGAGCAAAAACCg                                       >  2:92788/1‑251 (MQ=11)                                                                                                                                                                                  tGAATACCAAGTCTCAAGCGTGAACACGTAATTCATTACGAAGTTTACTTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGTGGAATGTATGGGAAACTGCCTGATGGGGGGGGATAACTACTGGAAACGGGAGATTATACCGCATa                                   >  2:469718/1‑250 (MQ=11)                                                                                                                                                                                   gAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATaa                                  <  1:482265/250‑1 (MQ=14)                                                                                                                                                                                    aaTACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACg                                <  1:27063/251‑1 (MQ=31)                                                                                                                                                                                     aTACCAAGTCTCAAGCGTGAACACGTAATTCATTAAGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGt                               >  1:124680/1‑251 (MQ=14)                                                                                                                                                                                     aTACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGt                               <  2:408523/251‑1 (MQ=14)                                                                                                                                                                                            gTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGa                        <  1:153413/251‑1 (MQ=31)                                                                                                                                                                                                cAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGAcc                      <  1:310989/249‑1 (MQ=14)                                                                                                                                                                                                     gTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGTCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGGAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACAGGAGAGTAATGTCTGGGAAACTGCCTGGTGGAGGGGGATACATACAGGAAACGGTAGCTAATCACGACTAACGGCGCCAAACCAAAGAgg               >  2:159456/1‑251 (MQ=11)                                                                                                                                                                                                           acGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGAcc         >  1:530460/1‑251 (MQ=14)                                                                                                                                                                                                               aaTTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTc      <  1:478697/250‑1 (MQ=31)                                                                                                                                                                                                                  tCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAATCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGc  >  1:131694/1‑251 (MQ=14)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      CGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACA‑TGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGC  >  NC_000913/4207901‑4208360

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