Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,841,822	G→A	5.9%	Q75* (CAG→TAG)	ldhA ←	D‑lactate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,841,822	0	G	A	5.9%	249.4 / 4.1	85	Q75* (CAG→TAG)	ldhA	D‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (42/38);  new base A (3/2);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GCCACGGCATGCTCGGCAACGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGAAGCTTTCCTGA  >  NC_002947/1841674‑1841970                                                                                                                                                     |                                                                                                                                                      gcCACGGCATGCTCGGCAACGGCGTGCGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTgc                                                                                                                                                      <  1:35386/150‑1 (MQ=255)  ccACGGCATGCTCGGCAACGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACAGCGTGCCAGCGGCCGCCAGGCGCTgca                                                                                                                                                     >  2:1183518/1‑150 (MQ=255)      ggCATGCTGGGCAACCGGGTGTGGCGAGTAGGCCGGCAGCCGCACCACGGCCAGCCGCAAGCGCTGAGCGGCGGCCAGGTCGACGTGGTTTTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCAc                                                                                                                                                 >  2:832473/1‑150 (MQ=255)       gCATGCTCGGCAACGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGTGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCAcc                                                                                                                                                <  2:1969196/150‑1 (MQ=255)          tGCTCGGCAACGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCC‑GGTGCTgcagcgcctggc                                                                                                                                            <  2:350058/150‑8 (MQ=255)              cGGCAACGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGt                                                                                                                                         <  2:1080029/150‑1 (MQ=255)               ggCAACGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTc                                                                                                                                        <  2:26246/150‑1 (MQ=255)                gCAACGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTAGACGTGGTTGTAGCCGGCCGAAACCGCGGCGAGAAACCGCGTGCAAGCGGCGGCAAGTCTCTGCAGCACCGGCGCGttg                                                                                                                                       >  2:733267/1‑148 (MQ=255)                  aaCGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAg                                                                                                                                     <  1:484859/150‑1 (MQ=255)                      gCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCg                                                                                                                                 >  2:1200611/1‑150 (MQ=255)                        gtgtGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTc                                                                                                                               >  1:1479260/1‑150 (MQ=255)                            ggCGAGTAGGCCGGCACCCGCACCACTGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATtg                                                                                                                           <  1:2124030/150‑1 (MQ=255)                            ggCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATtg                                                                                                                           <  1:1420885/150‑1 (MQ=255)                            ggCGAGTAGGCCGGCACCCGCACCACGGCCAACCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCCAGCGCCGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATtg                                                                                                                           >  2:309330/1‑150 (MQ=255)                             gCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATtga                                                                                                                          >  2:563520/1‑150 (MQ=255)                                 gTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGaa                                                                                                                      <  2:710910/150‑1 (MQ=255)                                 gTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTACAGCACCGGCGCGTCGAGTTCGTCATTGATGaa                                                                                                                      >  2:590508/1‑150 (MQ=255)                                   aGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAgg                                                                                                                    >  1:1688269/1‑150 (MQ=255)                                    ggCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGgc                                                                                                                   >  2:987315/1‑150 (MQ=255)                                       cGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTACAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCa                                                                                                                >  2:711678/1‑150 (MQ=255)                                           aCCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGacc                                                                                                            >  2:535035/1‑150 (MQ=255)                                           aCCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGacc                                                                                                            >  1:863521/1‑150 (MQ=255)                                             ccGCACCACGGCCAGCCCCAAGCGCTTGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGaccac                                                                                                          >  2:1189709/1‑150 (MQ=255)                                             ccGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGaccac                                                                                                          >  2:210015/1‑150 (MQ=255)                                             ccGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGaccac                                                                                                          >  1:1952631/1‑150 (MQ=255)                                             ccGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGACGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGaccac                                                                                                          >  1:2234576/1‑150 (MQ=255)                                              cGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGaccacc                                                                                                         >  1:2053155/1‑150 (MQ=255)                                                caccacGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTc                                                                                                       <  2:1468336/150‑1 (MQ=255)                                                   cacGGCCAGTCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCaaa                                                                                                    <  1:1049743/150‑1 (MQ=255)                                                   cacGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCaaa                                                                                                    <  1:1538518/150‑1 (MQ=255)                                                   cacGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCaaa                                                                                                    <  2:1881596/150‑1 (MQ=255)                                                      ggCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑cc                                                                                                >  2:747510/1‑150 (MQ=255)                                                      ggCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑cc                                                                                                <  2:516886/150‑1 (MQ=255)                                                      ggCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCGCAAAG‑cc                                                                                                <  2:516491/150‑1 (MQ=255)                                                        ccAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCAt                                                                                              <  1:628551/150‑1 (MQ=255)                                                        ccAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCAt                                                                                              <  1:1200611/150‑1 (MQ=255)                                                         cAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATc                                                                                             >  2:213348/1‑150 (MQ=255)                                                         cAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATc                                                                                             >  2:691317/1‑150 (MQ=255)                                                            ccccAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGc                                                                                          >  1:682272/1‑150 (MQ=255)                                                             cccAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGcc                                                                                         >  2:1188704/1‑150 (MQ=255)                                                             cccAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCATCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGcc                                                                                         >  2:428449/1‑150 (MQ=255)                                                              ccAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCa                                                                                        >  2:474323/1‑150 (MQ=255)                                                                aaGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGaccac                                                                                                          <  1:1633016/131‑1 (MQ=255)                                                                aaGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGaccac                                                                                                          >  2:1633016/1‑131 (MQ=255)                                                                aaGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAgg                                                                                      >  2:1049662/1‑150 (MQ=255)                                                                     cttggcggcggcCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGggcggc                                                                                 <  1:1189709/147‑1 (MQ=255)                                                                     cTGGGCAGCGACCAGGTCCACGGGGTTAGAGCGGGCCGAGCGCAGGGCGATCAACCGCGTGCCACAGGGCTCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGggcggc                                                                                 <  1:514224/150‑1 (MQ=255)                                                                        ggcggcggcCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGgt                                                                              <  2:1479260/150‑1 (MQ=255)                                                                          cggcggcCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGgtgt                                                                            <  1:1146728/150‑1 (MQ=255)                                                                            gcggcCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCg                                                                          >  1:1720822/1‑150 (MQ=255)                                                                              ggcCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAg                                                                        <  2:7527/150‑1 (MQ=255)                                                                                ccAGGTCGCGGTGGTTGTAGCCGGCCGAGCGCAGGGCGAACAACCGCGTGCCTGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGCTCGCCACTGATGAAGGCGCTGACCACCTCAACG‑CCATCTGCCAGGGCGGCGGTGTCGgggg                                                                      >  2:595621/1‑150 (MQ=255)                                                                                    gTCAACGTTGTTGTATCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCCGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGCCGCAGACCACCTCAAAG‑CCATCTGCCAGGTCGGCGGTGTCGCGGGTCAg                                                                  >  2:1043815/1‑150 (MQ=255)                                                                                          gtcgTTGTAGCTGGCCGAGCGCAGGGCGATCAACCGCGTGCGAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAggcgggc                                                            >  2:686216/4‑150 (MQ=255)                                                                                          gTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAggcgggc                                                            >  1:585984/1‑150 (MQ=255)                                                                                             gTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAGTCGATGT‑CCAGGGCGGCGGTGTCGAGGGTCAGGcgggcggg                                                         >  1:1288198/1‑150 (MQ=255)                                                                                                   gccggccgAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTGATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGAGCTGGaa                                                   >  1:490012/1‑150 (MQ=255)                                                                                                         cgAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAg                                             <  2:41194/150‑1 (MQ=255)                                                                                                            gcgcAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTc                                          <  1:2278766/150‑1 (MQ=255)                                                                                                                  ggCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGc                                    >  1:1775849/1‑150 (MQ=255)                                                                                                                    cGATCAACCGCGTGCCTGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCgg                                  <  1:1981614/150‑1 (MQ=255)                                                                                                                        cAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTgcc                              <  2:585984/150‑1 (MQ=255)                                                                                                                         aaCCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTgccg                             >  1:1742634/1‑150 (MQ=255)                                                                                                                           ccGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTgccgcc                           >  2:1825249/1‑150 (MQ=255)                                                                                                                            cgcgTGCCAGCGGCGGCCAGGCGCTACAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAgg                                      >  1:1067597/1‑138 (MQ=255)                                                                                                                            cgcgTGCCAGCGGCGGCCAGGCGCTACAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAgg                                      <  2:1067597/138‑1 (MQ=255)                                                                                                                               gTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCa                       >  2:235949/1‑150 (MQ=255)                                                                                                                               gTGCCAGCGGCGGCCAGGCGCTGCAGCACCCGCGCGTCGAGTTCGTCATTAATGAAGGCGCAGACCCCCTCAAAG‑CCAGCTTCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCa                       <  1:832473/150‑1 (MQ=255)                                                                                                                                 gCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGc                     >  1:2268169/1‑150 (MQ=255)                                                                                                                                  ccAGCGGCGGCCAGGCGCTACAGATTCTGCTCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGcc                    <  1:590508/150‑1 (MQ=255)                                                                                                                                   cAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGt                                 >  2:1727579/1‑136 (MQ=255)                                                                                                                                   cAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGt                                 <  1:1727579/136‑1 (MQ=255)                                                                                                                                     gcggcgGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTg                 <  2:1688269/150‑1 (MQ=255)                                                                                                                                      cggcggCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTgg                <  1:213348/150‑1 (MQ=255)                                                                                                                                       ggcggcCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGt               <  1:747510/150‑1 (MQ=255)                                                                                                                                       ggcggcCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGt               <  1:1049662/150‑1 (MQ=255)                                                                                                                                          ggcCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGTTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGaa            <  2:2234576/150‑1 (MQ=255)                                                                                                                                          ggcCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGaa            <  2:863521/150‑1 (MQ=255)                                                                                                                                             cAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCa                                              <  1:2052483/113‑1 (MQ=255)                                                                                                                                             cAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCa                                              >  2:2052483/1‑113 (MQ=255)                                                                                                                                              aGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGAAGCtt        >  1:2041168/1‑150 (MQ=255)                                                                                                                                              aGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCCGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGAAGCtt        >  1:238218/1‑150 (MQ=255)                                                                                                                                                gcgcTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCTGCAGCCTTGGTGAAGCTTTc      <  2:968263/150‑1 (MQ=255)                                                                                                                                                gcgcTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGAAGCTTTc      <  2:2053155/150‑1 (MQ=255)                                                                                                                                                  gcTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGAAGCTTTCCt    <  1:309330/150‑1 (MQ=255)                                                                                                                                                    tGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGAGCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGAAGCTTTCCTGa  >  2:1690857/1‑150 (MQ=255)                                                                                                                                                    tGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCa                                        <  1:206810/112‑1 (MQ=255)                                                                                                                                                    tGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑ACATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCa                                        >  2:206810/1‑112 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCCACGGCATGCTCGGCAACGGCGTGTGGCGAGTAGGCCGGCACCCGCACCACGGCCAGCCCCAAGCGCTGGGCGGCGGCCAGGTCGACGTGGTTGTAGCCGGCCGAGCGCAGGGCGATCAACCGCGTGCCAGCGGCGGCCAGGCGCTGCAGCACCGGCGCGTCGAGTTCGTCATTGATGAAGGCGCAGACCACCTCAAAG‑CCATCTGCCAGGGCGGCGGTGTCGAGGGTCAGGCGGGCGGGCTGGAAATGCAGCTCCAGGGCGGTGCCGCCGGCAGCCTTGGTGAAGCTTTCCTGA  >  NC_002947/1841674‑1841970

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