Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,926,3660AC28.6% 77.9 / 11.4 42G30G (GGT→GGGmioCflavoprotein MioC
Reads supporting (aligned to +/- strand):  ref base A (24/6);  new base C (0/12);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.68e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACG  >  NC_000913/3926227‑3926500
                                                                                                                                            |                                                                                                                                       
tCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGCCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTaaaa                                                                                                                                         <  2:263419/140‑1 (MQ=255)
         gCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCCCCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCctct                                                                                                                               <  2:400535/140‑1 (MQ=255)
            aTGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCCCCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGCCCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTcc                                                                                                                            <  2:278827/140‑1 (MQ=255)
            aTGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTcc                                                                                                                            >  2:212677/1‑140 (MQ=255)
            aTGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTcc                                                                                                                            >  2:332661/1‑140 (MQ=255)
              gCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAg                                                                                                                          >  2:163115/1‑140 (MQ=255)
                 taaTAGAAAGGAGAAAGGTTGTCCGAAATATCCCGGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACGGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCtt                                                                                                                       <  2:403813/138‑1 (MQ=255)
                 tCAAAGAAGAGAGAAAGGTTGTCCGAAATACCCCCGCCACCGTGGGTGGAGCTGATACCCACCCAGATCCCTGAGCCAGGTAAATCTTCAAACAGCGGACGGTGCGGCGTTCGGGTGGTAAACCCCGCCTCTTCCAGCtt                                                                                                                       <  2:550424/140‑1 (MQ=255)
                       aaaggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGc                                                                                                                 >  1:734573/1‑140 (MQ=255)
                        aaggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGcc                                                                                                                <  1:319537/140‑1 (MQ=255)
                          ggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAg                                                                                                              <  1:163115/140‑1 (MQ=255)
                            agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAggg                                                                                                            >  2:396482/1‑139 (MQ=255)
                            agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgt                                                                                                            >  2:155269/1‑140 (MQ=255)
                            agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgt                                                                                                            >  2:690949/1‑140 (MQ=255)
                             gaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgtg                                                                                                           <  2:507473/140‑1 (MQ=255)
                              aaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgtgt                                                                                                          >  1:364023/1‑140 (MQ=255)
                               aaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGtt                                                                                                         <  1:690949/140‑1 (MQ=255)
                                 gttggTCCGGAATATGTCACCCCCCGTGGGGGGGGGGACTACCCAGCCAGTTCACTGACGCAGGAACATCTTCAAACAGCGGGCAGCCCCGCGTTGCGGTGGAAAGCCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCt                                                                                                       <  2:568041/136‑1 (MQ=255)
                                     gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                   >  1:177958/1‑140 (MQ=255)
                                      tCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGGAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCGTCCAGCTTTTCAGCCAGGTGTTCTGCTAc                                                                                                  >  2:567276/1‑140 (MQ=255)
                                        cGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACat                                                                                                >  2:682673/1‑140 (MQ=255)
                                         gaaaTATCCCCGGCCCCGGGGGTGGAGCTGATACCCAGCCAGACCCCGGGGGAAGGAAAATCTTCAAACAGCGGCCGGTGCGGCTTTTGGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACata                                                                                               <  2:494308/138‑1 (MQ=255)
                                                                 gggCTAATACCCACCCAGCCCCCGGAGGCAGAAAATTCTTCAAACGGGGGCCGGTGCGGTTTTTGGGTGAAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGcc                                                                       <  2:734573/138‑1 (MQ=255)
                                                                   gCTAATACCCAGCCAGATCCCTGGGGCAGGTAAATCTTCAAACGGCGGACGGTGCAGCGTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGc                                                                     <  2:290337/140‑1 (MQ=255)
                                                                              gccagATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTg                                                          >  1:755738/1‑140 (MQ=255)
                                                                                 agATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGata                                                       >  2:656190/1‑140 (MQ=255)
                                                                                   aTCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGatact                                                     >  1:777281/1‑138 (MQ=255)
                                                                                      ccTGAGGCAGGAAAATCTTCTAACAGCGGCCCGTGCAGCTTTTCGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGc                                                  <  1:567276/140‑1 (MQ=255)
                                                                                         gAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCAt                                               >  1:31928/1‑140 (MQ=255)
                                                                                         gAGGCAGGTAAATCCTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCAt                                               >  1:280448/1‑140 (MQ=255)
                                                                                         gAGGCAGGAAAATCTTCAAACAGCGGACCGTGCAGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCAt                                               <  1:159613/140‑1 (MQ=255)
                                                                                                gaaaaTCTTCTAACAGGGGCCCGTGCGGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                        <  2:280448/138‑1 (MQ=255)
                                                                                                gTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                        >  1:50209/1‑140 (MQ=255)
                                                                                                  aatttttcaAACAGCGGCCGGTGCAGCTTTTGGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa                                      <  2:364023/135‑1 (MQ=255)
                                                                                                  aaaTTTTCAAACAGCGGCCCGGGCAGGTTTCGGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa                                      <  1:458026/140‑1 (MQ=255)
                                                                                                  aaaTCTTCTAACAGGGGCCCGTGCAGCGTTTGGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa                                      <  1:261508/140‑1 (MQ=255)
                                                                                                     tcttcaAACAGGGGCCGGTCCGGCTTTTGGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCt                                   <  2:50209/140‑1 (MQ=255)
                                                                                                          tAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatt                              >  2:36572/1‑140 (MQ=255)
                                                                                                                cgggcagTGGAGCGTTTGGGGGTAAAACCCCGCCTCTTCCAGTTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGt                        <  2:755738/134‑1 (MQ=255)
                                                                                                                 ggcccGTGCAGCTTTTCGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg                       <  1:396482/137‑1 (MQ=255)
                                                                                                                 ggACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg                       >  1:115948/1‑140 (MQ=255)
                                                                                                                 ggACCGTGCAGCGTTTCGCTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGGGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg                       >  1:458473/1‑140 (MQ=255)
                                                                                                                   cccGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGc                     <  1:598262/139‑1 (MQ=255)
                                                                                                                       tGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTAt                 >  1:104807/1‑140 (MQ=255)
                                                                                                                       tGCAGCGTTTCGGCGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTAt                 >  2:285469/1‑140 (MQ=255)
                                                                                                                              tttCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGc          >  1:305971/1‑140 (MQ=255)
                                                                                                                                     gggaaaacCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTCCATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGaa    <  2:618927/136‑1 (MQ=255)
                                                                                                                                       gtAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAAcg  >  2:375577/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                       
TCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACG  >  NC_000913/3926227‑3926500

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: