Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3350TG53.1% 33.8 / 46.3 32G25G (GGT→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (3/12);  new base G (17/0);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 2.02e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-05

TGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCG  >  NC_000913/523193‑523467
                                                                                                                                              |                                                                                                                                     
tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGaa                                                                                                                                >  2:266766/1‑150 (MQ=255)
     gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTTgcg                                                                                                                          >  2:2197170/1‑151 (MQ=255)
          tttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGccc                                                                                                                     >  2:3801536/1‑151 (MQ=255)
            tattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGccccc                                                                                                                   >  1:2617270/1‑149 (MQ=255)
                            tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTc                                                                                                    >  2:2887057/1‑150 (MQ=255)
                              acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgt                                                                                                 <  1:255137/151‑1 (MQ=255)
                                 caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGGTCGgggtg                                                                                              >  2:1626029/1‑151 (MQ=255)
                                      atatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCcgg                                                                                         >  2:529013/1‑151 (MQ=255)
                                       tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCGGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCcggg                                                                                        >  2:307448/1‑151 (MQ=255)
                                         tGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCAGGGTTCGGCAGGGGTAAGAATTGGGGCGCCCACCGGCGTGGCGTGCCCGGGGGGGCcggg                                                                                        >  1:3546699/1‑149 (MQ=255)
                                          gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCGGTGTGTCcgggcgg                                                                                     >  2:1153095/1‑151 (MQ=255)
                                           aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTTCAGGGTTCGGCAGGGGGAAGAATTGGCGGGCCCACCGGCGTGGCGTGGCCCGTGTGTAcggggggg                                                                                    >  2:1185365/1‑151 (MQ=255)
                                           aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcggg                                                                                    <  1:2606203/151‑1 (MQ=255)
                                                 aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACt                                                                              <  2:749663/151‑1 (MQ=255)
                                                     aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCgg                                                                          <  2:1108097/151‑1 (MQ=255)
                                                       aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGc                                                                        <  1:2790707/151‑1 (MQ=255)
                                                             aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCg                                                                 <  2:3486044/151‑1 (MQ=255)
                                                                 tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCACCCGCGTGGCGTGCTCCGGGTGGCcgggcggggagaatttcgg                                                                         >  1:2268009/1‑129 (MQ=255)
                                                                         cGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAAAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGGGTCCGGGCGGGGTGACTTCGGGC‑AACCCCGTAAATcccc                                                       >  1:174058/1‑147 (MQ=255)
                                                                                 cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCACGGTTCTGACGGGGTAAGAATTTGCGCGCCCCCCGGCGTTGCGTGCACCGTGTGTTCCGGGGGGGTGACAGCGGGC‑AACCCCGTAAAAccccag                                                     >  2:3387320/1‑139 (MQ=255)
                                                                                         cgTCAGGGAGATATGACTCAGTATGGCGGTGCCAGTGTCCCGGGGTTGGGAGGGGTAAGGATTGGGGGGGCCAACAGCGGGGCGAGGCCGGTGTTTTCGGGGGGGGGGACATCGGGG‑AACACGGTAAAAACCCTGCTGGGGGGGGAGGtg                                      >  1:3044424/1‑150 (MQ=255)
                                                                                             aGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGGGGCGGGCTCGGTGTGTCCGGGCGGGGTGAATTCGGGC‑AACCCGGGAAAACCGCCGCTGGGGGGGGAgggg                                      >  2:3828908/1‑144 (MQ=255)
                                                                                              gggAGATATGACTCAGTATGGCGATCCCAATTCCCCCGGTTCGGCAGGTGTAAGGATTTGCGCGGCCACCCGGGAGGCGCGCCCGGCGGGTCCGGGGGGGGGGACTTCGGGCAAACCCGGAAAAT‑CGCCGCTGGGCGGGAAGGAGCTGtcc                                >  1:2791644/1‑151 (MQ=255)
                                                                                               tgagaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGGCCGGGGGGGGTGGCCTCGGGCAAACCCGGAAAATCC‑CCCCTGGGGGGGGAGGGGCTGccc                                >  1:2567610/2‑150 (MQ=255)
                                                                                               ggAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGGAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGGGGGGTGACTTCGGGG‑AAACCCGGAAATCCGCCGCTTGGGGCGGAGGGGCTGccc                                >  1:973274/1‑150 (MQ=255)
                                                                                               ggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCg                               <  2:2268230/151‑1 (MQ=255)
                                                                                               ggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCg                               <  2:2160190/151‑1 (MQ=255)
                                                                                                   atatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCga                           <  2:3272879/151‑1 (MQ=255)
                                                                                                   atatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCga                           <  1:615980/151‑1 (MQ=255)
                                                                                                     atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgaga                         <  2:2108480/151‑1 (MQ=255)
                                                                                                           cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGGATTTGCGCGCCCACCGGCGTGGCGTGGTCGGTGTGTCCCGGCGGGGTGGCTTCGGGG‑AAACCCGTAAATCCGCTGCTGGGGGCGGAAGTGCTGCCCGGGGAGACg                       >  2:3635336/1‑147 (MQ=255)
                                                                                                                             gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg  <  1:3246879/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                     
TGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCG  >  NC_000913/523193‑523467

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