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.7%	244.4 / NA	74	T452I (ACC→ATC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (39/34);  minor base G (0/1);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGT  >  NC_000913/4265430‑4265904                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              aaaGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTc                                                                                                                                                                                                                                    >  1:271671/1‑250 (MQ=255)         gAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTc                                                                                                                                                                                                                             <  1:134878/249‑1 (MQ=255)         gAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATGTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGGAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTcgc                                                                                                                                                                                                                           <  2:41255/251‑1 (MQ=255)          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTcgcg                                                                                                                                                                                                                          >  1:544954/1‑251 (MQ=255)             gTGCCGGAGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTACCGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCtt                                                                                                                                                                                                                       >  1:573486/1‑251 (MQ=255)                    gggtggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAAc                                                                                                                                                                                                                <  2:97724/250‑1 (MQ=255)                        ggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGAGCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCAGCGCGTAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAAGGGTCGCGCTTCGACAACTATg                                                                                                                                                                                                            <  2:279748/251‑1 (MQ=255)                        ggCGCTGTCCCAGGTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTAGCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTacg                                                                                                                                                                                                            <  1:31182/251‑3 (MQ=255)                          cgcTGTCCCCGTTGAACCGTTCTCTGGAACAACGTGACGACAACCGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTTTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCg                                                                                                                                                                                                          >  2:343953/1‑251 (MQ=255)                           gcTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgg                                                                                                                                                                                                         <  1:267769/251‑1 (MQ=255)                              gTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGCTCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGc                                                                                                                                                                                                      >  2:41605/1‑251 (MQ=255)                              gTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCggg                                                                                                                                                                                                        <  1:300464/249‑1 (MQ=255)                              gTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCGCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACATCGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCggg                                                                                                                                                                                                        <  1:385011/249‑1 (MQ=255)                                 cccgTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGCACGCCTGATCTTTAACGGTCAATGGGCGCACTTCGAAACCTATGCGGGGcagc                                                                                                                                                                                                   >  2:315824/4‑248 (MQ=255)                                   aGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc                                                                                                                                                                                                   >  1:407864/1‑249 (MQ=255)                                       gAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTac                                                                                                                                                                                              <  1:508674/250‑1 (MQ=255)                                                     aacaacGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataataa                                                                                                                                                                               <  1:489407/251‑1 (MQ=255)                                                     aacaacGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                                                                                                                                                                <  2:426912/250‑1 (MQ=255)                                                     aacaacGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                                                                                                                                                                >  1:426912/1‑250 (MQ=255)                                                           gTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGTCTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACGATGCGGGGCCGCAGTACGACGACGAATAATAagtattg                                                                                                                                                                         <  2:233143/251‑2 (MQ=255)                                                                gACAAACGCCCGGTCAACTCCGACCTTCGTGAATCTGGCTCGATCGAGCAGTGTGCGGAATTGAGCATGTTTAGCTATCGTGGTGAGTTGT‑ATCACGAAAACAGTGATTTAAAAGGTATCGCGGAAATTATTATCGGTAAACGACGTAGCGGCC‑CAATCGGGACGGTACGCTTGAGCTTTGACGGTCAATGGTCGCGCTTCGGCTACTATGCGGGGCCGCAGTACGACGACGAATAATAAGTATTTTATGa                                                                                                                                                                    <  1:56987/251‑1 (MQ=255)                                                                gACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATg                                                                                                                                                                     >  2:17185/1‑250 (MQ=255)                                                                  cAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATATAGCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTAACGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAt                                                                                                                                                                  <  1:11493/251‑1 (MQ=255)                                                                   aaaCGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAtt                                                                                                                                                                 >  1:334957/1‑251 (MQ=255)                                                                         ccGGTCAACTCCGCCCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTaa                                                                                                                                                          <  1:464603/251‑1 (MQ=255)                                                                          cGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAAt                                                                                                                                                         >  1:82860/1‑251 (MQ=255)                                                                           ggTCAACTCCGACCTGCTTGAATATTGCTCTATAGAGCAGGATGCGGACTTGATCATGTTTAACTATCGTGATGAG‑TGTGATCACGAAAACAGTTATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACTGTACGCCTGATCTTTAACGGGCAATGTTCGCGCTTTGACATCGATTCGGGGCCGCAGTACGACGAAGAATTATAattagttgaggatttagttaagt                                                                                                                                                         <  2:380922/251‑20 (MQ=255)                                                                               aaCTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTc                                                                                                                                                                                                                             <  1:583496/179‑1 (MQ=255)                                                                               aaCTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGAGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTc                                                                                                                                                                                                                             >  2:583496/1‑179 (MQ=255)                                                                                     gACCTGCGTGCATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATACTAATTATTTTATGAATT‑AGGTAATTAAAGCAAac                                                                                                                                               >  2:364211/1‑250 (MQ=255)                                                                                     gACCTGCGTGACTCTAGCTCTATCGAGCAGGATGCGGCCTTGACCATGTTTATCTATCGTGATGAGGTGA‑ATCACGACACAAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGGAACGGCC‑CACTCGGGACGGTACGCCTGATCTTTAACGGTCAATTGCCGCGCTTCGACAACTATGCGGGGCTGCAGTACGACGCCGAAGAATAATAATTTTATGAATT‑AGGCAATTAAAGCAAac                                                                                                                                               >  2:583738/1‑250 (MQ=255)                                                                                         tGCGTGAATCTCGCTCTATAGCGCAGGATGCGGACTTGATCATGTTGATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCATGATCTTTAACGGTCAATGGTCGCGCATCGACAACTATGAGGGGCCGCAGTACGTCGACGAATATTATTAATTTAAGGAATT‑‑GGTAGTTAAAGCAAACAATTac                                                                                                                                         >  2:21345/1‑250 (MQ=255)                                                                                             tGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGCAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAGCTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTATATGAATT‑AGGTAATTAAAGCAAACACTTATCa                                                                                                                                       >  2:136093/1‑250 (MQ=255)                                                                                                    gggtatATCGATCAGGATGCGGACTTGATCATGTGTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTAACGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAacac                                                                                                                               <  1:292292/246‑1 (MQ=255)                                                                                                    ggctcgttCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCACTTCGACATCTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAATGCAAACACTTATCaattaa                                                                                                                                  <  2:352878/241‑5 (MQ=255)                                                                                                    ggCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATGATTAGCGGTAAACAACGTAACGGCC‑CAAGCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCTACAACTATGCGGGGCTGCAGGACTACGTCGAATAATAATTATTTTAGGAATT‑AGTTAATTAAAGCAAGCACTTAGCAAGGaa                                                                                                                                  <  2:384598/248‑1 (MQ=255)                                                                                                            cGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCaa                                                                                                                       <  2:271671/251‑1 (MQ=255)                                                                                                            cGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCa                                                                                                                        <  1:36728/250‑1 (MQ=255)                                                                                                               ggagtaTGCGGAATTGATCATGTTTAGCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTGAATGGAACTGCGGAAACTATTAGCGTTAAACAACGAAACGGCCACAA‑CTGTACGGTTCGCGTGTTCGTGAACGGTCAATGGTAGCGCTTCGAGAACTATGCAGGGACACAGTACGACGACCAATAATAATTATTTTATGAATG‑AGGTAATTTAAGCAAACGCTTATCaatgatcaacaatgaaagcg                                                                                                                    <  2:27977/246‑19 (MQ=255)                                                                                                                aaGGATGCGGACTTGATCCTGTTTATCTATCGTGATGAGGCGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGCCATTATTATCGGTAAACAACGTAACGGCA‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATACTTATTTTATGAATT‑AGGTAATAAAAGCAAACCCTTATCAAGGAACACAAAAGAAAGCgg                                                                                                                   >  2:480077/2‑251 (MQ=255)                                                                                                                      gCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAGCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtg                                                                                                             <  1:126403/251‑1 (MQ=255)                                                                                                                         gACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGAACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtgttt                                                                                                          >  1:418985/1‑250 (MQ=255)                                                                                                                            ttGATCATGCTTATCTATCGTGATGAGGTGT‑CTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGAACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTtgtg                                                                                                        >  1:128096/1‑250 (MQ=255)                                                                                                                                   tGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAAccg                                                                                                >  1:166166/1‑251 (MQ=255)                                                                                                                                                tgatgaGGTGT‑ATCACGCAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCg                                                                                    >  2:34823/1‑250 (MQ=255)                                                                                                                                                     aGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACTCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATg                                                                                                                                                                     <  1:276403/165‑1 (MQ=255)                                                                                                                                                     aGGTGT‑ATCACGAAAAAAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATg                                                                                                                                                                     >  2:276403/1‑165 (MQ=255)                                                                                                                                                      ggTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAAc                                                                             >  1:28091/1‑251 (MQ=255)                                                                                                                                                          t‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGc                                                                         >  1:493046/1‑251 (MQ=255)                                                                                                                                                            aTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGCACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTCTGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTCTGACTT‑CGGTAATTAAAGCAAACACTTATCCAGGAACACAAATGCAAACGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGc                                                                         >  1:572829/1‑250 (MQ=255)                                                                                                                                                                cGAAACCAGTGATTAAAAAGGAATCGCGGAAATAATTATCGGTAAACAACGTAACGGCC‑CAACCGGGACGGTACGCCTGATCTTGAACGGTCAATGGTCGCGTTTCGACAACTATGCGGGGCGGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATAAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCGGCAACGt                                                                    <  2:200809/251‑1 (MQ=255)                                                                                                                                                                  aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCt                                                                  >  2:471745/1‑251 (MQ=255)                                                                                                                                                                    aaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACACCGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGAAAACTATGCGGAGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATACAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTc                                                                >  1:354568/1‑251 (MQ=255)                                                                                                                                                                          gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGaa                                                           <  1:364211/250‑1 (MQ=255)                                                                                                                                                                          gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAAc                                                          >  2:379262/1‑251 (MQ=255)                                                                                                                                                                           aTTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACt                                                         >  2:235103/1‑251 (MQ=255)                                                                                                                                                                              tAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTgg                                                       <  1:584993/250‑1 (MQ=255)                                                                                                                                                                                aaaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CTATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGAAGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGccc                                                    >  1:575082/1‑251 (MQ=255)                                                                                                                                                                                 acggCATCGCGGCAATTATTCTCGGTAACCAACGTAACGGCA‑CAATCGGGACGGTACGCCTGTTCTTTAACGGTCCCTGGTCGCGCTTCGACAACTATACGGGGCCGCAGTACGACGAAGACTAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACCAATGCAAGCGGAAACTGTTGGGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGTcccc                                                   >  2:365172/3‑251 (MQ=255)                                                                                                                                                                                   ggCATCGCGGAAATGATTATCGGTAAACAACGTAACGGCC‑CAGTCGGGACGGTACGCTTGATCTTTAACGGGCAATGGTCGCGCTTCGAGAACTATGCGGGTCCGCAGTACGACGACGAATAATAAGTATTTTATGAATT‑AGGTAGTTAAAGCAAACACTTATTAATGAACACAAATGCAAGCGGCAACTGTTGTGATGGACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTg                                                 <  1:573771/251‑1 (MQ=255)                                                                                                                                                                                   ggCATCGCGGAAATAATTATCGGTAAACACCGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAAGGGTCGCGCTTCGACAACGATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTg                                                 <  2:470539/251‑1 (MQ=255)                                                                                                                                                                                        cgcgGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGt                                            >  1:176360/1‑251 (MQ=255)                                                                                                                                                                                            gAAATTATTATTGGTAAACAACGTAACGGCC‑TCATCGGGACGGTACACCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGCCGAATAATAATTATTTTATGATTT‑AGGTACTTAAAGCAAACACTTATCAAGGCACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTTCAACGTCTTCGTGAACTGGCCCCTGCCCGTaaa                                         >  2:71172/1‑250 (MQ=255)                                                                                                                                                                                               attattATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATgg                                     >  2:448010/1‑251 (MQ=255)                                                                                                                                                                                                  attatCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGtt                                   <  1:330083/250‑1 (MQ=255)                                                                                                                                                                                                   ttatCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCCGTAAAATGGTTGc                                 >  1:23460/1‑251 (MQ=255)                                                                                                                                                                                                      gcggaaaaCACCGTAACGGCC‑CAAGCGGGTCGTTCCGCTTGATCTTTAGCGGTCAATGGTCGCTTTTCGACAACTGTGTGGGGCCGCATTACGAGGACGAAAAATAATTATTTGATGAATT‑AGGTAATTAAAGCAAACCTTTATCAAGAAACACAAAAGCAAGGGAGAACGGTTGGGATGAACCGCCGCGCGCTGCGACACAACCTTCACCGTCATCGTGAGCTGGCCCCTGCCAGTAAAATGGTTGCggt                              <  2:418985/246‑1 (MQ=255)                                                                                                                                                                                                       cGGTAAACAACGTAACGGCC‑CAATCGGGACGTTACGCCTGATCTTTAACGGTCAAGGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCggt                              <  1:111241/250‑1 (MQ=255)                                                                                                                                                                                                              caacGTAACGGCC‑CACTCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGaa                        >  1:143890/1‑249 (MQ=255)                                                                                                                                                                                                              aaacgaaccGGCC‑CAATCGGGGCGGTACCCTTGGTCTTTAACGGTCAAGTGTCGCGTTTCGAAAACTAGGCGGGGCCGCAGTAAGAGTACAAATAATAAT‑ATTTTTAGAATTTGGGAAATTAAAGCAAACACTTACAAAGGACGAAAAATGAGAGGGGCAACTGTTGTGATAAACCTCCGCGCTCTGGGTCACAACCGGCAGCGTCTTCGTGAACTGGCCCCTGCCAGTAATATGGGTGCTGTGGTGAAAg                      <  2:244410/243‑1 (MQ=255)                                                                                                                                                                                                               aacGTAACGGCC‑CAATCGGAACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCCGTAAAATGGTTGCGGTGGTGAAAGc                     >  2:370119/1‑251 (MQ=255)                                                                                                                                                                                                                acGTAACGGCC‑CAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGACCTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCg                    >  1:20827/1‑251 (MQ=255)                                                                                                                                                                                                                                 gggACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGACAGCGAACGCTTATGGTCACg    >  1:354417/1‑251 (MQ=255)                                                                                                                                                                                                                                   gACGGTCCGCCTGATCTTAAACGGTCAATGGTCGCTCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTTTCAAGGAACACAAATGCACGCGGCAACTGTTGTGATAAACCGCCGCGCTCTGCGACACACCCTGCAACTTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGt  <  1:577738/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              AAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGT‑ATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCC‑CAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATT‑AGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGT  >  NC_000913/4265430‑4265904

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