Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,256,0690GC37.0% 37.7 / 10.1 27P133R (CCG→CGG) plsBglycerol‑3‑phosphate 1‑O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (5/12);  new base C (10/0);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 3.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGG  >  NC_000913/4255827‑4256293
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                  
gACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCGGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACGGACGGCGAGAAACGCACAAAACGGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGcgcgcg                                                                                                                                                                                                                             <  1:38640/250‑1 (MQ=255)
                                aGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTTGCCTTTTTTCCCCCCC‑CGCGCGCGA‑CCCAACATCAACC‑ACACTTGCACCATCttccc                                                                                                                                                                                          >  1:342894/1‑247 (MQ=255)
                                     ttGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGCCCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATcc                                                                                                                                                                                      <  1:185313/250‑1 (MQ=255)
                                                           cccACGCGCTCCATTTTCTGGGCGATAGTTTTATCCGTTCCTTGTTCATCCGCCATACGGCTCAGCGAAACTGACGGCGAGAAACGCACAAACCTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGGTAAGCATACGCATCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCaa                                                                                                                                                               <  1:245134/250‑1 (MQ=255)
                                                                    gcCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAAAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATa                                                                                                                                                      <  1:300312/251‑1 (MQ=255)
                                                                              gAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAAc                                                                                                                                            <  1:100986/251‑1 (MQ=255)
                                                                                gCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCCAACATAACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGGTGCTACGGTGCCAAACGAGATAGTCGTgggac                                                                                                                                            >  1:81376/1‑246 (MQ=255)
                                                                                  gATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTC‑CTTTTCCCTTTTCTCTCCCC‑CGCGCCCGG‑CCAAAAAACCCCC‑ACACCGGCACCATCTTCACATCCAGCCTTGGGTTGGTGAGGTGCaaaaagaaaaaatagtgggacagcc                                                                                                                                        >  2:327642/1‑229 (MQ=255)
                                                                                       ttttATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCACTTCCCCTTTTTTACGCCCC‑CGCGCGCGC‑CCAAACATCACCC‑ACACTGGCACCATCTGCACATCCAGATTTTGGTTGCTACGGGGCAAAACGAGATAGTCGGGGAACAGCTTAAt                                                                                                                                    >  1:19734/1‑250 (MQ=255)
                                                                                           aTCCGTGCCGTTTCCACCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGAc                                                                                                                                <  2:331534/250‑1 (MQ=255)
                                                                                             ccgtgccgtgTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCACTTCGCCCTTTTTTCCCCCC‑GGCGCGCGAACCAAA‑A‑CACCCCACCCTGGCACCACCTGCACATCCACATTTTGGTTGTTACCGGGCCAAAAGAGCTAGACGTGGAACACCCTAATAGActct                                                                                                                             >  2:241246/1‑251 (MQ=255)
                                                                                                           cccGCCCTACGGCGCACCGAAACTTACGGCGCGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCCCTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCgt                                                                                                                  <  2:110568/247‑1 (MQ=255)
                                                                                                           cccGCCATACGGCGCAGCGAACCGGACGGCGAGAAACGCACAAAACTGTCGCGCCCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTa                                                                                                               <  1:83107/250‑1 (MQ=255)
                                                                                                              gCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACCCCCC‑CGCGCGCGA‑CCAAAAATCACCC‑ACACTGGCCCCATCTGCACATCCCGATTTGGGTTGCTACCGTGCAAATCGAGAAAGTCGTGGAACAGCTTAATAGACTCTTCTTTCCGCGTGTaaa                                                                                                             >  1:275805/1‑249 (MQ=255)
                                                                                                                       cgcAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTAACTTCCACTTTTTTAAGCCCC‑CGCGCGCGA‑CCAAACATCAACG‑ACACTGGGAACCCCTGAAAAAACAGATTTGGGGTGCTAAGGGGGAAATAGAGATAGTGGGGGGAAAGCTTAATAGACTCTTTTTTTGGCGTGTAATAAGTGAaca                                                                                                    >  1:232073/1‑250 (MQ=255)
                                                                                                                                     cGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCACTTCGCCTTTTTTACCCCCC‑CGCCCGCGA‑CCAAACATCACCG‑ACCCTGGCACCCTCTGCACCTCCCCATTTTGGTTGCTACGGTGCCAAACGAGATAGTCGTGGGACAACTTAATAAACTCTTCTTTCTGCGTGTAATAGGTGAACACACCCGGCccgcccc                                                                                     >  2:74223/1‑249 (MQ=255)
                                                                                                                                            aaaCGCAAAAAACGGCCGCGCCCG‑GCCCCAGTACAAGAGAAAATTTTTTGTACGCCTTTAAGCATACGCGGGGGCGGG‑T‑CCCTTTCCCTTTTTCCCCCCCCCGGCGCACGC‑AAAAACAGCCCCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATATTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAAGAGGTGAACACCCGCGCCCCGCCGtgaatgaa                                                                              <  2:224119/251‑1 (MQ=255)
                                                                                                                                                        cTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTTGCCCTTTTCACACCCC‑CGCGCGCGA‑CCAAAAAATACCG‑AAAACGGGACCCCATTCACATCCAGATTTTGGGTGCTACGGGGCAAATCGAGAAAGGCGTGGGAAAACTTAATAAACTCTTTTTTCTGCGTGGAATAAGTGAACAAACACGGCCCGCCGtgaatgaaacaataacgcgg                                                                  >  1:94077/1‑239 (MQ=255)
                                                                                                                                                                aCCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCACCTTTCCTTTTTTACGCCCC‑CGCGCGCGA‑CCCAACATCACCC‑ACACTGGCCCCATCTCCACATCCAGATTTGGGGTGCTACGGGGCAAAACGAGATAGGCGTGGGACAGCATTATAGACTCTTCTTTCGGCGGGGAATAGGCGAACACACCCCGCCCCCCGCGAATGAaccactagcggggcaggagac                                                          >  2:67790/1‑232 (MQ=255)
                                                                                                                                                                       cacaGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTTACGCCCC‑GGCGCGCGA‑CCAAACATCACCC‑ACACTGGCCCCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAAAAGGTGAAAACACGCGGCCCGCCGTGAATGAACACATAACGCGGCAGAAGCGTGCCGTc                                                   >  1:97039/1‑251 (MQ=255)
                                                                                                                                                                       cacaGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑GTCACTTCGCCTTTTTTACGCCCC‑CGCGCGCGA‑CCAAACATCACCG‑ACACTGGCCCCCTCTGCACATCCAGATTTTGGGTGCTACGGGGCAAATCGAGATAGTCGTGGAAAAGCTTAATAGAATCTTCTTTTGGCGTGTAATAGGGGAACACACCCGGCCCCCCGTGAATGAAAAAAAAACGCGGCAGTAACGTGcccgc                                                   >  1:91762/1‑248 (MQ=255)
                                                                                                                                                                       cacaGTACA‑GAAAAAAATTTTTGTCCGCAGTAAAGCATCCGCGGGGGGGGT‑TTCCTTTCGCTTTTTCCCCCCCCC‑GGCGCGCGA‑CAAAACACCACCG‑ACACGGGCACCATCTGCACATCCAGATTGGGTTTGCTACGGTGCAAAGCGAGATAGTCGTGGAACAGCTTAATAGCCTCTTCTTTCGGCGTGTAATAGGTGAACACGCGGGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTc                                                   <  1:151523/251‑1 (MQ=255)
                                                                                                                                                                       caaagtccA‑GCAAAAAATTTCTGTCCGCCTTTAAGAATCCGGAGGGGAGGG‑TTCCCTTCCCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGTTACGTTGCAATCCGAGATGGTCGGGGAACAGCTTAAAAGCCTCTTCTTTCGGCGTGTAATAGGTGAAAACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTc                                                   <  1:324812/244‑1 (MQ=255)
                                                                                                                                                                                   gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCCTTTTTACGCCCC‑CGCGCGCGC‑CCCAACATCACCC‑ACACTGGCACCCTCTGCACCTCCAGATTTTGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCCGCGGGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCg                                        >  1:294654/1‑251 (MQ=255)
                                                                                                                                                                                      aaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCCCCTTTTTTACGCCCC‑GGCGCGCGA‑CCAAACATCACCC‑ACACCGGCACCATCTGCACATCCAGATTTTGGTTGCTACCGTGCAAAACGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGGGAATGAAAACATAGCGCGGCCGTAGCGGGCCGGCGATTTCCAACGGcc                                     >  1:236014/1‑250 (MQ=255)
                                                                                                                                                                                              gtccGCCGTAAAGCATACGCAGCGGCGGT‑TTCACTTCGCTTTTTCCCCCCCCC‑GGCGCGCGA‑CAAAACATCACCG‑ACACGGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAAc                                <  1:74223/245‑1 (MQ=255)
                                                                                                                                                                                                                         gg‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTgg  >  1:190630/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                  
GACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACA‑GCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGG‑TTCACTTCGCCTTTTTCACGCCCC‑GGCGCGCGA‑CCAAACATCACCG‑ACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGG  >  NC_000913/4255827‑4256293

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