Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,394,579	0	G	C	40.0%	91.9 / 21.8	50	G172R (GGC→CGC)	nnr	bifunctional NAD(P)H‑hydrate repair enzyme; C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase
Reads supporting (aligned to +/- strand):  ref base G (8/22);  new base C (20/0);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 1.26e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGCATTTTGACTCACTGGGGCTGGATA  >  NC_000913/4394432‑4394726                                                                                                                                                    |                                                                                                                                                      aaTATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGcc                                                                                                                                                    <  1:110122/151‑1 (MQ=255)      tGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATTCCTTTCCGcc                                                                                                                                                    >  1:426201/1‑146 (MQ=255)      tGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCct‑gct                                                                                                                                               <  2:983346/150‑1 (MQ=255)        tCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGGTATTCCCTCCCGcc                                                                                                                                                    >  2:493078/1‑144 (MQ=255)           ggCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑gctggctg                                                                                                                                          <  1:718440/150‑1 (MQ=255)                  aTCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGCTGTTTTTCCTTCCCGCCTGGCTGCCGGAAACTg                                                                                                                                    >  2:885614/1‑150 (MQ=255)                    cGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGGTATCCCCTCCGGCCC‑GGCGGCTGcaaaagg                                                                                                                                   >  2:585153/1‑144 (MQ=255)                    cGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGgcgc                                                                                                                                <  1:178372/151‑1 (MQ=255)                         gATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTTGTATTCCCTCCCGCCC‑TCCGGCTGaaaactgcgcctccc                                                                                                                           >  1:470796/1‑139 (MQ=255)                             tGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTTGTATTCCTTTCCGCCC‑GCTGGCTGCAAcgcgggcttcgc                                                                                                                           >  1:604223/1‑136 (MQ=255)                              gATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAg                                                                                                                        <  1:859599/149‑1 (MQ=255)                              gATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCCTCCGGCCC‑GCTTGCTGAAACTAGCGCCACGcc                                                                                                                          >  2:485082/1‑147 (MQ=255)                                ttgttgATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTTATATCCCCTCCCGCCC‑CCTGGCGGGAACCGgcgcctctccccgcg                                                                                                                     >  1:600693/1‑139 (MQ=255)                                 tgttgATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGgcgc                                                                                                                    <  1:22049/150‑1 (MQ=255)                                   ttgATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGGTATACCTTCCCGCCG‑CCTGGCTGGAACAGGCCCTGCGCCCGGCGCgg                                                                                                                  >  1:670711/1‑150 (MQ=255)                                     gATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTg                                                                                                                <  2:977717/150‑1 (MQ=255)                                        gcgcTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCa                                                                                                            <  1:958102/151‑1 (MQ=255)                                         cgcTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGTTATTCCTTCCCGCCC‑TCCGGCCGAAACTTGCCCTTCGCCCAGCGCGGGGGTCa                                                                                                            >  1:840111/1‑150 (MQ=255)                                         cgcTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCaa                                                                                                           <  1:833567/151‑1 (MQ=255)                                          gctgctCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAAc                                                                                                          <  2:424630/151‑1 (MQ=255)                                                 ggTACCGGTTTGCGGCGAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATc                                                                                                   <  1:697669/151‑1 (MQ=255)                                                    aCCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGGTATCCCCTCCGGCCC‑GCCGGCCGCAACTGGGGCCACGCCCCGCGCGGGGGTCAACCCCCATc                                                                                                   >  2:966919/1‑148 (MQ=255)                                                        gTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGTTATTCCTTCCCGCCC‑GCTCGCGGGAAAAGGCGCCACGCCCAGCGCGGGGGTCAACGCCgcaccca                                                                                                >  1:538604/1‑141 (MQ=255)                                                         tttGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTTATATTCCTTCCCGCCC‑GCTTGCTGAAACAGGCGCCTCGCCCGGCGCGGGGGTCAACGCCCAACACAcc                                                                                              >  1:101334/1‑148 (MQ=255)                                                               gCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATAACCCTTCCCGCCT‑GCTGCCGGAAACAGGGGCAACGCCCAGCGCGGGGGATAACGCCGAGCACAcccccac                                                                                         >  1:32997/1‑142 (MQ=255)                                                                 aaGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTAt                                                                                  <  2:293947/151‑1 (MQ=255)                                                                   gcgcCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTg                                                                                <  2:147129/151‑1 (MQ=255)                                                                           cgAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGaaa                                                                        <  1:546673/151‑1 (MQ=255)                                                                            gAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCCTCCCGCCC‑GCTCGCTGAAAACGGCGCCACGCCCGGGGCGGGGATCTACCCCGCTCACACCACCA‑CCTTTTTTTCGCTGaaa                                                                        >  2:591628/1‑150 (MQ=255)                                                                              aTCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATTCCCTCCCGCCC‑GCCGGCTGGAACCGGCGCCACGCCCCGCGCGGGGGTCAACGCCCACCACACCCTCA‑CCTTTTTTTCCCTGgaa                                                                        >  1:59277/1‑148 (MQ=255)                                                                              aTCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGATTGTAATCCTTCCCGCCT‑GCCCGGTGCAACATGGGCCACACCCGGCGCGGGGGTTAACCCCCGACACACCCTCA‑CTTTTTTTtcggcgtaacc                                                                      >  1:41036/1‑142 (MQ=255)                                                                               tCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGGTATTCCCTCCCGCCC‑GCTGGCTGGAAATGGCGCCACGCCCGGCGCGGGGGTCAAAGCCCAACCAACCCTCA‑CCTTTTTTTCGCTGaaaac                                                                      >  1:200155/1‑147 (MQ=255)                                                                                ccATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAg                                                                    <  1:762451/150‑1 (MQ=255)                                                                                 cATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGc                                                                  <  2:864857/151‑1 (MQ=255)                                                                                       ccAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTc                                                            <  1:330051/151‑1 (MQ=255)                                                                                         aGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCCTCCCGCCT‑GCCGGCTGAAACTGGCGCTAAGGCAGGCGCGGGGGTTAAAGCCGCGCACCCCCTCACCTTTTTTTGC‑CTGaaaaacagcgtggtcac                                                          >  2:313489/1‑137 (MQ=255)                                                                                             aaTCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTgg                                                       <  1:104239/150‑1 (MQ=255)                                                                                                 gACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTTCCGCCT‑GCTGGCTGAAACAGGGGCCACGCCCGGCGCGGGGGTCAACGCCGATCACACCATCA‑CTTTTTTTTTGCTGGAAACCGGCTTGCTCCCTGGgaaa                                                   >  1:737557/1‑150 (MQ=255)                                                                                                   ccACGCTAATTCCCATCCTGCGCCGATTGTGGCGGCTGATATCCCTTCCCGCCC‑GCTGGCCGCAACTAGGGCCACCCCAGGCGCGGGGGTCAAAGCCCGACACACCAACA‑CTTTTTTTTTGCCGAAAAAACGCCTGCTCACCGGgaaa                                                   >  2:947171/1‑148 (MQ=255)                                                                                                         tAATTCCCATCCTGCGCCGATTGTGGCGGGTGATTTCCCCTCCCCCCC‑CCTGGCTGGAACTGGCGCTCCGCCACGGGCGGGGGTCAACCCCCCTCACACCCTCC‑CCTTTTTTTCGCTGGAACCCCGCTTGCTCCCCGGGAAAACGCggg                                            >  2:266174/1‑149 (MQ=255)                                                                                                          aaTTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGAt                                          <  2:748273/150‑1 (MQ=255)                                                                                                             tCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGtt                                       <  1:839299/150‑1 (MQ=255)                                                                                                                cATCCTGCGCCGATTGTGGCGGTTTATATCCCCTCCCGCCC‑CCTTGCTGGCACTGGCGCCACGCCCAGCGCGGGGGTTAACGCCCCTCACACCCCCA‑CTTTTTTTTCGCTGTAACCAGGCGTGCTCACCGGAAAAACGCGGGGTGttttc                                    >  2:531668/1‑147 (MQ=255)                                                                                                                         ccGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGc                          <  1:128873/151‑1 (MQ=255)                                                                                                                                  gCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGCATTTTGACt                 >  2:207948/1‑151 (MQ=255)                                                                                                                                       tGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGCATTTTGACTCACTg            >  2:613922/1‑151 (MQ=255)                                                                                                                                           atCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGCATTTTGACTCACTGGGGc        >  2:355535/1‑151 (MQ=255)                                                                                                                                            tCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGCATTTTGACTCACTGGGGCt       >  2:818264/1‑151 (MQ=255)                                                                                                                                             cccTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGCATTTTGACTCACTGGGGCTg      >  1:568193/1‑151 (MQ=255)                                                                                                                                                  ccGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGCATTTTGACTCACTGGGGCTGGATa  <  2:252166/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AATATTGTCTGGCCCGAATCGGTAGATCTGATTGTTGATGCGCTGCTCGGTACCGGTTTGCGGCAAGCGCCCCGCGAATCCATTAGCCAGTTAATCGACCACGCTAATTCCCATCCTGCGCCGATTGTGGCGGTTGATATCCCTTCCGGCCT‑GCTGGCTGAAACTGGCGCTACGCCAGGCGCGGTGATCAACGCCGATCACACCATCA‑CTTTTATTGCGCTGAAACCAGGCTTGCTCACTGGAAAAGCGCGGGATGTTACCGGACAACTGCATTTTGACTCACTGGGGCTGGATA  >  NC_000913/4394432‑4394726

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