Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3350TG53.3% 23.9 / 38.4 30G25G (GGT→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base G (16/0);  total (24/6)
Fisher's exact test for biased strand distribution p-value = 5.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-04

GCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTC  >  NC_000913/523185‑523480
                                                                                                                                                      |                                                                                                                                                   
gCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGgt                                                                                                                                                     <  2:460157/151‑1 (MQ=255)
         gtgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCCGGGGGAAGAAtt                                                                                                                                            >  2:185507/1‑151 (MQ=255)
         gtgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAAt                                                                                                                                             >  1:49017/1‑150 (MQ=255)
                attttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGGAGGGGTAAGAATTTgcg                                                                                                                                        >  2:574148/1‑148 (MQ=255)
                       ttttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTGCCATTGTCCAAGGGTCGGCAGGGGTAAGAATTGGGGCCGCCACCgg                                                                                                                              >  1:388168/1‑151 (MQ=255)
                           ttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCggcgtg                                                                                                                          >  1:668707/1‑151 (MQ=255)
                            tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtg                                                                                                                          <  1:369526/150‑1 (MQ=255)
                                ttgCTAACACAGGGCTATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGAAATGACTCAGTATGGCGGTCCCACTGTCCAGGGTTTGGGAGGGGTAAGAATTTGCGGGCCCACCggcgcggtg                                                                                                                       >  2:67090/1‑147 (MQ=255)
                                ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCcgcgtggcgt                                                                                                                      >  2:614312/1‑150 (MQ=255)
                                         caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtgtg                                                                                                            <  2:433082/151‑1 (MQ=255)
                                               tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCCCGGTGTGTGcggg                                                                                                      >  2:507296/1‑151 (MQ=255)
                                                     cAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAt                                                                                                <  1:352183/151‑1 (MQ=255)
                                                          aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGGGTCCGGGCGGGGTG‑ACtt                                                                                          >  1:422355/1‑151 (MQ=255)
                                                             aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTc                                                                                         <  2:477264/149‑1 (MQ=255)
                                                              aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAAGGTTCGGCAGGGGTAAGAAATGGCGCGCCCCCCGGCGTGGCGGGCGCCGTGGGTCCCGGGGGGGTG‑GCCTCggg                                                                                      >  2:550457/1‑151 (MQ=255)
                                                              aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTCCGGGCGGGATG‑ACTTCggg                                                                                      >  1:646947/1‑151 (MQ=255)
                                                                      cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGGG‑ACTTCGGGCAAccc                                                                                >  2:172532/1‑149 (MQ=255)
                                                                       tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGGGG‑ACCTCGGGCAACCCgg                                                                              >  2:371951/1‑150 (MQ=255)
                                                                         tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGTCCTTTGCCTAGGGGTAAGGAGGGGTTAGAAATGGTGGGCCCCACGGCGTGGGGTGGTCGGGGTGGACCGGCGGGATGAACTTCGGCAAAccg                                                                                >  1:512674/1‑146 (MQ=255)
                                                                          atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCACTGTCCCGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCACCCGCGTGGCGGGCTCGGTGTGTCCGGGCGGGGTG‑ACTTCGGGGAACCCGGGaa                                                                           >  2:419416/1‑150 (MQ=255)
                                                                                          aGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGGCGGGGTAAGAATTGGCGCGCCCCCCCGCGTGGGGTGGCCGGTGTGGTCGGGCGGGGTG‑ACTTCGGGGAACCCGGGGAAACCC‑CTGCTggggg                                                           >  2:228762/1‑149 (MQ=255)
                                                                                              gcgcGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGGGGGTCCGGGGGGGGGG‑GCTTCCGGCAACCCCGGAAAACCGCCTGCTggggcggagg                                                      >  2:485093/1‑145 (MQ=255)
                                                                                                   tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGTGCTCGGTGTGTCCGGGCGGGGTG‑ACTTCCGGCAAACCCGGAAATCCG‑CTGCTGGGGGCGGAGGggct                                                 >  1:40857/1‑150 (MQ=255)
                                                                                                     aGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGGGTGTCCGGGGGGGGTG‑ACTTCGGGCAACCCGGGAAAACCG‑CTGCTGGGGGCGGAGGGGCTGcc                                              >  2:46740/1‑151 (MQ=255)
                                                                                                          gaTATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCGGGGTGTCCGGGGGGGGGG‑ACTTCGGGCAAACCGGGAAAACCG‑CCGCTGGGGGGGAAGGGGCTGCCCcgcg                                         >  1:35984/1‑151 (MQ=255)
                                                                                                                            ggTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgc                       <  1:203452/151‑1 (MQ=255)
                                                                                                                               cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                     <  2:388158/150‑1 (MQ=255)
                                                                                                                                 cATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggc                  >  2:375333/1‑151 (MQ=255)
                                                                                                                                     gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGc              >  1:84293/1‑151 (MQ=255)
                                                                                                                                              tCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTc      >  2:652942/1‑150 (MQ=255)
                                                                                                                                                 gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtc  >  1:170445/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATG‑ACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTC  >  NC_000913/523185‑523480

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