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	97.5%	121.7 / NA	40	T452I (ACC→ATC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/19);  minor base G (0/1);  total (20/20)
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.68e-01

GCACTGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTT  >  NC_000913/4265421‑4265910                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   gCACTGGCGAAAGACCTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGAACGCCTGATCtt                                                                                                                                                                                                                                                 >  1:123091/1‑251 (MQ=255)   aCTGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCttt                                                                                                                                                                                                                                                >  1:13893/1‑250 (MQ=255)        cGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATTAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGt                                                                                                                                                                                                                                          <  1:367772/251‑1 (MQ=255)          aaaGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACTCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCCTCGCGGAAATTATTATCGTTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGAGCTTTAACGGTCa                                                                                                                                                                                                                                        <  2:216291/251‑1 (MQ=255)           aaGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTACGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGCCGGTACGCCTGATCTTTAACGGTCa                                                                                                                                                                                                                                        >  1:198799/1‑250 (MQ=255)            aGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAACCGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAAt                                                                                                                                                                                                                                      <  2:381644/251‑1 (MQ=255)                    aCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCgcgc                                                                                                                                                                                                                              >  1:363846/1‑251 (MQ=255)                       ggCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCATGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTTTGCCTGATCTTTAACGGTCAATGGTCGCGCTTc                                                                                                                                                                                                                           <  1:281496/250‑1 (MQ=255)                          cGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGAc                                                                                                                                                                                                                        >  1:218043/1‑251 (MQ=255)                                            aGTTGACCCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGACAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGAGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc                                                                                                                                                                                                        >  2:302338/1‑249 (MQ=255)                                               tGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTa                                                                                                                                                                                                    <  1:312915/250‑1 (MQ=255)                                                gAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACGTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCTGGGCCGCAGTacg                                                                                                                                                                                                  <  1:333345/251‑1 (MQ=255)                                                  aCCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCACTGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgac                                                                                                                                                                                                >  2:90556/1‑251 (MQ=255)                                                  aCCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgac                                                                                                                                                                                                <  1:90556/251‑1 (MQ=255)                                                                  acGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAattat                                                                                                                                                                                <  1:172803/251‑1 (MQ=255)                                                                    gTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGATTCGACAACTATGCGGGGCAGCAGTACGACGACGAATAATAATTAttt                                                                                                                                                                              <  2:249388/251‑1 (MQ=255)                                                                           cAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAt                                                                                                                                                                       <  1:162805/251‑1 (MQ=255)                                                                                           tCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGCCGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGACTTAGGTAATTAAAGCa                                                                                                                                                        >  1:219016/1‑250 (MQ=255)                                                                                               aCCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAacac                                                                                                                                                   <  2:109548/251‑1 (MQ=255)                                                                                                ccTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGATAATTATGCGGGTCCGCATTACGACGACGAATAATAATTATTTTATGTATTGGGTAATTAAAGCAAACACt                                                                                                                                                  <  2:47064/251‑1 (MQ=255)                                                                                                                   aTCGAGCAGGATGCGGACGTGAGCATGTTTATCTAGCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATTCGGGGCCGCAGTACGACGACGAATAATAAGTAGTTTATGAATTAGGTAATTAAAGCAAACGCTTATCAAGGAACACAAATg                                                                                                                                <  2:198799/250‑1 (MQ=255)                                                                                                                             aTGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAAc                                                                                                                     >  1:234585/1‑251 (MQ=255)                                                                                                                                           aTGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAAcc                                                                                                       <  2:155591/251‑1 (MQ=255)                                                                                                                                                     atcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTAGCGCTTCGACATCTATGCGGGTCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGGAATTAAAGCAAACACTTATCAAGGATCACAAATGCAAGCGGCATCTGTGGTGATTAACCGCCGCGctct                                                                                             <  2:234177/251‑1 (MQ=255)                                                                                                                                                       cGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGc                                                                                           >  1:316409/1‑251 (MQ=255)                                                                                                                                                       cGTGATGAGGTGTATCACGAAAACAGGGATTTAAAAGGCATCGCGGAAATTATTAGCGGAAAACAACTTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTg                                                                                            <  2:293095/250‑1 (MQ=255)                                                                                                                                                        gTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCg                                                                                          >  1:128828/1‑251 (MQ=255)                                                                                                                                                              aGGTGTATCACGAAAACAGGGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACCGTCACTGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTCCGACGCCGAATAATCATTATTTTCTGAATTAGGTCATTAAATCAAACACTTATCAAGTACCACTAATGCAAGCGGCAACTGTTGTGATTAACCTCCGCGCTCTACGacacc                                                                                     >  1:264/1‑249 (MQ=255)                                                                                                                                                                          aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAAAGTAAAGGCCCAATCGGGACGGTAAGCCTGATCTTTAACGGTCAATGGTCGCGATTCGACAACTATGCGGGGCCGCAGTAAGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCAGCGCGCTGCGACACAACCTGCAACGtat                                                                        >  1:238773/1‑249 (MQ=255)                                                                                                                                                                                 ggATTAAAAAGGCAGCGCGGAAATTATTACCGGTAAACAACGTAACGGCCCAAGCGGGACGGTTCGCCTGATCTTTAACGGCCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGCCGAATAAGAATTATTTTATGAATTAGGTAATTAAAGCAACCACTGATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGAGTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGaa                                                                 <  2:19639/250‑1 (MQ=255)                                                                                                                                                                                        aaaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGccc                                                          >  1:293662/1‑251 (MQ=255)                                                                                                                                                                                         aaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGcccc                                                         >  1:235357/1‑251 (MQ=255)                                                                                                                                                                                              aTCGCGAAAATTATTATCGGTAAACAACGTAACGGACCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACAATGCGGGGCCGCAGTACGACGACGAATAATCATTATTTTATGAATTAGGTAATTAAAGCAAAGACTTATCAAGGAACACAAATGCACGCGGCAACTGATGTGATTAACAGCCGCGCTCTGAGACACAACCTGCAACGTCTTCGAGAAATGGACCCTGcac                                                    >  2:159415/1‑249 (MQ=255)                                                                                                                                                                                                       attattATCGGAAAACAAAGTAACGGCCCAATCGGGAAGGTACGCCTGATCTTTAACGGTCAATGATCGCGCTTCGACAAATATGCGGGGACGCAGTACGACGACGAATAATAATTCTTTTATGAATTAGGTAATTAAAGCAAACCCTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACATGCAACGTATTCGTGCACTGGCCCCCGCCAGTAAAATgg                                           >  2:192907/1‑251 (MQ=255)                                                                                                                                                                                                          attatCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTg                                        >  2:125928/1‑251 (MQ=255)                                                                                                                                                                                                                                      tCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTAt                 >  1:285514/1‑246 (MQ=255)                                                                                                                                                                                                                                      tCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTc             >  1:143248/1‑250 (MQ=255)                                                                                                                                                                                                                                        gggACGGTACGCCTGATCTTTACCGTGCAATGTGCGCGCTTCGACAACTATGCGGGGCCGCATTACGACGAAGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAAGGCAAGCGGCAACTTTTGGGATTATCCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGGCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACg          <  1:265048/251‑1 (MQ=255)                                                                                                                                                                                                                                              gtccGCCTGATCTTTAACGGCCAATGGTCGCGCTTCGAAAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGtcttc    <  2:316409/248‑1 (MQ=255)                                                                                                                                                                                                                                                aCGCCTGATCTTTAACGGGCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTACCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTcttctt  <  1:80595/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   GCACTGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTT  >  NC_000913/4265421‑4265910

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