Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,665,195 (T)7→8 intergenic (+75/‑120) dmsD → / → clcB twin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,665,1881.T90.3% 201.4 / 15.2 62intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base T (30/26);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281
                                                                                                    |                                                                                               
ttCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                                                                  <  2:948319/100‑1 (MQ=255)
    tGTCGCGGTTAAACCGCTGGTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGGTACAGTGTCACGtttt                                                                                            >  2:809205/1‑100 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtttttttta                                                                                      <  2:1640838/100‑3 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtttttttta                                                                                      <  2:1668122/100‑3 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttatc                                                                                    <  1:587211/100‑5 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttatc                                                                                    >  1:737308/1‑96 (MQ=255)
            ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttatct                                                                                   >  2:1361023/1‑95 (MQ=255)
            ttAAACCGCTGTTTCGATAGGATCACCGTAATATCGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttatct                                                                                   <  2:1080368/100‑6 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCtt                                                                                 <  2:733477/100‑1 (MQ=255)
                aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTa                                                                                <  2:1619164/100‑1 (MQ=255)
                aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTaa                                                                               >  1:1494555/1‑100 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCGGGATGGTGACGATTACGCGAAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaa                                                                               <  2:64567/100‑1 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTaaa                                                                              <  2:1494546/100‑1 (MQ=255)
                      gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAg                                                                             >  1:827997/1‑97 (MQ=255)
                       tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacg                                                                         <  1:520440/100‑1 (MQ=255)
                       tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAgcacgc                                                                        <  1:1841843/100‑1 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgc                                                                        <  2:1768148/100‑1 (MQ=255)
                         tCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgca                                                                       >  2:1846665/1‑100 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGCCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                      <  2:1794126/100‑1 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACt                                                                     <  1:1120005/100‑1 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACt                                                                     >  1:31195/1‑100 (MQ=255)
                            aTAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTg                                                                    >  1:2040248/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:1494478/1‑100 (MQ=255)
                                gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCttt                                                                >  2:1520767/1‑100 (MQ=255)
                                   cACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                             >  1:621724/1‑100 (MQ=255)
                                    aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            <  2:1842932/100‑1 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  1:1710695/1‑100 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  2:136482/1‑100 (MQ=255)
                                      cGTAATATGGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                          <  1:1403948/100‑1 (MQ=255)
                                       gTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCt                                                         <  2:1965979/100‑1 (MQ=255)
                                       gTAATATTGCCGGATGGTGACGATTACGCGCAATCCAGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        >  1:1441400/1‑100 (MQ=255)
                                        tAATATTGCCGGTTGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        <  2:529694/100‑1 (MQ=255)
                                        tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        <  1:1853352/100‑1 (MQ=255)
                                         aaTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                       >  1:840407/1‑100 (MQ=255)
                                          atatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc                                                      <  2:448122/100‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     <  2:488124/100‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     >  2:1968475/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                 >  2:1134467/1‑100 (MQ=255)
                                                   gATGGTGACGATTACGCGCCATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  1:809208/100‑1 (MQ=255)
                                                   gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  1:1672611/100‑1 (MQ=255)
                                                    aTGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc                                            >  2:1435070/1‑100 (MQ=255)
                                                     tGGTGACGATTACGCTCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                           <  1:1428305/100‑1 (MQ=255)
                                                     tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                           <  2:1318790/100‑1 (MQ=255)
                                                      ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                          <  2:840404/100‑1 (MQ=255)
                                                      ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                          <  2:730759/100‑1 (MQ=255)
                                                      ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                          >  2:1591008/1‑100 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         >  2:724931/1‑100 (MQ=255)
                                                        tGACGATTACGCGGAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGAACGCACTTCTTTTGCGGCTGTCCTCTTTTGCCGCa                                         >  2:1259684/1‑99 (MQ=255)
                                                                 cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                               >  2:1087901/1‑100 (MQ=255)
                                                                  gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                              <  1:890713/100‑1 (MQ=255)
                                                                  gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTTATTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                              <  1:1435079/100‑1 (MQ=255)
                                                                    gcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                            >  2:154454/1‑100 (MQ=255)
                                                                       aTCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                         >  2:791439/1‑100 (MQ=255)
                                                                          cGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTTTTTTTTTAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                      >  2:1300191/1‑100 (MQ=255)
                                                                          cGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                      >  1:88952/1‑100 (MQ=255)
                                                                               ataataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                   >  2:1721850/1‑98 (MQ=255)
                                                                               ataataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                 >  2:1721933/1‑100 (MQ=255)
                                                                                taataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCTTTGcc                >  2:1362092/1‑100 (MQ=255)
                                                                                taataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                >  1:1781172/1‑100 (MQ=255)
                                                                                cttccgatctcAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                <  1:1068015/90‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc               <  2:1494469/100‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                >  2:1310905/1‑99 (MQ=255)
                                                                                      gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCCCTGCTTTTGCGGCTGGCCTGTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          <  2:1116570/100‑1 (MQ=255)
                                                                                      gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          <  2:630408/100‑1 (MQ=255)
                                                                                          cAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccagatcggaag      >  2:1068010/1‑90 (MQ=255)
                                                                                           aGTGTCACGTTTTATTTATCTCTTAAAGCACGCACTGCGTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAATGCCCATTTCTGCTc     >  2:1261946/1‑100 (MQ=255)
                                                                                            gtgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa    >  2:196407/1‑100 (MQ=255)
                                                                                             tgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAt   >  2:1350725/1‑100 (MQ=255)
                                                                                              gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  1:313738/100‑1 (MQ=255)
                                                                                              gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  >  2:662400/1‑100 (MQ=255)
                                                                                                    |                                                                                               
TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281

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