Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	83,579	1	.	A	78.6%	101.0 / 22.0	42	intergenic (‑50/+43)	leuA/leuL	2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (6/3);  new base A (14/19);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑A‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCG  >  NC_000913/83336‑83816                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              aCCTGGCGGGCGATGGGTTGCACCGAGTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACGCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGGTTGGGGCCTTGTCTCTTTTAGAGCGCCTCGCGTCGGGCAGA‑A‑A‑A‑A‑AA‑A‑CCcgcg                                                                                                                                                                                                                                         <  2:375569/251‑1 (MQ=255)       cgggcgATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTGTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATgg                                                                                                                                                                                                                                   <  1:232398/251‑1 (MQ=255)         ggcgATGGTTTGCACCGAGTCAAAATCGCCCGTGGAAGAGACGGGGAAACCGACTTCCATCACTTCAACACCCATACGCGCAAGGGCCAGCGCAATTTGCAGTTGTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCTAAAAGAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCGCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGgcg                                                                                                                                                                                                                                 <  2:216150/251‑1 (MQ=255)          gcgATGGATTGCACCGATTCAAAATGGCCCGGCGAAGAGTCGGGGAAACCGACTTCCATCACGTCAACACCAATACGCTCAAGCGCCAGCGCAATTTGCAGTTTTTCTATCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGACACTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑ACA‑AC‑C‑CCCACGCAATGgcgc                                                                                                                                                                                                                                >  2:417694/1‑251 (MQ=255)            gATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTTTCTCGTTTAGAGCGCCGCGCTTCTGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCgg                                                                                                                                                                                                                              <  2:571264/251‑1 (MQ=255)               ggTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGtt                                                                                                                                                                                                                           >  2:48852/1‑251 (MQ=255)               ggTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGtt                                                                                                                                                                                                                           >  2:307340/1‑251 (MQ=255)                         gATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACGTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGa                                                                                                                                                                                                                 <  2:155127/251‑1 (MQ=255)                          aTTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCGCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGAc                                                                                                                                                                                                                <  1:49984/251‑1 (MQ=255)                                 aTCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGc                                                                                                                                                                                                         <  1:66346/251‑1 (MQ=255)                                    cccgatctaaGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCGCAAGGGCCAGCGCAAGTTGCAGTGTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAAGGGCTTGCTGGCTCATGGTTTGGGTCCTGGTCTCGTTTAGAGCGCCCCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCTCGCAATGGCTCGGGTTTTTTGTTTGACTGCGTGCTgg                                                                                                                                                                                                      <  1:235770/243‑1 (MQ=255)                                     cccGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGGAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTGAGAGCGCCTCGCTACGGGCA‑T‑A‑A‑A‑A‑AAGC‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGc                                                                                                                                                                                                     <  2:532485/251‑1 (MQ=255)                                            aaGAGACGGGGAAACCGCCTTCCATCCCGTCAACACCCCTACGCTCAAGGGCCAGCGCCATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGCAACTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑ACAA‑A‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGCCTGCGTGCTgg                                                                                                                                                                                                      >  2:235770/1‑243 (MQ=255)                                                   ggggAAACCGACTTACATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGc                                                                                                                                                                                       >  1:463465/1‑251 (MQ=255)                                                    gggAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGGGCTGGCTTAATGCTGGATGcc                                                                                                                                                                                      >  2:329848/1‑251 (MQ=255)                                                                              acCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaaga                                                                                                                                                              >  1:169798/1‑249 (MQ=255)                                                                                  aTACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGc                                                                                                                                                        <  2:426301/251‑1 (MQ=255)                                                                                  aTACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTAACACTAACGCTCGACTGTAACCACTGTTCACCGCCGCGCAATGTGGAATCGAAAATAATGACTTGCTGGATCATGGTATGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGCCC‑T‑A‑A‑A‑A‑AA‑AACCAGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTACTGCTGGATGCCGCTCACTCGTCTCCCGCGCAACGACGATGc                                                                                                                                                        >  1:567866/1‑251 (MQ=255)                                                                                      gCTCAAGGGCCAGCGAAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                                                                                                    >  2:105301/1‑251 (MQ=255)                                                                                         cAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt                                                                                                                                                 <  1:48852/251‑1 (MQ=255)                                                                                                   cgcAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc                                                                                                                                       >  2:402333/1‑251 (MQ=255)                                                                                                     cAATTTGCAGTTTGTCTTTCACACTCAAGCTTGCCTGTAACGCCGGTTCACCGGCGCGCAATGTTGTATCGAAAATAATGACTTGCTAGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGTGCGGGTTTTTTGTTTGACTGCGTGCGGGCTTAATGGTGGATTCCGCACACTCGTCTAGCTCGCGAAGAAGAGGCGTTGAGTAGTAGTAGGCCg                                                                                                                                     <  2:243042/251‑1 (MQ=255)                                                                                                                 tttctttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAGGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                                                                                                          <  1:307340/250‑1 (MQ=255)                                                                                                                     tttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgt                                                                                                                     >  1:125802/1‑251 (MQ=255)                                                                                                                       tcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgtg                                                                                                                    <  2:68008/250‑1 (MQ=255)                                                                                                                        cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGa                                                                                                                   <  1:402333/250‑1 (MQ=255)                                                                                                                         ccacTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTACTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGt                                                                                                                 <  2:463465/250‑1 (MQ=255)                                                                                                                              cAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCGCTTTTAGAGCGCCTCGCTGCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTa                                                                                                            <  2:529699/251‑1 (MQ=255)                                                                                                                                       cTGTAACGCCTGTTCACCGTAGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTc                                                                                                   >  1:227892/1‑251 (MQ=255)                                                                                                                                               ccTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa                                                                                           <  1:80058/251‑1 (MQ=255)                                                                                                                                                     cACCGTCGCGCAATGTGGTATCGAAAATAATGACGTGCTGGCTCATGGTTTGGGTCCTTGTCTCGTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGat                                                                                     <  2:510623/251‑1 (MQ=255)                                                                                                                                                           cgcgcggtgtgGTATCGAAAATAATGATTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑TCA‑A‑A‑A‑AA‑A‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGAATGCGTGCTGGGTTGATGCTGGTTGCCGCTCACTCGTCTACCGCGCCAAGTAGACGCGTTTAGTGGTAGTAGGCCGATAAAGCGACCGATGTGAGTCATTAAATCAGCTCCAGATGATTGCGCTGTGCtt                                                                               >  2:85412/8‑251 (MQ=255)                                                                                                                                                             cgcAATGTGGTAGCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCtttt                                                                             <  1:180597/251‑1 (MQ=255)                                                                                                                                                                 aTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑ACA‑A‑AA‑A‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGCTGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTag                                                                           >  1:541636/1‑249 (MQ=255)                                                                                                                                                                      gtagcGAAAATAATGACTTGCTGGCTCATTGGTTGGGGCCTTGTCTCTTTGAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑A‑CCCGCGCAATGGCGCGGGTTTTTTGTTGGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGAGAAAGCGACCGATGTTAGTCAGTAAATCAGCTCCAGACGAATGCGATATGCTTTTAGAGGTACtc                                                                   <  2:348098/247‑2 (MQ=255)                                                                                                                                                                        aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTgg                                                                  >  2:110605/1‑251 (MQ=255)                                                                                                                                                                                            ggCTCATGGTTTGGGACCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGACTACGATATGCTTTTAGGGTTAATGGATACAAAAACGGATGTCacc                                              >  2:9383/1‑249 (MQ=255)                                                                                                                                                                                             gCTCATGGTTTGGGTCCTTGGCTCTTTTAGAGCGCCTAGAGCCGGGCA‑TCA‑A‑A‑C‑AA‑A‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCCAAGAAGATGCGTTTAGTAGTAGTCGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAcc                                             >  1:37165/1‑251 (MQ=255)                                                                                                                                                                                                    gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGc                                      >  1:140201/1‑251 (MQ=255)                                                                                                                                                                                                       tggttCCTTGTCTCTTTTAGAGCGGCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCCGCTTCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCaa                                    <  1:417694/246‑1 (MQ=255)                                                                                                                                                                                                                                     ttCGGGCA‑T‑A‑A‑A‑A‑AA‑AACCCGCGCAATGGCGCGGGTTTTTTGATTGACTACGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCAt     >  1:358902/1‑251 (MQ=255)                                                                                                                                                                                                                                        gggCA‑T‑ACA‑A‑A‑AA‑A‑CCCGCGCCATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAACACGTCCCGCCAGCGTGAGTTCTGCATCCg  >  1:497338/1‑251 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              ACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCA‑T‑A‑A‑A‑A‑AA‑A‑CCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCG  >  NC_000913/83336‑83816

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