Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913691,3900GC23.3% 57.6 / 16.2 30A132G (GCC→GGC) ybeXputative ion transport
Reads supporting (aligned to +/- strand):  ref base G (13/10);  new base C (3/4);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGT  >  NC_000913/691249‑691531
                                                                                                                                             |                                                                                                                                             
aCCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAg                                                                                                                                         >  4:63570/1‑148 (MQ=255)
aCCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCCACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAg                                                                                                                                         >  5:126463/1‑148 (MQ=255)
  cAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCAGGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAt                                                                                                                                      <  8:144744/149‑1 (MQ=255)
       cGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTg                                                                                                                                 >  5:91453/1‑149 (MQ=255)
                cACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTATCACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAAc                                                                                                                        <  6:91453/149‑1 (MQ=255)
                  ccGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCCACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACgg                                                                                                                      <  6:126463/149‑1 (MQ=255)
                               tAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgga                                                                                                                                   >  6:48/1‑120 (MQ=255)
                               tAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgga                                                                                                                                   <  5:48/123‑4 (MQ=255)
                                aaCGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCg                                                                                                        >  4:99008/1‑149 (MQ=255)
                                    aTCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCAt                                                                                                    <  3:99008/149‑1 (MQ=255)
                                        ccATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACg                                                                                                                       >  2:109989/1‑126 (MQ=255)
                                        ccATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACg                                                                                                                       <  1:109989/126‑1 (MQ=255)
                                               gTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgga                                                                                                                                   <  6:21456/107‑4 (MQ=255)
                                               gTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgga                                                                                                                                   >  5:21456/1‑104 (MQ=255)
                                                     cTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgg                                                                                                                                    <  1:41419/100‑3 (MQ=255)
                                                     cTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgg                                                                                                                                    >  2:41419/1‑98 (MQ=255)
                                                                 ctctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCtt                                                                       <  3:63570/149‑1 (MQ=255)
                                                                    tttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                    >  3:29096/1‑149 (MQ=255)
                                                                        aGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGAGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAAt                                                                <  6:24129/149‑1 (MQ=255)
                                                                                             cTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGc                                           >  5:155261/1‑149 (MQ=255)
                                                                                              tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCaa                                                                                      >  5:119499/1‑105 (MQ=255)
                                                                                              tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCaa                                                                                      <  6:119499/105‑1 (MQ=255)
                                                                                                tCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGa                                                        >  4:52340/1‑133 (MQ=255)
                                                                                                tCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGa                                                        <  3:52340/133‑1 (MQ=255)
                                                                                                    gAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTc                                     >  2:54650/1‑148 (MQ=255)
                                                                                                                 cTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCaa                       <  3:86750/149‑1 (MQ=255)
                                                                                                                            cTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTcc            <  6:155261/149‑1 (MQ=255)
                                                                                                                                gTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGTAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCg        >  6:161988/1‑149 (MQ=255)
                                                                                                                                 tCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCCAGACATTCGTCCAGCGt       >  6:28774/1‑149 (MQ=255)
                                                                                                                                      gCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGt  >  1:180145/1‑149 (MQ=255)
                                                                                                                                             gCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGt                                               >  2:120602/1‑97 (MQ=255)
                                                                                                                                             gCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGt                                               <  1:120602/97‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
ACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGT  >  NC_000913/691249‑691531

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑