Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,622,796 T→G L1202R (CTG→CGG)  rhsB → 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_0009133,622,7960TG94.1% 168.5 / 3.7 68L1202R (CTG→CGG) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base G (39/25);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATA  >  NC_000913/3622699‑3622892
                                                                                                 |                                                                                                
tatGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGt                                                                                                <  2:1562938/100‑1 (MQ=38)
   gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                             <  2:1168712/100‑1 (MQ=18)
    atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa                                                                                            <  1:1892552/100‑1 (MQ=18)
     tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTAc                                                                                           >  2:93300/1‑100 (MQ=18)
     tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTAc                                                                                           >  2:510661/1‑100 (MQ=17)
        aTGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAAc                                                                                        <  2:1529877/100‑1 (MQ=18)
        aTGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAAc                                                                                        <  2:1530020/100‑1 (MQ=18)
            ggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGcc                                                                                    <  1:391021/100‑1 (MQ=18)
               aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                 <  1:1539421/100‑1 (MQ=18)
               aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                 >  1:777584/1‑100 (MQ=18)
                      tgaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatta                                                                          >  2:1646390/1‑100 (MQ=18)
                       gaatgaaGGGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat                                                                         >  1:1728743/1‑100 (MQ=255)
                       gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat                                                                         >  2:485370/1‑100 (MQ=18)
                         atgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGa                                                                       >  2:1896193/1‑100 (MQ=18)
                           gaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAcc                                                                     >  2:1431741/1‑100 (MQ=18)
                            aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                    <  2:1386319/100‑1 (MQ=18)
                             agagAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATAATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                    >  2:720670/1‑99 (MQ=255)
                              gagaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGc                                                                  >  1:607880/1‑100 (MQ=18)
                                gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTg                                                                >  2:1108576/1‑100 (MQ=18)
                                gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTg                                                                >  1:571321/1‑100 (MQ=18)
                                   cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                             >  1:1094432/1‑100 (MQ=33)
                                   cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                             >  2:1924285/1‑100 (MQ=18)
                                      gCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgggg                                                          <  1:215381/100‑1 (MQ=18)
                                      gCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgggg                                                          >  1:225839/1‑100 (MQ=18)
                                       cATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGc                                                         <  2:1625753/100‑1 (MQ=18)
                                        aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                        <  2:1485348/100‑1 (MQ=18)
                                            gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                    >  2:438193/1‑100 (MQ=18)
                                            gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                    >  1:1406775/1‑100 (MQ=17)
                                              tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatat                                                  >  2:692630/1‑100 (MQ=18)
                                               gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatat                                                  >  2:1377520/1‑99 (MQ=18)
                                               gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                 <  1:454037/100‑1 (MQ=17)
                                               gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                 >  2:1095015/1‑100 (MQ=18)
                                               gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                 <  2:1032655/100‑1 (MQ=17)
                                                   cagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTc                                             >  2:387884/1‑100 (MQ=18)
                                                   cagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTc                                             <  1:1414589/100‑1 (MQ=18)
                                                    agcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCa                                            <  2:1437116/100‑1 (MQ=18)
                                                     gcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAg                                           >  1:1557454/1‑100 (MQ=18)
                                                        tATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                        >  2:774616/1‑100 (MQ=18)
                                                        tATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                        >  2:462831/1‑100 (MQ=18)
                                                        tATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                        >  2:1417301/1‑100 (MQ=18)
                                                           ccGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCg                                     >  2:859144/1‑100 (MQ=17)
                                                           ccGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCg                                     <  1:692629/100‑1 (MQ=18)
                                                             ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                   <  1:297526/100‑1 (MQ=14)
                                                                 gCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGa                               >  1:979345/1‑100 (MQ=11)
                                                                 gCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGa                               >  1:249042/1‑100 (MQ=11)
                                                                  ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              <  2:1358408/100‑1 (MQ=14)
                                                                  ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              >  2:125211/1‑100 (MQ=14)
                                                                  ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              <  2:607882/100‑1 (MQ=14)
                                                                   cGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACt                             >  2:1654138/1‑100 (MQ=32)
                                                                       gcagcaGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAg                         <  1:1187423/100‑1 (MQ=35)
                                                                       gcagcaGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGaan                         <  2:1761921/100‑2 (MQ=18)
                                                                            agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                    <  2:683920/100‑1 (MQ=255)
                                                                            agTATGATGAGGAGTCCGGCCGGTATAACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                    >  2:1186780/1‑100 (MQ=255)
                                                                             gTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGa                   <  1:1924272/100‑1 (MQ=255)
                                                                                tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATga                >  1:1070629/1‑99 (MQ=255)
                                                                                 gatgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACCGAAGGGGGGATGGa               >  1:549553/1‑100 (MQ=25)
                                                                                   tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc             >  2:17144/1‑100 (MQ=255)
                                                                                    gaggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc            >  1:682371/1‑100 (MQ=255)
                                                                                    gaggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc            >  2:1274167/1‑100 (MQ=255)
                                                                                      ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg          >  1:1302545/1‑100 (MQ=255)
                                                                                      ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGACATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg          >  2:1881173/1‑100 (MQ=255)
                                                                                       gagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGt         >  1:351367/1‑100 (MQ=255)
                                                                                       gagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGt         >  2:851874/1‑100 (MQ=255)
                                                                                        agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa        >  2:1839373/1‑100 (MQ=255)
                                                                                         gTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt       >  2:1756613/1‑100 (MQ=255)
                                                                                           ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGTATCCGCTTGGACTGCAGGGGGGATGGAACCTGTATgg     >  1:1308045/1‑100 (MQ=17)
                                                                                           ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg     >  1:1095392/1‑100 (MQ=255)
                                                                                           ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg     >  2:1042933/1‑100 (MQ=255)
                                                                                             ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat   <  1:112548/100‑1 (MQ=255)
                                                                                             ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat   <  1:972128/100‑1 (MQ=255)
                                                                                              gCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGata  <  2:443622/100‑1 (MQ=255)
                                                                                                 |                                                                                                
TATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATA  >  NC_000913/3622699‑3622892

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