Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,474,822	0	A	G	21.9%	75.3 / 17.9	41	R48R (CGT→CGC)	tusB	sulfurtransferase complex subunit TusB
Reads supporting (aligned to +/- strand):  ref base A (18/14);  new base G (5/4);  total (23/18)
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.28e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCT  >  NC_000913/3474573‑3475069                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        cTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTAc                                                                                                                                                                                                                                                        >  1:225756/1‑251 (MQ=255)         tCGTAAACACGTAATAAACCGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGTTGATATACCTTAATGGGGGCATTACGCAGAc                                                                                                                                                                                                                                                  <  1:214462/249‑1 (MQ=255)               acacGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTGCGCAGACTTTCAAgg                                                                                                                                                                                                                                          <  1:128904/251‑1 (MQ=255)                cacGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGCTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGt                                                                                                                                                                                                                                         >  2:236110/1‑251 (MQ=255)                  cGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTGCGCAGACTTTCAAGGt                                                                                                                                                                                                                                         >  2:160851/1‑249 (MQ=255)                                  cacacaATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTACGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTg                                                                                                                                                                                                                       >  1:251889/1‑251 (MQ=255)                                  cacacaATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACt                                                                                                                                                                                                                        <  2:10868/250‑1 (MQ=255)                                    aacaATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGGCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTGCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCg                                                                                                                                                                                                                     <  2:251889/250‑1 (MQ=255)                                                                 gATGCCAAAAAGGTGGCAAGAAATATACACCGATCTCGCCATCCCCAGGCCATCTGTCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTGCGAAATTTGTCCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAg                                                                                                                                                                                        <  2:71107/251‑1 (MQ=255)                                                                  aTGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGt                                                                                                                                                                                       >  1:211599/1‑251 (MQ=255)                                                                        aaaaGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTGCGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGtct                                                                                                                                                                                 <  1:96164/251‑1 (MQ=255)                                                                               gtCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTGCGTCTCCTTCa                                                                                                                                                                           <  2:36161/250‑1 (MQ=255)                                                                                                      gCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAgcagcgcagc                                                                                                                                                    >  2:169612/1‑250 (MQ=255)                                                                                                       ccATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCgcagca                                                                                                                                                   >  2:138785/1‑250 (MQ=255)                                                                                                           cACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTGCGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATc                                                                                                                                              >  2:123037/1‑251 (MQ=255)                                                                                                             ccAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTGCGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGTAGCGCAGCAAAATCCg                                                                                                                                            <  1:50598/251‑1 (MQ=255)                                                                                                                                      aCCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGGTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAAtgt                                                                                                                   <  1:149777/251‑1 (MQ=255)                                                                                                                                         gTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAAtgtgtg                                                                                                                <  1:216038/251‑1 (MQ=255)                                                                                                                                          tAAGTCTGACGCAATCCGTATAGTCTATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACCCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCTGCCGCAGCAGCGCAGCACAATCCGTCCGCCAGGGTGCGCGATGACATGTGTgc                                                                                                               >  2:21692/1‑251 (MQ=255)                                                                                                                                                 gACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTgcagca                                                                                                           <  1:59707/248‑1 (MQ=255)                                                                                                                                                    gAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACg                                                                                                     <  1:80493/251‑1 (MQ=255)                                                                                                                                                       aTCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTGCGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGcct                                                                                                  >  2:112769/1‑251 (MQ=255)                                                                                                                                                       aTCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGCCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGcct                                                                                                  >  1:130948/1‑251 (MQ=255)                                                                                                                                                             ataGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGcctcaaaa                                                                                             <  1:57613/250‑1 (MQ=255)                                                                                                                                                                  cAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGcctcaaaacctca                                                                                        >  1:116311/1‑250 (MQ=255)                                                                                                                                                                       agaATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATa                                                                                   >  1:189447/1‑250 (MQ=255)                                                                                                                                                                          atgatgTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACAt                                                                               >  1:47785/1‑251 (MQ=255)                                                                                                                                                                           ggatgTCGTGAGAAATTTCACCAGCAAAACCGCGGGAAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCACAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATc                                                                              <  2:194519/250‑1 (MQ=255)                                                                                                                                                                              tgTCGTTCGAAATTTGACCAGTCAAACCGCGGGCACTAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTACTGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGt                                                                            >  2:270611/1‑250 (MQ=255)                                                                                                                                                                              tgTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTGCGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGt                                                                            >  1:21719/1‑250 (MQ=255)                                                                                                                                                                                  gTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTg                                                                       >  2:26965/1‑251 (MQ=255)                                                                                                                                                                                     cGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTGCGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGc                                                                     <  1:160851/250‑1 (MQ=255)                                                                                                                                                                                      gAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTAAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCg                                                                    >  1:123639/1‑250 (MQ=255)                                                                                                                                                                                           ttGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTc                                                               <  2:123720/250‑1 (MQ=255)                                                                                                                                                                                                          gcgGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCt                                               <  1:144078/251‑1 (MQ=255)                                                                                                                                                                                                                   aaGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAgc                                      <  2:30442/251‑1 (MQ=255)                                                                                                                                                                                                                           tgTTCACGGCATAGACGTTAATGGTGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGAGGCAGTTACGCCATCTTGCAATACTAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTTAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATATGCTTCGCGCGGCGTGGc                              >  2:116256/2‑251 (MQ=255)                                                                                                                                                                                                                            gTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCTAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTACGGCATCTGCTTCGAGCGGCGTGGCt                             >  2:258615/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGaca                       >  1:237293/1‑250 (MQ=255)                                                                                                                                                                                                                                             tAATGGGGGCATTGCGAAGACTTTCAAGGTAGCGGTTACCGTCACATGCGGCAGTTACGCCATCTTGCAATAATAGCAGATCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCACAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGgtgt            >  1:150331/1‑251 (MQ=255)                                                                                                                                                                                                                                                        ttACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCt  >  2:94879/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        CTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCT  >  NC_000913/3474573‑3475069

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