Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,265,668	C→T	T452I (ACC→ATC)	dnaB →	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,265,668	0	C	T	98.2%	169.1 / NA	54	T452I (ACC→ATC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/28);  minor base G (0/1);  total (25/29)
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.75e-01

ACTGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACG  >  NC_000913/4265423‑4265902                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           aCTGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCttt                                                                                                                                                                                                                                        >  2:97437/1‑250 (MQ=255)     gCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACgg                                                                                                                                                                                                                                   <  2:247461/251‑1 (MQ=255)     gCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTGTAACgg                                                                                                                                                                                                                                   <  2:25028/251‑1 (MQ=255)          agaacagcaCGTGCCGGTGGAGGCGATGTCCCAGTTGAACCGTTCTCTGGAACAACGGGCCGACAAACGCCAGGTCCACATCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCAAGGTTATCTATCGTGATGAGGAGTATCAAGAAAACAGTGATTTAAAAGGCATCGAGGAAATTATTATCGGTCCACACCGTAACGGCCCAATCAGGACGGTACGCATGATCATTAACGGTCAAt                                                                                                                                                                                                                              >  2:587678/9‑251 (MQ=255)                 aaCGTGCCGGTTGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGCCAACTCCGTCCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTTATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTGAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTcgcg                                                                                                                                                                                                                       <  1:393688/251‑1 (MQ=255)                               ggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTAt                                                                                                                                                                                                          <  1:354125/250‑1 (MQ=255)                                    tGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAACTTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgggg                                                                                                                                                                                                    >  1:97392/1‑251 (MQ=255)                                           gTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGCTCTTTAACGGTCAATGGTCGCGCTTCGCCAACTATGCGGGGCCGCAGt                                                                                                                                                                                             >  2:365453/1‑251 (MQ=255)                                            ttGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGt                                                                                                                                                                                             >  2:154517/1‑250 (MQ=255)                                                        ctGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaa                                                                                                                                                                                  <  2:85838/249‑1 (MQ=255)                                                        ctGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaa                                                                                                                                                                                  <  2:259576/249‑1 (MQ=255)                                                        ctGGAAAAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaa                                                                                                                                                                                  <  2:17055/249‑1 (MQ=255)                                                           gAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTTCGTCGACGaataata                                                                                                                                                                             <  1:510760/251‑1 (MQ=255)                                                            aacaacGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATCATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                                                                                                                                                             >  1:189405/1‑250 (MQ=255)                                                             acatcgggacGACAAACGCCCGGTCAACTCCGATCTGCGTGAATCTGTCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACCCCTGATCTTTAACGGTCAATGGTCTCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataataat                                                                                                                                                                           <  1:391948/242‑1 (MQ=255)                                                               aacGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCTGTAAACAACGTAACGGCCCAATCGGGATGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACTAATAATAatta                                                                                                                                                                         <  1:74188/251‑1 (MQ=255)                                                                       gACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGa                                                                                                                                                                 <  1:32459/251‑1 (MQ=255)                                                                        aCAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGa                                                                                                                                                                 >  2:215801/1‑250 (MQ=255)                                                                                  ggTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAAtt                                                                                                                                                      >  1:343025/1‑251 (MQ=255)                                                                                   gTCAACTCCGACCTGCGTGAAACTGGCTCTAGCGAGCAGGATGCGGACTTGATCATGTATATCTATCGTGATGAGGTGTATCACGAAAACAGTGATGTAAAAGGCATCGAGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCATGAACGGGCAATGGTCGCGCTTCGACAACTAGGCTGGGACGCAGTACGACGACGAAGAATAAGTATTGGATAAATTAGGTAATTa                                                                                                                                                     <  2:278654/251‑1 (MQ=255)                                                                                              ccTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACt                                                                                                                                          >  1:112155/1‑251 (MQ=255)                                                                                                 gCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTa                                                                                                                                        <  1:287251/250‑1 (MQ=255)                                                                                                 gCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATGTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGTGACGGTACGCCTGATCTTTAACGGGCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTAGGAATTAGGTGAGTAAAGGAAACTCTTAt                                                                                                                                       <  2:554241/251‑1 (MQ=255)                                                                                                                tATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCTGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATg                                                                                                                        >  1:457107/1‑251 (MQ=255)                                                                                                                 aTCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAAt                                                                                                                         <  1:154517/249‑1 (MQ=255)                                                                                                                 aTCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATg                                                                                                                        <  1:132831/250‑1 (MQ=255)                                                                                                                 aTCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGAATTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATg                                                                                                                        <  1:9464/250‑1 (MQ=255)                                                                                                                   cGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCGTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCaa                                                                                                                     <  1:278920/251‑1 (MQ=255)                                                                                                                             gCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACt                                                                                                            <  2:261605/250‑1 (MQ=255)                                                                                                                             gCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTAGTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtg                                                                                                           <  2:549727/251‑1 (MQ=255)                                                                                                                               ggACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtgtt                                                                                                         >  1:497718/1‑251 (MQ=255)                                                                                                                                    tGATCCTGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGAt                                                                                                    >  2:290394/1‑251 (MQ=255)                                                                                                                                    tGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAAGCGGGACGGTACTCCTGAGCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCTGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGa                                                                                                     <  2:345515/250‑1 (MQ=255)                                                                                                                                              tatctatcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAAAGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAACAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACcgccgc                                                                                          >  1:401756/1‑251 (MQ=255)                                                                                                                                                   atcGTTATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGTACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGctct                                                                                     <  1:574646/251‑1 (MQ=255)                                                                                                                                                   atcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCcgcg                                                                                         <  1:577363/247‑1 (MQ=255)                                                                                                                                                        gatgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCACTGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGac                                                                                >  1:527127/1‑251 (MQ=255)                                                                                                                                                         atgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGac                                                                                >  1:220862/1‑250 (MQ=255)                                                                                                                                                         atgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAACCGGGACGGTACCCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTAATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATGAACCGCCGCGCTCTGCGaca                                                                               <  2:497718/251‑1 (MQ=255)                                                                                                                                                             ggTGTATCACGAAAACAGGGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGGCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTAGGAATTAGTTAATTAAAGCAATCACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAAc                                                                           <  1:567120/251‑1 (MQ=255)                                                                                                                                                                    cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAAc                                                                    <  2:400282/251‑1 (MQ=255)                                                                                                                                                                        aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCt                                                                >  1:411212/1‑251 (MQ=255)                                                                                                                                                                                gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAAc                                                        >  1:229272/1‑251 (MQ=255)                                                                                                                                                                                     aaaaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAAc                                                                                                             <  2:506538/193‑1 (MQ=255)                                                                                                                                                                                     aaaaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAAc                                                                                                             >  1:506538/1‑193 (MQ=255)                                                                                                                                                                                                cgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGCCGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTa                                         >  2:405393/1‑250 (MQ=255)                                                                                                                                                                                                   aaaTTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAAt                                     >  2:507325/1‑251 (MQ=255)                                                                                                                                                                                                   aaaTTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTAAGCCTGATCTTTAACGGTAAATAGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATGTTATGAATTAGGGAATTAAAGAAAACACTTATCAAGGAACACAAATGCAAGCGGCACCTGTTGTGATTAACCGCCGCGCTATGCGACACAACCTCCAACGTCATCGTGAACTGGCCCCGGCCAGTAAAAt                                     >  2:135022/1‑251 (MQ=255)                                                                                                                                                                                                    aaTTATTATCGGTAAACAACGTAACGGCCCAATAGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATg                                    >  2:419811/1‑251 (MQ=255)                                                                                                                                                                                                     ataattatCGGAAAACAACGTAACGGCCCAATCGGGACGGTACGCCGGATCTTTAACGGTAAAGGGTCGCGCTTCGAAAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTAATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATgg                                   <  2:278589/248‑1 (MQ=255)                                                                                                                                                                                                                        gTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGa                 >  1:278408/1‑250 (MQ=255)                                                                                                                                                                                                                            cGGCCCAATCGGGAAGGTAGGCCTGATCTTTAACGGCCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGCATAATAATTATTTTATGACTTCGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACAGCCGCGCTCTGCGCCACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGAGGCGGTGAAAGCGAActct            >  2:364989/1‑248 (MQ=255)                                                                                                                                                                                                                                      gggACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACg  >  1:75256/1‑251 (MQ=255)                                                                                                                                                                                                                                      gggACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAAGGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACg  <  1:365453/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           ACTGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACG  >  NC_000913/4265423‑4265902

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