Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	47.2%	46.4 / 12.8	36	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (17/2);  new base C (0/17);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 2.00e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATA  >  NC_000913/2421294‑2421580                                                                                                                                                        |                                                                                                                                          tCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTa                                                                                                                                           <  2:438940/150‑1 (MQ=255) tCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTa                                                                                                                                           <  1:298933/150‑1 (MQ=255)  cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTAc                                                                                                                                          <  2:601031/150‑1 (MQ=255)        aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                     <  2:479752/149‑1 (MQ=255)               acaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgc                                                                                                                            >  2:226577/1‑151 (MQ=255)                 aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                          >  1:244322/1‑151 (MQ=255)                      agaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                     >  2:227848/1‑151 (MQ=255)                               aaggccacaCCCTGCCGTTTTTATTCGTAAAAAGACCCTTCTTTTGCTTCCTTTTATCGAAATA‑AAGGAAAAAAAAATAAACACCCG‑TAAGAATTTTGTTATAATTGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                            <  2:534225/146‑1 (MQ=255)                                 tGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                          >  1:324765/1‑151 (MQ=255)                                       acccTGCCGTTATTATTCATAAAAAGACCCTTCTTTTGCTTACTTTTATCGAAATA‑AAGAAAAAAAAAATAAACACCCG‑TAAGGATATTGTTATCAAGGTGACGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                      <  1:365543/147‑1 (MQ=255)                                             ccGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                              >  1:241376/1‑151 (MQ=255)                                              cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                             >  2:611270/1‑151 (MQ=255)                                                    tattCGTAAAAAGACCCTTCGTTTGCGTACTTTTATCGAAATA‑AAGGAAAAAAAAATAAACACCCG‑TCAGGATATTGTTATAAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                        <  1:453722/150‑1 (MQ=255)                                                       tCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                    >  1:638546/1‑151 (MQ=255)                                                        caaaaaaaGACCCTTCTTTGGCGTACTTTTATCGAAATA‑AGGAAAAAAAAAATAAACACCCG‑AAAGAATTTTTTTACAAATGTTACGCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                    <  1:412572/145‑1 (MQ=255)                                                           ataAAGACCCTTCTTTGGCTTATTTTTATGGAAATA‑AAGGAAGAAAAAATAAACACCCG‑TCAGAATATTGTAATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                   <  1:504926/148‑1 (MQ=255)                                                           aaaaaaacccTTTGTTTGCGTATTTTTATCGAAATA‑AAGAAAAAAAAAATAACCACCCG‑TCAGGATTTTTTTATAAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                    <  1:397592/139‑1 (MQ=255)                                                            taAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                               >  1:404601/1‑151 (MQ=255)                                                              aaaacccTGCGTTTGCTTACTTTTATGGAAATA‑AAGAAAAAAAAAATAAACACCCG‑ACAGAATTTTGTTAAAAATGTGACTACCCATCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                             <  1:288118/146‑1 (MQ=255)                                                              aaGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                              >  2:90325/1‑150 (MQ=255)                                                                aaccctttttttGCGTCCTTTTATCGAAATT‑AAGAAAGAAAAATTAAACACCCG‑TAAGTATATTTTAAAAAATGCGACGTTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                           <  1:275326/142‑1 (MQ=255)                                                                     ttCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                       >  1:151372/1‑150 (MQ=255)                                                                      ttttttGCGTCCGTTTAT‑GGAATTAAAGAAAAAAAAAATAAACACCCG‑TAAGAATTTTTTTATAAATGTCCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaaa                                                                      <  1:234520/147‑1 (MQ=255)                                                                                  ttttATC‑GAAAT‑AAAGAAAAAAAAAATAACCACCGATCAGAATTTTTTTATCAATGTGCCAATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa                                                          <  1:588167/149‑1 (MQ=255)                                                                                   tttactgaaaTT‑AAGAAATAAGAAATCAACAACCG‑GAAGGATTTGTTTATCAGTGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                            <  1:279068/139‑1 (MQ=255)                                                                                      aTCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                       >  1:493720/1‑149 (MQ=255)                                                                                           aaTA‑AAGAAAAAAAAAATAAACACCCG‑TAAGGATATTTTTATAAATGTGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                 <  1:35347/150‑1 (MQ=255)                                                                                                    aagaaaaaATAAACACCCG‑GCAGAATGTTGTTATAAAGGGGCCCCCCCATTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                            <  2:205624/147‑1 (MQ=255)                                                                                                          aaaTTAACAACCG‑TAAGGATATTGTTATCAATGTGACGACCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                  <  1:292988/151‑1 (MQ=255)                                                                                                            aTTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGc                                >  1:618736/1‑151 (MQ=255)                                                                                                              tAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                              >  2:656479/1‑151 (MQ=255)                                                                                                               aacaccCT‑TAAGGTTATTTTTAAAAATGTGACGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                               <  1:389305/149‑1 (MQ=255)                                                                                                               aacaacCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                             >  1:2178/1‑151 (MQ=255)                                                                                                                acaacCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                             >  2:111093/1‑150 (MQ=255)                                                                                                                  acccG‑TAAGGATATTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                           <  2:224201/148‑1 (MQ=255)                                                                                                                                      tCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt       >  2:630336/1‑151 (MQ=255)                                                                                                                                          tgtgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGat   >  1:444969/1‑151 (MQ=255)                                                                                                                                           gggacgacCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGata  <  1:550970/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                          TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCG‑TCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATA  >  NC_000913/2421294‑2421580

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