Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,926,366	0	A	C	32.6%	73.8 / 14.6	43	G30G (GGT→GGG)	mioC	flavoprotein MioC
Reads supporting (aligned to +/- strand):  ref base A (24/5);  new base C (0/14);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.48e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAA  >  NC_000913/3926231‑3926498                                                                                                                                        |                                                                                                                                     ttCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCg                                                                                                                                  <  1:300865/140‑1 (MQ=255)            cTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGc                                                                                                                       >  1:238671/1‑140 (MQ=255)             ttCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCt                                                                                                                      >  1:262032/1‑140 (MQ=255)                aTAGAAAGGAGAATGGTTGCCGGGAATTTCTCCGGCCCCGTGGGTGGGGGTGTTAACCAGCAACCTCCCGGAGGCAGTAAAATCTTCTAACAGCGGCCCGTGCGGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCtttt                                                                                                                   <  2:238671/140‑1 (MQ=255)                     aaggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGACCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGcc                                                                                                              <  2:164062/140‑1 (MQ=255)                         agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCa                                                                                                                 >  2:558210/1‑133 (MQ=255)                         agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCa                                                                                                                 <  1:558210/133‑1 (MQ=255)                           aaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgtgt                                                                                                        >  1:239703/1‑140 (MQ=255)                            aaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAACCAGGTGtt                                                                                                       <  1:448289/140‑1 (MQ=255)                                  gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCt                                                                                                  >  2:531916/1‑139 (MQ=255)                                  gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCt                                                                                                  <  1:531916/139‑1 (MQ=255)                                  gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                 >  1:58547/1‑140 (MQ=255)                                                              gAGCTGATACCCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGcc                                                                     <  2:377481/140‑1 (MQ=255)                                                                  tGATACCCAGCCAGATCCCGGAGGCAGTAAAATTTTCTAACAGGGGCCCGGGCAGTGTTTGGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTg                                                                 <  2:262032/140‑1 (MQ=255)                                                                    aTACCCAGCCAAACCCCGGGGGCAGGTAAATCTTCTAACAGCGGCCCGTGCGGCTTTTCGGTGTAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAt                                                               <  2:259049/140‑1 (MQ=255)                                                                    aTAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAt                                                               >  1:53840/1‑140 (MQ=255)                                                                     agacaggccagATCCCTGAGGCAGGAAAATCTTCTAACAGCGGACCGTGCAGGTTTTCGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGTCGCTGATa                                                              <  2:300809/134‑1 (MQ=255)                                                                     aaaCCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCTTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATa                                                              <  2:198603/139‑1 (MQ=255)                                                                               gaccccGGAGGCAGGTAAATCTTCAAACAGCGACCCGTGCGGCTTTTCGGGGGAAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGatat                                                    <  1:101251/137‑1 (MQ=255)                                                                               gaccccGGAGGAGGGAAAATCTTCTAACAGGGGCCGGGGCAGCTTTTGGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGatat                                                    <  2:58547/137‑1 (MQ=255)                                                                                      gAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGTCGCTGATAAGAGTGATATcctgtct                                             >  1:300809/1‑134 (MQ=255)                                                                                       aGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATa                                            >  2:215409/1‑140 (MQ=255)                                                                                             taaaaCCTTCAAACAGCGGACCGGGCAGCGTTTCGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                      <  1:11667/138‑1 (MQ=255)                                                                                               aatttttcaAACAGCGGCCCGGGCGGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa                                    <  2:218525/135‑1 (MQ=255)                                                                                                     tctAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTtta                              >  2:72108/1‑140 (MQ=255)                                                                                                       tAACAGGGGCCGGTGCAGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatt                            <  1:370722/140‑1 (MQ=255)                                                                                                       tAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGCGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatt                            >  2:128336/1‑140 (MQ=255)                                                                                                       tAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatt                            >  1:218525/1‑140 (MQ=255)                                                                                                       tAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatt                            >  2:505369/1‑140 (MQ=255)                                                                                                           ggCGGCCCGTGCAGGTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTCATTAaga                        <  2:464869/139‑1 (MQ=255)                                                                                                            gCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAagag                       >  2:121944/1‑140 (MQ=255)                                                                                                             cGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGt                      >  1:341258/1‑140 (MQ=255)                                                                                                              ggccggTGCGGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg                     <  2:228862/135‑1 (MQ=255)                                                                                                              ggACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttattactgtct                     >  1:188565/1‑134 (MQ=255)                                                                                                               gACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTgg                    >  1:228862/1‑140 (MQ=255)                                                                                                                aCCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGc                   >  2:173753/1‑140 (MQ=255)                                                                                                                    tGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTAt               >  2:546295/1‑140 (MQ=255)                                                                                                                    tGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTAt               >  2:448646/1‑140 (MQ=255)                                                                                                                    tGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTAt               >  2:481725/1‑140 (MQ=255)                                                                                                                          gTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACg         <  1:173753/140‑1 (MQ=255)                                                                                                                          gTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACg         <  1:546295/140‑1 (MQ=255)                                                                                                                          gTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACg         >  2:31856/1‑140 (MQ=255)                                                                                                                            ttCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCt       >  1:241300/1‑140 (MQ=255)                                                                                                                            ttCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCt       >  1:77371/1‑140 (MQ=255)                                                                                                                                 tggtAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGaa  <  1:72108/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     TTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAA  >  NC_000913/3926231‑3926498

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