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.

TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGC  >  NC_000913/1769339‑1769597
                                                                                                                                     |                                                                                                                               
tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATgggg                                                                                                                             >  1:1221528/1‑137 (MQ=255)
tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATgggg                                                                                                                             <  2:1221528/137‑1 (MQ=255)
                      aaTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                    <  2:473023/139‑1 (MQ=255)
                             gCCATTCCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGCATGCAGGCTGGAACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCt                                                                                              <  1:1006916/139‑1 (MQ=255)
                                 tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGt                                                                                                     <  2:1605482/127‑1 (MQ=255)
                                 tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGt                                                                                                     >  1:1605482/1‑127 (MQ=255)
                                     cAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGc                                                                                     <  1:479030/139‑1 (MQ=255)
                                      aGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGca                                                                                    <  1:611545/139‑1 (MQ=255)
                                             aCATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacca                                                                              >  2:730677/1‑139 (MQ=255)
                                                               tGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGc                                                            <  1:153387/139‑1 (MQ=255)
                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                    >  1:337463/1‑96 (MQ=255)
                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                         >  1:1254551/1‑139 (MQ=255)
                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACCTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                    <  2:337463/96‑1 (MQ=255)
                                                                              ccATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCgg                                            <  1:679092/139‑1 (MQ=255)
                                                                                 ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgc                                          >  1:299129/1‑139 (MQ=255)
                                                                                  tCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgcg                                         <  2:1254551/139‑1 (MQ=255)
                                                                                            ggCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTg                               >  1:1259763/1‑139 (MQ=255)
                                                                                                               gCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacc                                                                            >  1:700172/1‑75 (MQ=255)
                                                                                                               gCTGGCACAACTTGCTGGATATGGGGGGGTACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacc                                                                            <  2:700172/75‑1 (MQ=255)
                                                                                                                cTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtctc          >  1:558363/1‑139 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGt                                <  1:1452875/116‑1 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGt                                >  2:1452875/1‑116 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCtt        >  2:294227/1‑139 (MQ=255)
                                                                                                                        aCTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGc  <  1:294227/139‑1 (MQ=255)
                                                                                                                         cTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                       >  2:1126011/1‑102 (MQ=255)
                                                                                                                         cTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                       <  1:1126011/102‑1 (MQ=255)
                                                                                                                         cTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGc  <  2:299129/139‑1 (MQ=255)
                                                                                                                         cTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTg                                                                                       <  1:485707/53‑1 (MQ=255)
                                                                                                                         cTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTg                                                                                       >  2:485707/1‑53 (MQ=255)
                                                                                                                             cTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg   >  2:187783/1‑133 (MQ=255)
                                                                                                                             cTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg   <  1:187783/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGC  >  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: