Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,404,2180TG32.6% 70.1 / 15.0 43G245G (GGT→GGGyhdHacrylyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (4/25);  new base G (14/0);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 3.90e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATT  >  NC_000913/3404083‑3404352
                                                                                                                                       |                                                                                                                                      
gTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTa                                                                                                                                    <  2:521129/140‑1 (MQ=255)
   tGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTtctc                                                                                                                                 >  2:621324/1‑140 (MQ=255)
     ccGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTg                                                                                                                               <  1:821163/140‑1 (MQ=255)
       gAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGcc                                                                                                                             <  2:144437/140‑1 (MQ=255)
          tCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGGCTGGCGGGGGGGTTTTTTCTGCCAAc                                                                                                                          >  2:31524/1‑140 (MQ=255)
               tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACgg                                                                                                                     >  1:775105/1‑140 (MQ=255)
                 ctctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTc                                                                                                                   <  1:666850/140‑1 (MQ=255)
                   ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                 <  2:546918/140‑1 (MQ=255)
                   ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                 <  1:31524/140‑1 (MQ=255)
                       aaaaaCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCa                                                                                                             <  2:377105/140‑1 (MQ=255)
                            cAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAt                                                                                                        >  1:826310/1‑140 (MQ=255)
                            cAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGTCATGCCAttttt                                                                                                        >  1:675890/1‑138 (MQ=255)
                                tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTg                                                                                                    <  1:621324/140‑1 (MQ=255)
                                tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTg                                                                                                    <  1:197039/140‑1 (MQ=255)
                                        ggggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGt                                                                                            >  1:308412/1‑140 (MQ=255)
                                          ggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTcc                                                                                          <  2:357025/140‑1 (MQ=255)
                                                tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                    <  1:886074/140‑1 (MQ=255)
                                                tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                    <  1:730276/140‑1 (MQ=255)
                                                   cacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTTTTCTGCGGTATGTCCGTTTGGAAg                                                                                 >  2:568135/1‑140 (MQ=255)
                                                       gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                             <  1:568135/140‑1 (MQ=255)
                                                       gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                             <  1:642996/140‑1 (MQ=255)
                                                         tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTgg                                                                           >  1:729631/1‑140 (MQ=255)
                                                         tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCCACCACGGGCATGCCATTTATTTTGCGTAATGTCCGTTTTCAAGgggggg                                                                           >  2:174262/1‑140 (MQ=255)
                                                         tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCCCGGGCATGCCCTTTTTTTTGCGTAATGTCCGTTTGCAAGgggggg                                                                           >  2:730276/1‑140 (MQ=255)
                                                         tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCCCGGGCATGCCCTTTTTTCTGCGGAATGTCCGTTTTCAAGgggggg                                                                           >  1:935165/1‑140 (MQ=255)
                                                         tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTTCCAACCACGGTCATGCCATTTATTCTGCGTAATGGCCGTTTGCAAGgggggg                                                                           >  1:721155/1‑140 (MQ=255)
                                                          ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGa                                                                          <  1:797374/140‑1 (MQ=255)
                                                          ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGa                                                                          <  1:756524/140‑1 (MQ=255)
                                                          ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGa                                                                          <  1:105144/140‑1 (MQ=255)
                                                          ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGa                                                                          <  1:174262/140‑1 (MQ=255)
                                                          ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTATTCTGCGTAATGTCCGGTTGCAAGggggggg                                                                          >  2:850340/1‑139 (MQ=255)
                                                           gCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTATTCTGCGGAATGTCCGTTTTGAAGgggggggt                                                                         >  2:756524/1‑138 (MQ=255)
                                                           gCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTATCTGCCCACCCCGGGCCTGCCATTTTTTTTGCGTAATGTCCGTTTGGCAGgggggggt                                                                         >  2:48143/1‑138 (MQ=255)
                                                            cGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTACCAACCACGGGCATGCCCTTTTTTTTGCGTAATGTCCGTTTTCAAGGGGGGGAtt                                                                        >  2:732739/1‑140 (MQ=255)
                                                                              aGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCCACCCCGGGCCTGCCCATTTTTCTGCGTAATGTCCGGTTTCAAGGGGGGGGTTCAGTAAAGACGCCACCAg                                                      >  2:786130/1‑140 (MQ=255)
                                                                                   tGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCCACCACGGGCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTTAGTAATGACGCCACCAGAAcgc                                                 >  1:721408/1‑140 (MQ=255)
                                                                                    ggCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAAcgcc                                                <  2:826310/140‑1 (MQ=255)
                                                                                        cAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCgcgc                                            <  1:550468/140‑1 (MQ=255)
                                                                                          aaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGcac                                          <  1:850340/140‑1 (MQ=255)
                                                                                                   cggcggCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTTTTTTGCGTAATGTCCGTTTTCAAGGGGGGGGTTTAGTAATGACGCCACCAGAACGCCGCGCCCAAACCCGGc                                 >  2:550468/1‑140 (MQ=255)
                                                                                                     gcggcTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTTCCCACCACCGGCATTCCCTTTATTCTGCGGAATGTCCGTTTTCAAGGGGGGGATTTAGTAATGACGCCCCCCGAACGCCCCGCCCAAACCTTGccg                               >  2:686291/1‑138 (MQ=255)
                                                                                                        gcTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTTCTCTGCCAACCACGGGCATGCCCTTTTTTTTGCGTAATGTCCGTTTGGAAGGGGGGGGTTTAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGa                            >  2:174643/1‑140 (MQ=255)
                                                                                                        gcTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTTATCTGCCAAACCCGGGCATGCCCTTTATTCTGGGGAAAGTCCGTTTGGAAGGGGGGGATTTAGTAATGACGCCCCCAGAAAGCCGCGCACAAGCCTGGCAGCgg                            >  2:835407/1‑139 (MQ=255)
                                                                                                                    ccTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGAt                <  1:441788/140‑1 (MQ=255)
                                                                                                                     cTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGAtt               <  1:933502/140‑1 (MQ=255)
                                                                                                                              ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      >  2:44708/1‑140 (MQ=255)
                                                                                                                              ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  2:675890/140‑1 (MQ=255)
                                                                                                                              ggCGGGTGGTTTTACTCTACCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  1:732739/140‑1 (MQ=255)
                                                                                                                                  ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt  <  1:44708/140‑1 (MQ=255)
                                                                                                                                  ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt  <  1:140355/140‑1 (MQ=255)
                                                                                                                                  ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt  >  1:122617/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATT  >  NC_000913/3404083‑3404352

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