Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,511,772	0	A	C	42.2%	82.9 / 10.5	45	V312G (GTG→GGG)	mntH	manganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base A (26/0);  new base C (0/19);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 4.10e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGT  >  NC_000913/2511624‑2511922                                                                                                                                                       |                                                                                                                                                       agCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCaccacc                                                                                                                                                      >  1:217604/1‑151 (MQ=255) agCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCaccacc                                                                                                                                                      >  1:347060/1‑151 (MQ=255)        ccAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCg                                                                                                                                               >  1:195052/1‑151 (MQ=255)         cAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGc                                                                                                                                              >  2:857911/1‑151 (MQ=255)          aaaaCTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGcc                                                                                                                                             >  1:276162/1‑151 (MQ=255)                            aCTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGtt                                                                                                                           >  1:372133/1‑151 (MQ=255)                                aTAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGa                                                                                                                         >  1:333709/1‑149 (MQ=255)                                aTAACAAAAACCCGTGTCGAATCTAACCCCACCAGAATAACAAAAAAGGACGCCAACAGGGTGCTTGTCCGCCCCCCCCCCGGGGGGTTATAAAAGC‑G‑AAGAAACCCCTCCATCCCCCCCGGCCCCGCCAGTGTCCCCACCACCGtt                                                                                                                           <  2:814755/147‑1 (MQ=255)                                 tAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGa                                                                                                                         >  2:581895/1‑148 (MQ=255)                                  aaCCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGAc                                                                                                                     >  1:354820/1‑151 (MQ=255)                                   aCCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGAc                                                                                                                     >  2:835701/1‑150 (MQ=255)                                    ccaaactCCGTGTCGGATCAACTCCCAACAAAAAAAAAAAAAATGGCGACAACAGTGGCACTGTACGCCCCCCCAACAGCGGTATATAAAAGC‑G‑AATGAACCCCGGTATCCCCCCCCCCCCCGCCAGTGTCCCCACCACCGTTGAAg                                                                                                                       <  2:486299/141‑1 (MQ=255)                                      aaaaTCCGTTGCGCATCAAACCCCAGCAAAATGAAAAAAAATGACGAAAACATGGTGACTGTACGACGCACCCCCAGCGGGATA‑GAAAGCGG‑AATGAGGCCGGGCACCACCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAg                                                                                                                   <  2:441520/147‑1 (MQ=255)                                      aGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGcc                                                                                                                 >  2:588822/1‑151 (MQ=255)                                        aacccGGGTCGGATCAAACCCCATCAGAATAAAAAAAAATGCGGGCAACATGGGGCCTGTCCGACGCCCCCCCGGGGGGATATG‑AAGC‑GAAATGAGGCCCGGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGcc                                                                                                                 <  2:885715/146‑1 (MQ=255)                                                        aaTCCCACAAGAATGACAAAAAATGACGCCAACTTGGTGACTGTACGACGCACCCACAGCGGGATAAGAAAGC‑G‑AATAAAGCCCTGCTCCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCtt                                                                                                <  1:495707/150‑1 (MQ=255)                                                         aTCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTa                                                                                               >  1:252380/1‑150 (MQ=255)                                                          tCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAAt                                                                                             >  2:384061/1‑151 (MQ=255)                                                             cATCAGAATAAAAAAAAATGACGGCAACATGGTGACTGTACGCCGCCCCCCCGGGGGGATAGAAAAGC‑G‑AATAAACCCCTCCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCa                                                                                          <  1:583548/151‑1 (MQ=255)                                                              aTCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCa                                                                                          >  1:782506/1‑150 (MQ=255)                                                                cAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCaa                                                                                         >  1:485331/1‑149 (MQ=255)                                                                 aGAATGACAATAAATGACGGCAACAGGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCaaa                                                                                        <  1:41859/149‑1 (MQ=255)                                                                   aaggACAATAAATGACGGCAACATGGTGACTGTACGACGCCCCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGAc                                                                                     <  1:617104/147‑1 (MQ=255)                                                                   aaTGACAAAAAATGCCGGCAACAGGGTGACTGTACGCCGCCCCCACAGCGGGATATGAAAGC‑G‑AATGAACCCCTGCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGAc                                                                                     <  1:248172/150‑1 (MQ=255)                                                                             aaTGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGcc                                                                              >  2:614941/1‑147 (MQ=255)                                                                               tGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCg                                                                        >  2:694389/1‑151 (MQ=255)                                                                                   ggCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGc                                                                    >  2:580639/1‑151 (MQ=255)                                                                                      aaCCTGGTGCCTGTCCGCCGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCt                                                                   <  1:700989/149‑1 (MQ=255)                                                                                           ggTGACTGTACGACCCCCCCACGGCGGGATATGAAAGC‑A‑AATGAACCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACa                                                              <  2:544039/149‑1 (MQ=255)                                                                                           ggTCACTGTCCAACGCCCCAACAGCGGAATATGAAAGC‑A‑AATGAAGCCCTGCAACCCCCCCAACCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAg                                                             <  1:450376/150‑1 (MQ=255)                                                                                              gccTGTACGCCCCCCCCACGGCGGGATAGAAAAGC‑G‑AATGACCCCCTGCACCCCCCCCGCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCgg                                                          <  2:364737/148‑1 (MQ=255)                                                                                                 tGTACGCCGCCCCCACAGCGGAATATAAAAGC‑G‑AATGAACCCCTCCATCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTaa                                                                <  1:559263/141‑1 (MQ=255)                                                                                                      gacgCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCg                                                  >  2:449494/1‑150 (MQ=255)                                                                                                         gCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTc                                                >  1:845436/1‑149 (MQ=255)                                                                                                          ccccccCAGGGGGATATAAAAGC‑G‑AATGACCCCCGCCCCCCCCCCCGCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAg                                                    <  2:743285/142‑1 (MQ=255)                                                                                                          cccccACAGCGGGATATGAAAGC‑A‑AATGAACCCCTGCCCCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGa                                             <  2:68573/149‑1 (MQ=255)                                                                                                             ccACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATa                                           >  2:283952/1‑150 (MQ=255)                                                                                                                            aaaGC‑G‑AATGAAGCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGAt                             <  1:123621/149‑1 (MQ=255)                                                                                                                            aaaGC‑A‑AATGAACCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATc                            <  1:30928/150‑1 (MQ=255)                                                                                                                              aGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGc                         >  2:254306/1‑151 (MQ=255)                                                                                                                               gC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGc                         >  2:53281/1‑150 (MQ=255)                                                                                                                                  g‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAac                      >  1:56597/1‑151 (MQ=255)                                                                                                                                           cccTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCa                  <  2:373354/147‑1 (MQ=255)                                                                                                                                            ccTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATg             >  2:835978/1‑151 (MQ=255)                                                                                                                                             cTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATg             >  2:217883/1‑150 (MQ=255)                                                                                                                                               gCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGAcc          >  1:238575/1‑151 (MQ=255)                                                                                                                                                  tCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCgg        >  2:385454/1‑150 (MQ=255)                                                                                                                                                    accaccTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGaa      >  1:368441/1‑150 (MQ=255)                                                                                                                                                       accTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGt  >  1:239866/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AGCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGC‑G‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGT  >  NC_000913/2511624‑2511922

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