Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3350TG41.9% 68.0 / 12.0 44G25G (GGT→GGGrhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (7/18);  new base G (18/0);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 8.66e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523188‑523483
                                                                                                                                                   |                                                                                                                                                      
ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGgtg                                                                                                                                                       <  2:509995/149‑1 (MQ=255)
     tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGaa                                                                                                                                                 >  1:708829/1‑150 (MQ=255)
            cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGgcgc                                                                                                                                          <  2:736013/150‑1 (MQ=255)
                      ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggc                                                                                                                                <  2:545807/150‑1 (MQ=255)
                           tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCCGGGGTCGGCCGGGGGAAGAATTGGCGCGCCCACCcgcg                                                                                                                               >  1:518060/1‑146 (MQ=255)
                           tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtgg                                                                                                                            <  1:394742/149‑1 (MQ=255)
                            gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCCGGGTTCGGCAGGGGTAAGAAATTGCGGCCCCACCcgcgtgggg                                                                                                                          >  2:94206/1‑148 (MQ=255)
                            gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGTAAGAAATTGCGCGCCCACCcgcgggggg                                                                                                                          >  2:422493/1‑148 (MQ=255)
                                  aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTc                                                                                                                     <  1:423537/149‑1 (MQ=255)
                                           atatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGGAGGGGTAAGAAATTGCGCGCCCCCCCGCGTGGGGGGGTCGgggtggtccg                                                                                                          >  2:679754/1‑146 (MQ=255)
                                             atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcggg                                                                                                         <  1:370788/150‑1 (MQ=255)
                                                aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcggg                                                                                                     <  1:667101/151‑1 (MQ=255)
                                                     tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGGAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTcggggggg                                                                                                     >  1:641005/1‑146 (MQ=255)
                                                             tgggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGGGGAAGAATTTGCGCGCCCCCCCGCGTGGGGGGCTCGGGGTGTCCGGGCGGGGTGACtttcgg                                                                                         >  2:311291/2‑145 (MQ=255)
                                                                 cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCCGGGTGTCCCGGCGGGGGGACTTCGGGCAAccc                                                                                   >  1:327074/1‑151 (MQ=255)
                                                                 cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGTGCTCGGTGGGGCCGGGGGGGGTGGCTTTGGGCAAacc                                                                                   >  2:473111/1‑151 (MQ=255)
                                                                   cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCAGGGGTCGGCAGGGGTAAGAAATGGGGCGCCCCCCGGCGTGGGGTGCTCGGTGTGTCCGGGGGGGGTGACGTCGGGCAAccc                                                                                   >  2:133169/1‑149 (MQ=255)
                                                                     ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCACTGTCCAGGGGTCGGGAGGGGGTAGAATTTGCGCGCCCACCGGCGGGGCGGGGTCGGGGGGTCCCGGCGGGGGGACCTCCGGCAAACCgg                                                                                 >  2:465600/1‑149 (MQ=255)
                                                                      tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTTGGCCGGGGTAAGAAATGGGGCGCCCACCCGCGTGGGGGGCTCGGGGTGGCCGGGGGGGGTGGATTCGGGCAAacc                                                                                   >  2:635907/1‑146 (MQ=255)
                                                                      tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTTAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCTGGGTGTCCCGGCGGGGTGGATTCGGGCAACCCggga                                                                               >  2:197234/1‑148 (MQ=255)
                                                                         atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCCGGGGTAAGAATTTGGGCGCCCACCCGCGTGGCGGGCTCGGGGGGTCCCGGGGGGGTGGCTTCGGGGCACCCCGGaaa                                                                             >  2:478967/1‑149 (MQ=255)
                                                                         atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAAt                                                                            <  1:413637/150‑1 (MQ=255)
                                                                         atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAAt                                                                            <  2:42167/150‑1 (MQ=255)
                                                                                gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctgct                                                                   <  1:559488/151‑1 (MQ=255)
                                                                                 aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTTCGGGCGGGGTGGCCTTGGGCCACCCCGGAAAATccccgc                                                                      >  2:219937/1‑145 (MQ=255)
                                                                                 aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGCGGCGTGCTCGGGGTGGCCGGGGGGGGTGAACTCGGGGAAACCGGGAAAACCgc                                                                       >  2:509160/1‑146 (MQ=255)
                                                                                    aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgggg                                                               <  2:224673/151‑1 (MQ=255)
                                                                                      cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAAGGGTCGGCGGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGGTCGGTGGGGCCGGGGGGGGTGGATTCTGGCAACCCGGGAAATTCGCTCCTggggg                                                              >  1:705472/1‑150 (MQ=255)
                                                                                         cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGaa                                                          <  1:569539/151‑1 (MQ=255)
                                                                                            cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGGGTGGCTGGGGTGTCCCGGCGGGGGGACTTCCGGGAACCCGGGGAATCCCCCGCTGGGGGCgccggg                                                       >  1:613763/1‑146 (MQ=255)
                                                                                              cgTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCCCCCGCGTGGCGGGCCCGGTGTGTCCGGGGGGGGTGGCCTCGGGCAACCCCGGAAAACCGCCGCTGGGGGCGGAg                                                         >  2:61051/1‑147 (MQ=255)
                                                                                              cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGTGTGTTCGGGGGGGGTGGCTTCGGGCAACCCCGGAAATCCCCTGCCGGGGGCGGAgggg                                                      >  2:303561/1‑148 (MQ=255)
                                                                                                    ggggaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGTGCCGACCGGCGTGGCGTGCTCGGTGTGTCCGGACGGGATGCCTTCGGGCAGCCCGGTAAATCCGCTGCTGGGGAAGAAGGTTCTACCCg                                               <  2:192076/148‑1 (MQ=255)
                                                                                                    ggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTCGGGCGGGGTGACTTTGGGCAACCCCGGAAAACCGCTGCTGGGGGGGAAGGGGCTGccc                                                >  1:589909/1‑150 (MQ=255)
                                                                                                         tatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGCGTAAGAATTGGCGCCCCCACCGGCGTGGCGTGCTCGGTGTGGCCGGGCGGGGTGACTTCGGGCAACCCCGCAAAACCGCTGGTGGGGGGGAAGGGGCCGCCCGGCg                                            >  2:384122/1‑149 (MQ=255)
                                                                                                          atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgag                                          <  2:312354/150‑1 (MQ=255)
                                                                                                                cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGa                                     <  1:568731/149‑1 (MQ=255)
                                                                                                                cAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGGAAGAATTGGGGCGGCCACCCGCGTGGCGTGCGCGGTGTGTCCGGGGGGGGGTACTTTGGGGAACCCCGGAAATCCGCCTCCTGGGGGGGAGGTGCTGTCCCGCGGGGCgg                                      >  2:285285/1‑148 (MQ=255)
                                                                                                                  gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                                 <  1:50633/151‑1 (MQ=255)
                                                                                                                      ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcgc                             >  2:557184/1‑151 (MQ=255)
                                                                                                                        cGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTg                           >  1:308134/1‑151 (MQ=255)
                                                                                                                           tCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                        >  1:156446/1‑151 (MQ=255)
                                                                                                                            cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg                       <  1:362501/151‑1 (MQ=255)
                                                                                                                                 tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                   <  2:768900/150‑1 (MQ=255)
                                                                                                                                    ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTg               <  1:59193/151‑1 (MQ=255)
                                                                                                                                                  ctgtAAGAATTCGCGCGCCCACCGGCGTGTCGTGCTCGGGGAGTCCGAGCGGGATGACTTCCAGAAACCCCGCAAATCCGCTGCTGGGGGTGCAGGTGATGCCCGGCGAGACGGACCTTGCCCTGCCCGGCCCGCCGCCGATCCTtctct  >  2:656688/2‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523188‑523483

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: