Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,376 G→T 7.6% intergenic (+36/‑207) ybfA → / → rhsC DUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913729,3760GT 7.6% 172.2 / 3.0 52intergenic (+36/‑207)ybfA/rhsCDUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (18/30);  new base T (2/2);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATA  >  NC_000913/729230‑729524
                                                                                                                                                  |                                                                                                                                                    
tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGaa                                                                                                                                                    <  3:620156/149‑1 (MQ=255)
  gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                                  <  3:272685/149‑1 (MQ=255)
    aaGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGc                                                                                                                                                >  3:394595/1‑149 (MQ=255)
       aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTgg                                                                                                                                             <  4:394424/149‑1 (MQ=255)
        ccGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAATCCGCCAGCGGAACCGCTGGc                                                                                                                                            <  4:1293440/149‑1 (MQ=255)
          gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCgg                                                                                                                                          <  3:665681/149‑1 (MQ=255)
          gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCgg                                                                                                                                          >  2:217781/1‑149 (MQ=255)
                gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                                    >  4:739980/1‑149 (MQ=255)
                          cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa                                                                                                                          >  2:716384/1‑149 (MQ=255)
                          cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa                                                                                                                          <  1:716384/149‑1 (MQ=255)
                             ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCACCTGGCGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAt                                                                                                                       <  3:1376479/149‑1 (MQ=255)
                             ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAt                                                                                                                       >  3:291477/1‑149 (MQ=255)
                                  aTCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                                                  >  1:719412/1‑149 (MQ=255)
                                  aTCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                                                  <  4:394595/149‑1 (MQ=255)
                                    cgcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                                                <  1:217781/149‑1 (MQ=255)
                                     gcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCagag                                                                                                               <  3:695876/149‑1 (MQ=255)
                                               gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                        <  3:430770/114‑1 (MQ=255)
                                               gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                        >  4:430770/1‑114 (MQ=255)
                                                      aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGt                                                                                                                                         >  3:1062768/1‑106 (MQ=39)
                                                      aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGt                                                                                                                                         <  4:1062768/106‑1 (MQ=39)
                                                      aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                                              <  1:1082624/149‑1 (MQ=255)
                                                        cGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                            <  3:739980/149‑1 (MQ=255)
                                                           tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAc                                                                                         <  1:724550/149‑1 (MQ=255)
                                                               ccGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                            <  2:32621/142‑1 (MQ=255)
                                                               ccGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                            >  1:32621/1‑142 (MQ=255)
                                                                      ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                >  2:1464903/1‑147 (MQ=255)
                                                                      ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                <  1:1464903/147‑1 (MQ=255)
                                                                       gatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                             >  1:971584/1‑149 (MQ=255)
                                                                        atggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACata                                                                            <  2:198288/149‑1 (MQ=255)
                                                                         tggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatat                                                                           <  4:1150890/149‑1 (MQ=255)
                                                                           gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                      >  3:40075/1‑120 (MQ=255)
                                                                           gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                      <  4:40075/120‑1 (MQ=255)
                                                                            atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATAGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg                                                                                                                      >  1:967573/1‑103 (MQ=255)
                                                                            atCAGGCCGAAAAGGCAACTGGTGATTTTTAGTGATTGCTGAGTTTTGAATACATTGAAAACCTCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg                                                                                                                      <  2:967573/103‑1 (MQ=255)
                                                                             tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTCGTATCAGGAAGACGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTc                                                                       <  4:1070864/149‑1 (MQ=255)
                                                                               aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                >  2:909872/1‑138 (MQ=255)
                                                                               aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                <  1:909872/138‑1 (MQ=255)
                                                                               aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                                     <  2:971584/149‑1 (MQ=255)
                                                                                       aaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatata                                                             <  3:239195/149‑1 (MQ=255)
                                                                                        aGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatat                                                            <  4:291477/149‑1 (MQ=255)
                                                                                          gCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCaa                                                                                      >  4:993631/1‑121 (MQ=255)
                                                                                          gCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCaa                                                                                      <  3:993631/121‑1 (MQ=255)
                                                                                                               ttgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCt                                     >  2:1236252/1‑149 (MQ=255)
                                                                                                                   tGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc                                 >  1:79624/1‑149 (MQ=255)
                                                                                                                      gTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATc                                                                                    <  3:754366/95‑1 (MQ=255)
                                                                                                                      gTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATc                                                                                    >  4:754366/1‑95 (MQ=255)
                                                                                                                           gAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGat                         >  1:413798/1‑149 (MQ=255)
                                                                                                                             aTACATTGAAACCCGCCAGCGGTAACGCTGGCGGTTTTTTGTATCAGGAACATGTTCTGTGGAATATTTAGATATCGTTATACCTATCGAAACTTATATTCATGCAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGatat                       <  4:308323/149‑1 (MQ=255)
                                                                                                                                   aaaaaaCCGCCAGCGTTTCCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGt                                        >  2:573791/3‑126 (MQ=39)
                                                                                                                                   aaaaaaCCGCCAGCGTTTCCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGt                                        <  1:573791/124‑1 (MQ=39)
                                                                                                                                        aCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAgg            <  1:1236252/149‑1 (MQ=255)
                                                                                                                                             cAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTc       <  2:413798/149‑1 (MQ=255)
                                                                                                                                                  gAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTtata  >  3:233957/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATA  >  NC_000913/729230‑729524

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: