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.4%	237.2 / NA	80	T452I (ACC→ATC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (40/36);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAG  >  NC_000913/4265425‑4265913                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      tGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCCGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAAc                                                                                                                                                                                                                                                >  1:155425/1‑251 (MQ=255)   gCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATGATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACgg                                                                                                                                                                                                                                              <  1:517477/251‑1 (MQ=255)     gAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGCGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTc                                                                                                                                                                                                                                            >  1:361664/1‑251 (MQ=255)       aaGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATATGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCaa                                                                                                                                                                                                                                          <  1:468222/251‑1 (MQ=255)         gAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATg                                                                                                                                                                                                                                        >  2:57714/1‑251 (MQ=255)                      cGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAACTTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGAc                                                                                                                                                                                                                           >  1:317896/1‑251 (MQ=255)                        gtggtgGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                                                                                                                                                                          >  1:385467/1‑250 (MQ=255)                          ggtggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACt                                                                                                                                                                                                                       >  1:86465/1‑251 (MQ=255)                          ggaggCGATGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGCGGTGTATCACGAAAACAGTGATTTAAAAGGCATAGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCAAATCGGTCAGGAACGCCTGATCTTTATCGGGAAATGGTCGCGATTAGACAAc                                                                                                                                                                                                                        >  2:303888/4‑250 (MQ=255)                                gcTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgg                                                                                                                                                                                                                 >  1:535570/1‑251 (MQ=255)                                   ggcccAGTTGAACCGTTCTCTGGAACGACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTAACGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgta                                                                                                                                                                                                                <  2:50715/247‑3 (MQ=255)                                   gTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGCAATTATTATCGGTAAACAACGTAACGGACCAATCGGGCCGGTACGCCTGTTCTTTAACGGTCAATGGACGCGCTTCGACAACTATACggggg                                                                                                                                                                                                              >  2:312507/1‑250 (MQ=255)                                   gTCCCAGTTGAACCGTTCTCTGGAAATACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCATGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgggg                                                                                                                                                                                                               <  2:389558/250‑1 (MQ=255)                                      ccAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGAAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCACTGATCGCGCTTCGACAACTATGCGGGGCCGc                                                                                                                                                                                                           >  2:157202/1‑251 (MQ=255)                                       cAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCa                                                                                                                                                                                                          <  1:527301/251‑1 (MQ=255)                                            gAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTac                                                                                                                                                                                                      <  1:191635/250‑1 (MQ=255)                                              aCCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgac                                                                                                                                                                                                   <  1:14773/251‑1 (MQ=255)                                                   tcgctGGAACAAAGTGCCGACAAACGCCCGGTCAACTCCGACCTGAGTGAATCTGGTGCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTAGTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacgacg                                                                                                                                                                                               <  2:421299/247‑1 (MQ=255)                                                        ggAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGGCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataa                                                                                                                                                                                          <  2:291137/250‑1 (MQ=255)                                                         gAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                                                                                                                                                                        <  1:560193/251‑1 (MQ=255)                                                                gTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTAttt                                                                                                                                                                                 >  1:120708/1‑251 (MQ=255)                                                                gTGCCGACAAACGCCCGGTCAACTACGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTAttt                                                                                                                                                                                 >  1:490867/1‑251 (MQ=255)                                                                 tGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTAtttt                                                                                                                                                                                >  1:171760/1‑251 (MQ=255)                                                                 tGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTAtttt                                                                                                                                                                                >  1:194079/1‑251 (MQ=255)                                                                 tGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTAtttt                                                                                                                                                                                >  1:480635/1‑251 (MQ=255)                                                                     gACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGa                                                                                                                                                                            <  1:88488/251‑1 (MQ=255)                                                                          aCGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAtt                                                                                                                                                                         >  1:425623/1‑249 (MQ=255)                                                                             cccgggaaatgccGACCTGAGTGAACCAGGGTCTAACGAGCAGGATGCGGACGTGAACATTATTAACTATAGTGATGAGGTGTATCACGAAAACAGTAATTTAAAAGGCATCGC‑GGAAATTATTAACGTTAAACAACGAAACGGCCCAATCGGGACGTTACGCCTGATCTTTAACGGTCAATAGTAGCAATTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAttcga                                                                                                                                                                      <  2:77049/238‑4 (MQ=255)                                                                                ggTCAACTCCGACCTGCGTGAATCTGGCTCTATCTAGCAGGCTGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGAAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGAAAATGGTCGCGCTTCGCCAACTATGCGGGCCCGCAGTACGAAGACGAATAATAATTATTTTATAATTTGGGTAAtt                                                                                                                                                                 >  2:386216/1‑251 (MQ=255)                                                                                    aaCTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAg                                                                                                                                                             <  2:535570/251‑1 (MQ=255)                                                                                        ccGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCa                                                                                                                                                           <  1:430353/249‑1 (MQ=255)                                                                                                 gTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCa                                                                                                                                                <  1:171815/251‑1 (MQ=255)                                                                                                    aaTCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAgg                                                                                                                                             <  2:120708/251‑1 (MQ=255)                                                                                                         ggCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATGTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGa                                                                                                                                            <  2:254694/247‑1 (MQ=255)                                                                                                          gCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAAcaca                                                                                                                                       <  1:181195/251‑1 (MQ=255)                                                                                                          gCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAAcaca                                                                                                                                       <  2:432736/251‑1 (MQ=255)                                                                                                              tATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTTATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAAt                                                                                                                                    <  2:134543/250‑1 (MQ=255)                                                                                                                   aGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGc                                                                                                                              >  1:188017/1‑251 (MQ=255)                                                                                                                    gtctgATGCGGACATGATCATGTTTATCATTTGTTAGGTGTTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATGATTATCGTTAAACAACGTACCGGCCCAATCGGGACGGTACGCCTTATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGAAGACGAATAATAATTATTGTATGAATTAGGTAATTAAAGCAAACACTTAGCAAGGAACCCAAATGCAAGCg                                                                                                                             <  2:335561/247‑1 (MQ=255)                                                                                                                    gCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCg                                                                                                                             >  1:80168/1‑251 (MQ=255)                                                                                                                      aGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGc                                                                                                                           <  1:350324/251‑1 (MQ=255)                                                                                                                             ggACTTGATCATGTTTAGCTATCGTGATGAGGTTTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtgtt                                                                                                                    <  2:91322/251‑1 (MQ=255)                                                                                                                              gACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtgttg                                                                                                                   <  1:549638/251‑1 (MQ=255)                                                                                                                                cTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCAGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATAAAGGAACACAAATGCAAGCGGCAACTGTtgtg                                                                                                                 >  2:292753/1‑251 (MQ=255)                                                                                                                                 ttGATCATGTTTATCTATCGTGATGAGTTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTAACGGTAAACAATGAAACGGCCCAATCGGGCCGGGACGCCTGTTATTTAACGGTCAATGGTCGCGCTTCGCCAACTATGCGGGGCAGAAGAACGACGAGGAATCAAAATTATTTTATTAATTAGGTAATTAAAGCAAACACTTATAAAGGAACAAAAATgaaatgggcaaatttagct                                                                                                                 >  2:33976/1‑232 (MQ=255)                                                                                                                                    aTCATGTTTATATATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCTTCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTAAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTCAGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAACCACTTATCAAGGAACACAAATGCAAGCGGCAACTGTtgtgctt                                                                                                              >  2:79738/1‑247 (MQ=255)                                                                                                                                           ttATCTATCGTTATGAGGTGTCTCACGAAAACATTGATTTAAAAGGCCTCGC‑GGAAATTATTATCGGTAAACACCGTAACGGCCCAATCGGGCCGGTACTCCTGTTCTTTAACGGTCAATGGTCGCGCTTCGACACCTATACGGCGCCCCATTACGTCGACGACTCATAATTATTTTATGAATTAGGTAATTAAAGCAACCACTTATCAAGTAACACAAATGCACGCGTCAACTGTtgtgtttcaccaccc                                                                                                      >  2:295816/1‑240 (MQ=255)                                                                                                                                                       atgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGac                                                                                           >  1:468098/1‑250 (MQ=255)                                                                                                                                                          gggTGTAGCACGAAAACAGTGATTTAAAAGGAAGCGA‑GGAAATTATTATCGGAAAACAACGTAACGGCCCAATGGGGACGGTGCGCCAGATCTTTAACGGTAAATGGGCGCGCTTCGACAACTATGCGGGGCCGAAGTACGACAACGAATAATAATAATTTGAAGAATTAGGTAATTAAAGGAAACACTAATAAAGGAACCCAAAGGCAAGCGGCAACTCTTGTGATTAACCGCCGCGCTCTGCGACACaa                                                                                       <  1:295816/250‑1 (MQ=255)                                                                                                                                                          aGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTCCGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGAcaca                                                                                        >  1:345411/1‑250 (MQ=255)                                                                                                                                                            gtgtATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAAcc                                                                                     >  1:16501/1‑251 (MQ=255)                                                                                                                                                                  cACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCaa                                                                                <  1:534732/250‑1 (MQ=255)                                                                                                                                                                    cGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAAt                                                                                                                                    >  1:210738/1‑196 (MQ=255)                                                                                                                                                                     aaaaaCAGGGATTTAAAAGGCAT‑GCATGAATTTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCGGCGACACAACCTGCAACGTc                                                                            <  2:171760/250‑1 (MQ=255)                                                                                                                                                                      aaaaCAGTGATTTACAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCACTCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCt                                                                           >  2:46216/1‑251 (MQ=255)                                                                                                                                                                       aaaCAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTATGATTAACCGCCGCGCTCTGCGACCCAACCGGCAACGTCtt                                                                          >  2:371326/1‑251 (MQ=255)                                                                                                                                                                       aaaCAGTGATTTAAAAGGCAGCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTAAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATAAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCtt                                                                          <  1:232147/251‑1 (MQ=255)                                                                                                                                                                            gTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGa                                                                     <  1:388006/251‑1 (MQ=255)                                                                                                                                                                                 taaaacgtcATCGC‑GGAAATTATTGTCGGTGACAAACGATACGTCCCAATCGGGCTGGTCCCCCTGATCTTTACCGTACAATGGTTGCGCTCTGACAACTATGCTGGGCCGCATTACGACGACGAATAATAATTATTTTATGAATAAGGTAATTAAAGCAAACACTTATCAAGGAACACAATTGCAAGCGGCAACTGTTGTGTTAAACCGCCGCGCTCTTCGATCCAACCTGCAACTTCTTCGTGAACTg                                                                 <  2:392456/242‑1 (MQ=255)                                                                                                                                                                                   aaaaGGCATCGC‑GGAAATTATTAGCGGTAAACAACGTAACGGCCCAATCGGGACGGTACCCCTGACCTTAAACGGTCAATGGTCGCGCTGCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGcc                                                              <  2:80168/251‑1 (MQ=255)                                                                                                                                                                                    aaaGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGccc                                                             <  1:312507/251‑1 (MQ=255)                                                                                                                                                                                      aggtaTCGC‑GGAAATTATTATCGGTAATCAACGTAACGGCCCAACCGGGACGGTACTCCTTATCTTTAACGGTCAATGGTCGCGCTTCGACACCTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTGGGTAATTAAAGCAATCCCGTATCAAGGAACAAAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGccccg                                                           <  2:385467/247‑2 (MQ=255)                                                                                                                                                                                      aGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCt                                                           >  1:539563/1‑251 (MQ=255)                                                                                                                                                                                       ggCATCGC‑GGAAATTATTATCTGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTg                                                          >  1:150982/1‑251 (MQ=255)                                                                                                                                                                                        gCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACCCAACCTGCAACGTCTTCGTGAACTGGCCCCTGc                                                         >  1:234864/1‑251 (MQ=255)                                                                                                                                                                                              c‑gGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGAACTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAGTAGGTAATTAAAGCAAACACTTATCACGGAACACAAATGAAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTAGCTCCGGCCAGTa                                                    >  2:108058/1‑250 (MQ=255)                                                                                                                                                                                                   aaTTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATg                                               >  1:483359/1‑251 (MQ=255)                                                                                                                                                                                                        tttgcGGTAAACAACTTACCGGCCCAACGGGGACGGTGCGCTTGACTGTAAACGGTCAATGTCCGCGTTTCGAAAACTATGGGGGGCCGAAGTACGACGACAAATAATAATTATTTTATGAATTAGGAAATAAAAGCAAACCCTTATCAAGTAACACAAATGCAAGGGACACGTGTTGTGATTAACCGCCGCGCTCGGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTg                                           <  2:194079/246‑1 (MQ=255)                                                                                                                                                                                                         tatCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACTTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCg                                         >  2:137090/1‑251 (MQ=255)                                                                                                                                                                                                            gggTAAAAAACGTAACGGCCAAAGGGGGAGGGTACGCGTGTACTTTAACGGTCAATGGTCGCGATTCGAAAACAATGCGGGGCCGAAGTACGACGAAGAATAAAAATTATGTTATAAATTAGTAAATAAAAGAAAACACTAAAAAAGAAACACAAAGGGAAGCGACAACTGGTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCggtg                                      <  1:273320/250‑1 (MQ=255)                                                                                                                                                                                                            cGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGTTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTACCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCggt                                       <  1:8178/250‑1 (MQ=255)                                                                                                                                                                                                             ggTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCggtgg                                     >  1:123642/1‑251 (MQ=255)                                                                                                                                                                                                              gTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCggtggg                                    >  1:493404/1‑250 (MQ=255)                                                                                                                                                                                                                aaaCAACGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGgtggtg                                   >  1:451395/1‑250 (MQ=255)                                                                                                                                                                                                                    aacGTCACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCCGTAAAATGGTTGCGGTGGTGaacgc                              >  2:362506/1‑248 (MQ=255)                                                                                                                                                                                                                    aacGTAACGGCCCAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGc                              >  1:120117/1‑251 (MQ=255)                                                                                                                                                                                                                               ccAATCGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATg                   >  1:252580/1‑251 (MQ=255)                                                                                                                                                                                                                                    cGGGACGGTACGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACCCTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCa               <  1:195583/250‑1 (MQ=255)                                                                                                                                                                                                                                              cGCCTGATCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAGGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTg    <  1:46216/251‑1 (MQ=255)                                                                                                                                                                                                                                                ccTGATCTTTAACGGCAAATGGTCGCGCTTCGACAACTAGGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTAGGAATTAGGTAATTAAAGCAAACTCTTATCAAGGAACACAAAGGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTgag  <  1:386216/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TGGCGAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGC‑GGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGACACAACCTGCAACGTCTTCGTGAACTGGCCCCTGCCAGTAAAATGGTTGCGGTGGTGAAAGCGAACGCTTATGGTCACGGTCTTCTTGAG  >  NC_000913/4265425‑4265913

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