Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 524,158 C→T L300L (CTG→TTG)  rhsD → Rhs protein with DUF4329 family 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_000913524,1580CT100.0% 212.2 / NA 67L300L (CTG→TTG) rhsDRhs protein with DUF4329 family putative toxin domain, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/32);  total (35/32)

GCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGCTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGT  >  NC_000913/524022‑524284
                                                                                                                                        |                                                                                                                              
gCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTg                                                                                                                              <  2:410812/139‑1 (MQ=255)
     cGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                         >  1:1312155/1‑139 (MQ=255)
         agagGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGAc                                                                                                                     >  2:941863/1‑139 (MQ=255)
         agagGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGAc                                                                                                                     <  2:1312155/139‑1 (MQ=255)
             gCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                 >  1:711957/1‑139 (MQ=255)
             gCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                 <  1:941863/139‑1 (MQ=255)
             gCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                 >  1:472187/1‑139 (MQ=255)
             gCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                 <  1:1405448/139‑1 (MQ=255)
                                      gtataagagacagCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACAACCCGGCATACCCGGAGAGCCTGCCCGCTgc                                                                                        >  1:1534055/13‑139 (MQ=255)
                                         aCAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCg                                                                                            >  1:1537053/1‑132 (MQ=255)
                                         aCAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCg                                                                                            <  2:1537053/132‑1 (MQ=255)
                                         aCAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGcc                                                                                     >  2:121599/1‑139 (MQ=255)
                                             ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTg                                                                                 >  2:606782/1‑139 (MQ=255)
                                             ttCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTg                                                                                 >  2:1657407/1‑139 (MQ=255)
                                                ccGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                              <  2:853163/139‑1 (MQ=255)
                                                ccGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                              <  2:472187/139‑1 (MQ=255)
                                                   cccTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCgg                                                                           <  1:1657407/139‑1 (MQ=255)
                                                      tctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTac                                                                        <  2:711957/139‑1 (MQ=255)
                                                       ctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                              >  1:458373/1‑132 (MQ=255)
                                                       ctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTg                                                                              <  2:458373/132‑1 (MQ=255)
                                                       ctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTaca                                                                       >  1:707317/1‑139 (MQ=255)
                                                       ctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTaca                                                                       >  2:1119189/1‑139 (MQ=255)
                                                       ctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTCCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTgg                                                                                >  1:1235494/1‑130 (MQ=255)
                                                       ctctcAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTCCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTgg                                                                                <  2:1235494/130‑1 (MQ=255)
                                                                  agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                 >  2:250490/1‑86 (MQ=255)
                                                                  agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCgg                                                                                                                 <  1:250490/86‑1 (MQ=255)
                                                                  agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGaa                                                            >  2:1266517/1‑139 (MQ=255)
                                                                  agcGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACAACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGaa                                                            <  2:1534055/139‑1 (MQ=255)
                                                                          ccGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGa                                                    <  1:1119189/139‑1 (MQ=255)
                                                                           cGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGaa                                                   <  2:1094595/139‑1 (MQ=255)
                                                                              cacacTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAActg                                                <  2:707317/139‑1 (MQ=255)
                                                                                cacTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                         >  1:1092357/1‑64 (MQ=255)
                                                                                cacTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                         <  2:1092357/64‑1 (MQ=255)
                                                                                 acTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAActgctg                                             >  2:607758/1‑139 (MQ=255)
                                                                                 acTGCCCGG‑ACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGAc                                                                                                                     >  1:1376474/1‑66 (MQ=255)
                                                                                 acTGCCCGG‑ACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGAc                                                                                                                     <  2:1376474/66‑1 (MQ=255)
                                                                                   tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCt                                                                                                 <  2:755646/85‑1 (MQ=255)
                                                                                   tGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCt                                                                                                 >  1:755646/1‑85 (MQ=255)
                                                                                     cccGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCgg                                         >  1:762659/1‑139 (MQ=255)
                                                                                        ggTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTa                                                                         >  2:1681634/1‑104 (MQ=255)
                                                                                        ggTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTa                                                                         <  1:1681634/104‑1 (MQ=255)
                                                                                          tACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTg                                                                                 <  1:1617488/94‑1 (MQ=255)
                                                                                          tACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTg                                                                                 >  2:1617488/1‑94 (MQ=255)
                                                                                          tACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTatat                                    <  1:606782/139‑1 (MQ=255)
                                                                                           aCCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATg                                   >  1:1141779/1‑139 (MQ=255)
                                                                                               aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                         >  2:945715/1‑49 (MQ=255)
                                                                                               aaTACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCa                                                                                                                         <  1:945715/49‑1 (MQ=255)
                                                                                                    ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCa                                                                                    >  2:1279787/1‑81 (MQ=255)
                                                                                                    ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCa                                                                                    <  1:1279787/81‑1 (MQ=255)
                                                                                                    ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTa                                                                         <  1:350805/92‑1 (MQ=255)
                                                                                                    ggCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTa                                                                         >  2:350805/1‑92 (MQ=255)
                                                                                                      ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATg                                                                                                                           >  2:1356135/1‑40 (MQ=255)
                                                                                                      ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATg                                                                                                                           <  1:1356135/40‑1 (MQ=255)
                                                                                                      ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAAt                        >  2:792240/1‑139 (MQ=255)
                                                                                                      ccccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAAt                        >  2:1255334/1‑139 (MQ=255)
                                                                                                       cccGACAGGGGTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATa                       >  2:541897/1‑139 (MQ=255)
                                                                                                              ggggTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGa                                                             <  2:1016956/94‑1 (MQ=255)
                                                                                                              ggggTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGa                                                             >  1:1016956/1‑94 (MQ=255)
                                                                                                               gggTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAActg                                                >  1:158505/1‑106 (MQ=255)
                                                                                                               gggTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAActg                                                <  2:158505/106‑1 (MQ=255)
                                                                                                                ggTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCAc                                                                                   <  2:680961/70‑1 (MQ=255)
                                                                                                                ggTATCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCAc                                                                                   >  1:680961/1‑70 (MQ=255)
                                                                                                                   aTCCGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTgcgcg           <  1:1255334/139‑1 (MQ=255)
                                                                                                                     ccGCCTTTCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCt         <  1:792240/139‑1 (MQ=255)
                                                                                                                         cTTTCGGCGGTGTGGTTGATGTACGACCCGGCATACCCGGAGAGCCTGCCCGCTGTGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCt         <  1:1112553/135‑1 (MQ=255)
                                                                                                                         cTTTCGGCGGTGTGGTTGATGTACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACTCGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCt         >  2:1112553/1‑135 (MQ=255)
                                                                                                                            tCGGCGGTGTGGTTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTTACGt  >  1:1701570/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                              
GCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCGGTGTGGCTGATGCACGACCCGGCATACCCGGAGAGCCTGCCCGCTGCGCCACTGGTGCGGTACACGTATACGGAAGCCGGTGAACTGCTGGCGGTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGT  >  NC_000913/524022‑524284

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