Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,404,2150TG30.0% 70.9 / 13.8 40G244G (GGT→GGG) yhdHacrylyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (1/27);  new base G (12/0);  total (13/27)
Fisher's exact test for biased strand distribution p-value = 2.33e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATT  >  NC_000913/3404079‑3404352
                                                                                                                                        |                                                                                                                                         
atgaGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGGCTGGCggggggg                                                                                                                                        >  2:14215/1‑139 (MQ=255)
   aGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGtttt                                                                                                                                     >  1:160406/1‑140 (MQ=255)
      ttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAct                                                                                                                                  <  1:410272/140‑1 (MQ=255)
      ttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAct                                                                                                                                  <  2:352917/140‑1 (MQ=255)
      ttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTct                                                                                                                                  >  2:410272/1‑140 (MQ=255)
          cGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGc                                                                                                                              <  1:823823/140‑1 (MQ=255)
                 cGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGAAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGGGGCCTGTGGTCTGGGGGGGGGGGTTTTTCTTCCAACCAc                                                                                                                       >  2:823823/1‑140 (MQ=255)
                   tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAGAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCAccg                                                                                                                     >  2:823896/1‑138 (MQ=255)
                   tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACgg                                                                                                                     <  1:263309/140‑1 (MQ=255)
                   tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTTCTCTGCCAACCACgg                                                                                                                     >  1:701870/1‑140 (MQ=255)
                     ctctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTc                                                                                                                   <  2:344509/140‑1 (MQ=255)
                       ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                 <  1:764667/140‑1 (MQ=255)
                        tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTTCTCTGCCCAACACCGTCATg                                                                                                                >  2:671531/1‑140 (MQ=255)
                        tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTTCCCACCCCCGGCATg                                                                                                                >  2:911663/1‑140 (MQ=255)
                         ggAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGGCTGGCGGGGGGGTTTTTTCTGCCAACCACCGGCATGc                                                                                                               >  2:808837/1‑140 (MQ=255)
                                cAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAt                                                                                                        <  1:911663/140‑1 (MQ=255)
                                    tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTg                                                                                                    <  1:219282/140‑1 (MQ=255)
                                          ctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTAATCTGCCAACCACGGGCATGCCATTTATTTTGCGGaaa                                                                                              >  1:448259/1‑139 (MQ=255)
                                             gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGGTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGGCATGCCATTTTTTCTGCGTAATGTc                                                                                           >  2:858135/1‑140 (MQ=255)
                                             gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGTCATGCCATTTATTCTTCGTAATGTc                                                                                           >  1:54754/1‑140 (MQ=255)
                                                  aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGttt                                                                                      >  1:328379/1‑140 (MQ=255)
                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                    <  2:160406/140‑1 (MQ=255)
                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                    <  1:14215/140‑1 (MQ=255)
                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                    <  2:11632/140‑1 (MQ=255)
                                                     gACACCGTTGGCGACAAAGTGCTGGTAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCa                                                                                   <  1:154443/140‑1 (MQ=255)
                                                      acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                  <  2:502761/140‑1 (MQ=255)
                                                      acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                  <  1:858135/140‑1 (MQ=255)
                                                      acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                  <  2:328379/140‑1 (MQ=255)
                                                       cacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGTCATGCCCTTTTTTTTGCGTAATGTCCGTTTGCAAg                                                                                 >  1:352683/1‑140 (MQ=255)
                                                           gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                             <  1:796904/140‑1 (MQ=255)
                                                            ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGATCTGGCGGGTGGTTTTTCTCTGCCAACCACGGGCATGCCATTTTTTCTGCGTAATGTCCGTTTTCAAgggggg                                                                            >  2:623788/1‑138 (MQ=255)
                                                             tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTgg                                                                           <  1:653214/140‑1 (MQ=255)
                                                                      aGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTTTCTGCCCACCACGGGCATGCCCTTTTTTTTTCGTAATGTCCGTTTTCAAGGGGGGGATTCAGTaa                                                                  >  2:449751/1‑140 (MQ=255)
                                                                       gTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGATCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAAt                                                                 <  1:623788/140‑1 (MQ=255)
                                                                            ggCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCCAACACGGGCATGCCCTTTTTTCTGCGGAATGTTCCTTTTCAAAGGGGGGATTCAGTAAAGACGc                                                            >  2:50387/1‑140 (MQ=255)
                                                                                      cTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAAcg                                                  <  1:93300/140‑1 (MQ=255)
                                                                                       tGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTTCCCACCCCCGGCATGCCATTTTTTTTGCGGAATGTCCGTTTGCCAGGGGGGGAGTTAGTAATGACGCCACCAGAAcgc                                                 >  2:25670/1‑140 (MQ=255)
                                                                                              aaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTTCCAACCCCGGGCATGCCATTTTTTTTGCGTAATGTCCGGTTGCAAGGGGGGGGTTCAGTAATGACGCCACCAGAACGCCGCGcac                                          >  2:776085/1‑140 (MQ=255)
                                                                                                             cTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGAc                           <  1:630081/140‑1 (MQ=255)
                                                                                                                   ggcggcCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCg                     <  1:264693/140‑1 (MQ=255)
                                                                                                                                 tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa       <  1:808837/140‑1 (MQ=255)
                                                                                                                                 tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa       <  1:823896/140‑1 (MQ=255)
                                                                                                                                 tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa       <  1:583332/140‑1 (MQ=255)
                                                                                                                                  ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  1:25670/140‑1 (MQ=255)
                                                                                                                                  ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  1:49317/140‑1 (MQ=255)
                                                                                                                                  ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  2:448259/140‑1 (MQ=255)
                                                                                                                                      ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt  >  1:355492/1‑140 (MQ=255)
                                                                                                                                      ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt  >  2:367322/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                         
ATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATT  >  NC_000913/3404079‑3404352

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