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

TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑TT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281
                                                                                                    |                                                                                               
ttCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                                                                  <  2:948319/100‑1 (MQ=255)
    tGTCGCGGTTAAACCGCTGGTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGGTACAGTGTCACG‑tt‑tt                                                                                            >  2:809205/1‑100 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑tt‑tt‑tttta                                                                                      <  2:1640838/100‑3 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑tt‑tt‑tttta                                                                                      <  2:1668122/100‑3 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑tt‑tt‑ttttatc                                                                                    <  1:587211/100‑5 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑tt‑tt‑ttttatc                                                                                    >  1:737308/1‑96 (MQ=255)
            ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑tt‑tt‑ttttatct                                                                                   >  2:1361023/1‑95 (MQ=255)
            ttAAACCGCTGTTTCGATAGGATCACCGTAATATCGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑tt‑tt‑ttttatct                                                                                   <  2:1080368/100‑6 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCtt                                                                                 <  2:733477/100‑1 (MQ=255)
                aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTa                                                                                <  2:1619164/100‑1 (MQ=255)
                aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑TT‑TT‑TTTATCTCTTaa                                                                               >  1:1494555/1‑100 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCGGGATGGTGACGATTACGCGAAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTaa                                                                               <  2:64567/100‑1 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑TT‑TT‑TTTATCTCTTaaa                                                                              <  2:1494546/100‑1 (MQ=255)
                      gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAg                                                                             >  1:827997/1‑97 (MQ=255)
                       tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAgcacg                                                                         <  1:520440/100‑1 (MQ=255)
                       tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑TT‑TT‑TTTATCTCTTAAAgcacgc                                                                        <  1:1841843/100‑1 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAgcacgc                                                                        <  2:1768148/100‑1 (MQ=255)
                         tCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAgcacgca                                                                       >  2:1846665/1‑100 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGCCACGTTT‑TT‑TTTATCTCTTAAAgcacgcac                                                                      <  2:1794126/100‑1 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACt                                                                     <  1:1120005/100‑1 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACt                                                                     >  1:31195/1‑100 (MQ=255)
                            aTAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTg                                                                    >  1:2040248/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:1494478/1‑100 (MQ=255)
                                gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCttt                                                                >  2:1520767/1‑100 (MQ=255)
                                   cACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGc                                                             >  1:621724/1‑100 (MQ=255)
                                    aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            <  2:1842932/100‑1 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  1:1710695/1‑100 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  2:136482/1‑100 (MQ=255)
                                      cGTAATATGGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                          <  1:1403948/100‑1 (MQ=255)
                                       gTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCt                                                         <  2:1965979/100‑1 (MQ=255)
                                       gTAATATTGCCGGATGGTGACGATTACGCGCAATCCAGCAATAATAGGTTACAGTGTCACG‑TT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        >  1:1441400/1‑100 (MQ=255)
                                        tAATATTGCCGGTTGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        <  2:529694/100‑1 (MQ=255)
                                        tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        <  1:1853352/100‑1 (MQ=255)
                                         aaTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                       >  1:840407/1‑100 (MQ=255)
                                          atatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc                                                      <  2:448122/100‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     <  2:488124/100‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     >  2:1968475/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                 >  2:1134467/1‑100 (MQ=255)
                                                   gATGGTGACGATTACGCGCCATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  1:809208/100‑1 (MQ=255)
                                                   gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  1:1672611/100‑1 (MQ=255)
                                                    aTGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc                                            >  2:1435070/1‑100 (MQ=255)
                                                     tGGTGACGATTACGCTCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                           <  1:1428305/100‑1 (MQ=255)
                                                     tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                           <  2:1318790/100‑1 (MQ=255)
                                                      ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                          <  2:840404/100‑1 (MQ=255)
                                                      ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                          <  2:730759/100‑1 (MQ=255)
                                                      ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                          >  2:1591008/1‑100 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         >  2:724931/1‑100 (MQ=255)
                                                        tGACGATTACGCGGAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGAACGCACTTCTTTTGCGGCTGTCCTCTTTTGCCGCa                                         >  2:1259684/1‑99 (MQ=255)
                                                                 cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                               >  2:1087901/1‑100 (MQ=255)
                                                                  gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                              <  1:890713/100‑1 (MQ=255)
                                                                  gcgcAATCCGGCAATAATAGGTTACAGTGTCACG‑TTATT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                              <  1:1435079/100‑1 (MQ=255)
                                                                    gcAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                            >  2:154454/1‑100 (MQ=255)
                                                                       aTCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                         >  2:791439/1‑100 (MQ=255)
                                                                          cGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTTTTTTTTTTAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                      >  2:1300191/1‑100 (MQ=255)
                                                                          cGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                      >  1:88952/1‑100 (MQ=255)
                                                                               ataataGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                   >  2:1721850/1‑98 (MQ=255)
                                                                               ataataGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                 >  2:1721933/1‑100 (MQ=255)
                                                                                taataGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCTTTGcc                >  2:1362092/1‑100 (MQ=255)
                                                                                taataGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                >  1:1781172/1‑100 (MQ=255)
                                                                                cttccgatctcAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                <  1:1068015/90‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc               <  2:1494469/100‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                >  2:1310905/1‑99 (MQ=255)
                                                                                      gTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCCCTGCTTTTGCGGCTGGCCTGTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          <  2:1116570/100‑1 (MQ=255)
                                                                                      gTTACAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          <  2:630408/100‑1 (MQ=255)
                                                                                          cAGTGTCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccagatcggaag      >  2:1068010/1‑90 (MQ=255)
                                                                                           aGTGTCACG‑TT‑TTATTTATCTCTTAAAGCACGCACTGCGTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAATGCCCATTTCTGCTc     >  2:1261946/1‑100 (MQ=255)
                                                                                            gtgtCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa    >  2:196407/1‑100 (MQ=255)
                                                                                             tgtCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAt   >  2:1350725/1‑100 (MQ=255)
                                                                                              gtCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  1:313738/100‑1 (MQ=255)
                                                                                              gtCACGTTT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  >  2:662400/1‑100 (MQ=255)
                                                                                                    |                                                                                               
TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑TT‑TT‑TTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  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: ‑