Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,905,761	G→A	G25D (GGT→GAT)	mntP →	Mn(2(+)) exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,905,761	0	G	A	97.3%	275.4 / NA	74	G25D (GGT→GAT)	mntP	Mn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/40);  minor base C (2/0);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATG‑AAGATGAAGAGCCGCGCCGTCGACACGGTTTC  >  NC_000913/1905521‑1906002                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   gCAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTc                                                                                                                                                                                                                                          <  1:8304/251‑1 (MQ=255)   aaaaCGTGGCACGTACGGACTGCATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                                                                                                                                                                                                         >  2:139900/1‑250 (MQ=255)     aaCGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                                                                                                                                                                                                      <  2:354901/251‑1 (MQ=255)     aaCGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTGGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                                                                                                                                                                                                      <  2:798851/251‑1 (MQ=255)        gTGGCACGTACGGACTGAATACTTTCAGTCAGGCGCGCCCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCa                                                                                                                                                                                                                                        >  2:126393/1‑246 (MQ=255)        gTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGTGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCa                                                                                                                                                                                                                                        <  1:126393/246‑1 (MQ=255)               gTACGGACTGACTACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCCTGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCCTAAACCGAAAtt                                                                                                                                                                                                                            >  2:814206/1‑251 (MQ=255)                  cGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCGTGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTc                                                                                                                                                                                                                         >  1:789136/1‑251 (MQ=255)                    gACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTg                                                                                                                                                                                                                       >  2:793801/1‑251 (MQ=255)                    gACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTg                                                                                                                                                                                                                       <  2:585514/251‑1 (MQ=255)                     aCTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCCCGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTg                                                                                                                                                                                                                       >  2:936217/1‑250 (MQ=255)                        gAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGaa                                                                                                                                                                                                                     <  1:936217/249‑1 (MQ=255)                        gAATACTTTCAGGCAGGCGAGACCATATGCACCTCAAGCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAg                                                                                                                                                                                                                    <  2:430213/250‑1 (MQ=255)                           tACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                                                                                                                                                                                                                >  1:206365/1‑251 (MQ=255)                                 cagtcagGCGAGACCATATGCCCATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAc                                                                                                                                                                                                          <  2:751736/251‑1 (MQ=255)                                 cagtcagGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAc                                                                                                                                                                                                          >  2:441446/1‑251 (MQ=255)                                  agtcagGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGCTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAcc                                                                                                                                                                                                         >  2:229858/1‑251 (MQ=255)                                  agtcagGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAcc                                                                                                                                                                                                         <  1:940869/251‑1 (MQ=255)                                     cagGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                                                                                                                                                                                                      <  1:139900/251‑1 (MQ=255)                                       ggCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCt                                                                                                                                                                                                    <  1:260403/251‑1 (MQ=255)                                               cATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTg                                                                                                                                                                                            <  2:5388/251‑1 (MQ=255)                                                          cAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGa                                                                                                                                                                                  >  2:18764/1‑250 (MQ=255)                                                                  aTGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAGTGCGAACCGGCCTGATTTTTGGTGCCGTCGAAACCCTGAc                                                                                                                                                                         <  2:449515/251‑1 (MQ=255)                                                                    gCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGAcgc                                                                                                                                                                       <  2:206365/251‑1 (MQ=255)                                                                     ccTGCATGAACGCGAAAAACATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGCCATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTCTGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGAcgcc                                                                                                                                                                      >  2:625153/1‑251 (MQ=255)                                                                         cATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGGTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTGGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTg                                                                                                                                                                  <  2:336882/251‑1 (MQ=255)                                                                            gAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTg                                                                                                                                                                  <  1:675646/248‑1 (MQ=255)                                                                            gAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTg                                                                                                                                                                  >  2:675646/1‑248 (MQ=255)                                                                             aaCGCGAAAACCATCACGGCTTGTTGTTCATGTATGCCTGTTTGCTGGGGGCCATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCg                                                                                                                                                              >  1:50510/1‑251 (MQ=255)                                                                              aCGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGCTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCgg                                                                                                                                                             >  2:161153/1‑251 (MQ=255)                                                                               cgcgAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGc                                                                                                                                                            <  1:728811/251‑1 (MQ=255)                                                                                gcgAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCTCTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGGTGCATCAATCGGTAAAGATGCCACCCTCCATAAGCCGAAAGGTTCTGGAGCATTGCGAACCGGc                                                                                                                                                                                                      <  2:105932/208‑1 (MQ=255)                                                                                gcgAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGCGCGAAGATGTGTTTTATGGATACCCCGATCAGGACATTGTCCTGAATCTCAATGCTACTGATCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                                                                                                                                                                                                      >  1:105932/1‑208 (MQ=255)                                                                                   aaaaCCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTggg                                                                                                                                                        >  1:752090/1‑251 (MQ=255)                                                                                             cGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCAt                                                                                                                                              <  1:793801/251‑1 (MQ=255)                                                                                              ggCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGCTGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATg                                                                                                                                             >  1:812817/1‑251 (MQ=255)                                                                                              ggCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTTCATCAATCGGGAAAGATGCCACCCTcc                                                                                                                                                                                                                                         <  1:326762/159‑1 (MQ=255)                                                                                               gCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGt                                                                                                                                            <  1:441446/251‑1 (MQ=255)                                                                                                   gttgttCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGAGGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagc                                                                                                                                        <  1:50286/251‑1 (MQ=255)                                                                                                     tgttCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagcc                                                                                                                                       <  1:9018/250‑1 (MQ=255)                                                                                                     tgttCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagcc                                                                                                                                       <  1:18764/250‑1 (MQ=255)                                                                                                        tCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGCAATGGGCATGTTagccagcg                                                                                                                                   >  2:677907/1‑250 (MQ=255)                                                                                                            gTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGtt                                                                                                                               <  1:381196/251‑1 (MQ=255)                                                                                                            gTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGtt                                                                                                                               >  1:17149/1‑251 (MQ=255)                                                                                                                gCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGTACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGt                                                                                                                            <  2:812817/250‑1 (MQ=255)                                                                                                                    gTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGTAATGGGCATGTTagccagc                                                                                                                                    >  2:284125/1‑238 (MQ=255)                                                                                                                    gTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagccagc                                                                                                                                    <  1:284125/238‑1 (MQ=255)                                                                                                                          ggtgggCGATGGTGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCTAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGa                                                                                                                 <  2:752090/248‑1 (MQ=255)                                                                                                                             gggCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATc                                                                                                                                                               <  2:242340/202‑1 (MQ=255)                                                                                                                             gggCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATc                                                                                                                                                               >  1:242340/1‑202 (MQ=255)                                                                                                                              ggCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAACCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGCATGGAAcc                                                                                                              >  1:416664/1‑250 (MQ=255)                                                                                                                               gCGATGATGTGTTTTATGGATTCCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCAc                                                                                                            <  1:161153/251‑1 (MQ=255)                                                                                                                                 gatgatGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTg                                                                                                          >  2:341085/1‑251 (MQ=255)                                                                                                                                       gtgtTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCCCTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGc                                                                                                    <  2:36573/251‑1 (MQ=255)                                                                                                                                                ggATACCCCGGTCCGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGGATGTCGAAGGCTGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTCTTTTTGGGGCCGTCGAACCCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTCGCCAGCCGGTTTGTCCTTGAATGGTACCACTGGATTGCGTTtga                                                                                              >  2:651077/1‑247 (MQ=255)                                                                                                                                                ggATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTtgt                                                                                              <  1:651077/248‑1 (MQ=255)                                                                                                                                                  aTACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGtgct                                                                                           >  1:669528/1‑249 (MQ=255)                                                                                                                                                  aTACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGtgct                                                                                           <  2:669528/249‑1 (MQ=255)                                                                                                                                                    aCCCCGGTCAGGACATTGTCATGAATATCACTGCTACGGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGGGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGctgctg                                                                                       <  2:17149/251‑1 (MQ=255)                                                                                                                                                         ggTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATAtt                                                                                  >  1:66504/1‑251 (MQ=255)                                                                                                                                                             aGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTc                                                                              >  1:355385/1‑251 (MQ=255)                                                                                                                                                              ggACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcg                                                                             <  1:188237/251‑1 (MQ=255)                                                                                                                                                               gACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTgggg                                                                                                                                                       <  2:727731/176‑1 (MQ=255)                                                                                                                                                               gACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTgggg                                                                                                                                                       >  1:727731/1‑176 (MQ=255)                                                                                                                                                                       cATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCggcgggcgg                                                                     >  1:894077/1‑249 (MQ=255)                                                                                                                                                                        aggAATATCTCGGCTAATGTTCATCTGGGGTTTGGTATGTCTATGGATGCATTTGCTTCAGCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGTGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGTACCACTGGATTGCGTTTGTGCTGCCGATATTCCTCGGCGGGCGAAt                                                                   <  2:416664/249‑1 (MQ=255)                                                                                                                                                                                 aCTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCCCCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGa                                                          <  1:229858/251‑1 (MQ=255)                                                                                                                                                                                           ttcttcTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGGTTTCGt                                                >  1:541627/1‑251 (MQ=255)                                                                                                                                                                                                     gTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGGTTTCGTGGCGCAgatg                                      >  1:836078/1‑251 (MQ=255)                                                                                                                                                                                                                          aTTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGTTTATTGAGGGTTTTTGTGGCGCCGATGATG‑AAGATGAAGAGCCGCGcc                >  1:639397/1‑251 (MQ=255)                                                                                                                                                                                                                              gctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGGGCATGATTATTTAGGGGTTTCTTGGGGCAGGTGATGGAAG‑TGGAGAGACGCGCCggcg            >  2:523592/1‑248 (MQ=255)                                                                                                                                                                                                                                     caagcGGTAAAGATGCCACCCTCCATAAACCGAAATGTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATGCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATG‑AAGATGAAGAGCCGCGCCGTCGACACGGt     <  1:89342/247‑1 (MQ=255)                                                                                                                                                                                                                                       aTCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATG‑AAGATGAAGAGCCGCGCCGTCGACACGGttt   <  1:341085/251‑1 (MQ=255)                                                                                                                                                                                                                                        gcGGTAAAGTGGCCACCCTACATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATG‑AAGATGAAGAGCCGCGCCGTCGACACGGTTTc  <  2:789136/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   GCAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATG‑AAGATGAAGAGCCGCGCCGTCGACACGGTTTC  >  NC_000913/1905521‑1906002

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