Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,938,255 T→A 7.2% L10Q (CTG→CAG)  rbsR → DNA‑binding transcriptional dual regulator RbsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,938,2550TA 7.2% 159.2 / 4.5 57L10Q (CTG→CAG) rbsRDNA‑binding transcriptional dual regulator RbsR
Reads supporting (aligned to +/- strand):  ref base T (22/30);  new base A (2/2);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

TTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAAC  >  NC_000913/3938122‑3938392
                                                                                                                                     |                                                                                                                                         
ttcgCCCAGGCTGCCGCTGCGATTGCTGTAACACGGAAAGGTGCACAACCTGACGTACCGTGGCGTGAAGAGATCGACGCATTGTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCAggcgggcaa                                                                                                                                  <  2:554741/140‑3 (MQ=255)
tgtGCCCAGGCTGCCGCTGCGATTGCCGTAACGCGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGAGCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCAggcgggcaa                                                                                                                                  <  2:135801/141‑3 (MQ=255)
     ccATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTa                                                                                                                             <  1:426670/143‑1 (MQ=255)
        tgcggcCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCt                                                                                                                          <  1:417886/143‑1 (MQ=255)
               gCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGt                                                                                                                   <  2:627477/143‑1 (MQ=255)
                 tGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGttt                                                                                                                 <  1:239730/143‑1 (MQ=255)
                 tGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGAAGAGGGGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGttt                                                                                                                 >  1:299512/1‑143 (MQ=255)
                     aTTGCCGTAACACGTCAGGGTGTACATCCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCa                                                                                                             <  1:499679/143‑1 (MQ=255)
                     aTTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCGACAATGAACGATGTTGCGCGCCTGGCGGGCGATTGTACATCAACAGTTTCTCa                                                                                                             >  2:635182/1‑143 (MQ=255)
                           gTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGGGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCCCGttct                                                                                                       >  1:247987/1‑141 (MQ=255)
                            tAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTCATc                                                                                                      >  1:27955/1‑143 (MQ=255)
                              acacGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCaa                                                                                                    >  2:223617/1‑143 (MQ=255)
                              acacGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCaa                                                                                                    >  1:713338/1‑143 (MQ=255)
                              acacGTAAAGGCGCAAAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCaa                                                                                                    >  2:680093/1‑143 (MQ=255)
                                 cGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATaa                                                                                                 >  2:799495/1‑143 (MQ=255)
                                  gTACAGGCGCACGACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACTCTTGGCTACAATGAAAGATGTTGCCCGCTTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATaaa                                                                                                <  2:28400/143‑1 (MQ=255)
                                  gTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATaaa                                                                                                >  1:460479/1‑143 (MQ=255)
                                    aaaGGCGCACAACTTGCCGGCCCGTGGCGTGAAGAGATCGAGGCATTTTTCGACTGGCAGAGGCGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGa                                                                                              <  2:276185/143‑1 (MQ=255)
                                     aaGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGAGGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                                                             >  1:260867/1‑143 (MQ=255)
                                      aGGCGCACAACCTTCCGCACCGTGGGGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTTGCTACAATGAAAGATGTTGCTCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATc                                                                                            <  1:589068/143‑1 (MQ=255)
                                          gcACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGCCAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCtt                                                                                        <  1:246351/143‑1 (MQ=255)
                                             caACCTTCCGGACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGt                                                                                     <  1:1050/143‑1 (MQ=255)
                                                ccTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAg                                                                                  <  1:304804/143‑1 (MQ=255)
                                                 cTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGGTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCggccgt                                                                                 >  1:504685/1‑138 (MQ=255)
                                                  gtCCGTACCGGGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTg                                                                                <  1:718478/142‑1 (MQ=255)
                                                    ccGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGac                                                                              >  1:356050/1‑142 (MQ=255)
                                                    ccGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGaa                                                                              >  2:454076/1‑143 (MQ=255)
                                                         ccGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGAt                                                                         <  2:636198/143‑1 (MQ=255)
                                                           gTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCACGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGaaccaatca                                                                       >  2:807567/1‑136 (MQ=255)
                                                             ggCGTGAAGAGATCGACGCAGTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTAcc                                                                     <  1:94687/143‑1 (MQ=255)
                                                              gCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCAGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCg                                                                    <  1:485455/143‑1 (MQ=255)
                                                                 tGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCACTAAGGATCGCTTCGTCAGTGAAGCGATTACCGCCa                                                                 >  2:24326/1‑143 (MQ=255)
                                                                   aaGAGTTCGATGCATTTGTGGACGGGCATAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGTTTCGTCAGTGAAGCGATTACCGCCaaa                                                               <  2:288941/143‑1 (MQ=255)
                                                                     tagaTCGACCTAGTTTTAGCCAGGCAGGGGTGTCGCTTGGCTACAATGCAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGt                                                             <  2:459954/142‑1 (MQ=255)
                                                                     gagaTCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGt                                                             <  1:798916/143‑1 (MQ=255)
                                                                      agaTCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGtt                                                            >  1:84816/1‑143 (MQ=255)
                                                                         tCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTgca                                                         >  1:177475/1‑141 (MQ=255)
                                                                         tCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGaa                                                         <  1:509778/143‑1 (MQ=255)
                                                                              gCATTTTTAGACAGGCAGAGGTGACGCTTGGCTCCAATGAAAGACGTCGCCCGCCTGGCGGGCGTTTCTACCGCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgcggc                                                    <  1:156812/143‑1 (MQ=255)
                                                                              cgaTTTTTAGGCAGGCAGAGATGAAGCTTGGCTACAATGAAAGATGGTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGGCAGTGAAGCGATTACCGCCAAAGTTGAAgcggc                                                    <  2:418876/141‑1 (MQ=255)
                                                                                  ttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTGTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGAtt                                                >  2:749819/1‑143 (MQ=255)
                                                                                       gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGGTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAgc                                           >  1:504667/1‑142 (MQ=255)
                                                                                        aCAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTCAAGaa                                          >  1:127077/1‑143 (MQ=255)
                                                                                                  gTGACGCTTGTCTACTATGAATGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGCTATCAATGAAGATGGTTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTAcg                                <  2:322609/143‑1 (MQ=255)
                                                                                                  gTGACGCTTCGCTACAATGTAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTAcg                                <  1:794737/143‑1 (MQ=255)
                                                                                                   tGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTAcgc                               <  2:239650/143‑1 (MQ=255)
                                                                                                   tGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACCCAATTAcgc                               <  1:410053/143‑1 (MQ=255)
                                                                                                         ttGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAAGAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATc                         <  1:412532/143‑1 (MQ=255)
                                                                                                          tGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTCAAGAACTCACTTACGCGCCATcc                        >  2:734564/1‑142 (MQ=255)
                                                                                                              tACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCACCAGctc                    >  2:270029/1‑143 (MQ=255)
                                                                                                                     aaaGCTGCTGGCCGCCCGGCGGGCGTGTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTCGAAGCAAACATTAAAGAACTCAATTACGCGCCATCAGCTCTGGcgcg             <  2:574935/143‑1 (MQ=255)
                                                                                                                        gctGTTGCCCGCCTGGCGGGCGTTTCTACGTTAACAGGTTGTCACGTTTTCAATGACGATCGCTTCGTCAGTGTAGCGATTCCCGGCAAAGTTGACGCGGCGTTTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAg          <  2:283773/141‑1 (MQ=255)
                                                                                                                        gATGTCGCCCGGCTGGTGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTAAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAg          <  2:615424/143‑1 (MQ=255)
                                                                                                                            acgCCCGCCAGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGCAGCGGCGATTAACGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTc      >  1:604755/3‑143 (MQ=255)
                                                                                                                              gCCCGCCAGGCGGGCATTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGAGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAAAACTCAATTAAGCGCCATCAGCTCTGGCGCGTACCCCCac    >  2:158853/1‑142 (MQ=255)
                                                                                                                               cccGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCaaa   <  1:779214/143‑1 (MQ=255)
                                                                                                                                ccGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAAc  <  1:454076/143‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                         
TTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAAC  >  NC_000913/3938122‑3938392

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