Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,4711.G51.7% 43.2 / 48.0 29coding (398/1113 nt)ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (5/9);  new base G (8/7);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-02
Rejected as consensus: Frequency below/above cutoff threshold.

TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGC  >  NC_000913/1769339‑1769597
                                                                                                                                     |                                                                                                                               
tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑gggg                                                                                                                             >  1:1221528/1‑137 (MQ=255)
tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑gggg                                                                                                                             <  2:1221528/137‑1 (MQ=255)
                      aaTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTc                                                                                                    <  2:473023/139‑1 (MQ=255)
                             gCCATTCCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGCATGCAGGCTGGAACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCt                                                                                              <  1:1006916/139‑1 (MQ=255)
                                 tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGt                                                                                                     <  2:1605482/127‑1 (MQ=255)
                                 tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGt                                                                                                     >  1:1605482/1‑127 (MQ=255)
                                     cAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGc                                                                                     <  1:479030/139‑1 (MQ=255)
                                      aGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGca                                                                                    <  1:611545/139‑1 (MQ=255)
                                             aCATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacca                                                                              >  2:730677/1‑139 (MQ=255)
                                                               tGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGc                                                            <  1:153387/139‑1 (MQ=255)
                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTc                                                                                                    >  1:337463/1‑96 (MQ=255)
                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                         >  1:1254551/1‑139 (MQ=255)
                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACCTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTc                                                                                                    <  2:337463/96‑1 (MQ=255)
                                                                              ccATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCgg                                            <  1:679092/139‑1 (MQ=255)
                                                                                 ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgc                                          >  1:299129/1‑139 (MQ=255)
                                                                                  tCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgcg                                         <  2:1254551/139‑1 (MQ=255)
                                                                                            ggCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTg                               >  1:1259763/1‑139 (MQ=255)
                                                                                                               gCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacc                                                                            >  1:700172/1‑75 (MQ=255)
                                                                                                               gCTGGCACAACTTGCTGGATAT‑GGGGGGGTACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacc                                                                            <  2:700172/75‑1 (MQ=255)
                                                                                                                cTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtctc          >  1:558363/1‑139 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGt                                <  1:1452875/116‑1 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGt                                >  2:1452875/1‑116 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCtt        >  2:294227/1‑139 (MQ=255)
                                                                                                                        aCTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGc  <  1:294227/139‑1 (MQ=255)
                                                                                                                         cTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                       >  2:1126011/1‑102 (MQ=255)
                                                                                                                         cTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                       <  1:1126011/102‑1 (MQ=255)
                                                                                                                         cTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGc  <  2:299129/139‑1 (MQ=255)
                                                                                                                         cTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTg                                                                                       <  1:485707/53‑1 (MQ=255)
                                                                                                                         cTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTg                                                                                       >  2:485707/1‑53 (MQ=255)
                                                                                                                             cTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg   >  2:187783/1‑133 (MQ=255)
                                                                                                                             cTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg   <  1:187783/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGC  >  NC_000913/1769339‑1769597

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