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.9% 98.6 / 4.9 33intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base T (13/17);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

GTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCAT  >  NC_000913/1665094‑1665284
                                                                                               |                                                                                                
gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttt                                                                                             >  1:4565/1‑100 (MQ=255)
 tCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtttttt                                                                                            <  2:746366/100‑1 (MQ=255)
     ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCCATAATAGGTTACAGTGTCCCGttttttt                                                                                           >  2:1392831/1‑97 (MQ=255)
     ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttat                                                                                        >  2:226783/1‑97 (MQ=255)
     ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGAGGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttat                                                                                        <  2:184802/100‑4 (MQ=255)
       ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTAtctc                                                                                      <  2:1268894/100‑1 (MQ=255)
          aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCtt                                                                                    <  2:293821/100‑1 (MQ=255)
           aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTa                                                                                   >  2:72965/1‑100 (MQ=255)
            ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaa                                                                                  >  2:1514417/1‑100 (MQ=255)
              gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAg                                                                                <  1:1514426/100‑1 (MQ=255)
               cTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgc                                                                               >  2:134773/1‑100 (MQ=255)
               cTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAgca                                                                              <  1:858764/100‑1 (MQ=255)
                  tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacg                                                                            <  1:1110948/100‑1 (MQ=255)
                      gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACt                                                                        >  1:1300097/1‑100 (MQ=255)
                           gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCttn                                                                   <  2:1134573/100‑2 (MQ=255)
                              cACCGTAATATTGCCGGATGGTGACGATTACGCGTATTCCGGTAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                                <  2:1107455/100‑1 (MQ=255)
                               aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                               >  2:1310383/1‑100 (MQ=255)
                                ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                              >  1:1519772/1‑100 (MQ=255)
                                ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                               >  1:872395/1‑99 (MQ=255)
                                 cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                             <  1:257529/100‑1 (MQ=255)
                                         tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctcn                                                     <  2:300208/100‑2 (MQ=255)
                                              gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                                <  1:281580/100‑1 (MQ=255)
                                                  gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                            <  1:929110/100‑1 (MQ=255)
                                                      cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                        >  2:1200325/1‑100 (MQ=255)
                                                        aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAg                                      <  2:1388373/100‑1 (MQ=255)
                                                             gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                                 >  1:948941/1‑100 (MQ=255)
                                                               gcAATCCGTCAATAATAGGTTACAGTGTCACTTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                               <  2:280930/100‑1 (MQ=255)
                                                                  aTCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                            >  2:659027/1‑100 (MQ=255)
                                                                     ctgCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                         >  2:1085021/3‑100 (MQ=255)
                                                                       gCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                       <  2:93483/100‑1 (MQ=255)
                                                                          ataataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                    <  1:1298214/100‑1 (MQ=255)
                                                                             ataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCTTTGCCCa                 <  1:828429/100‑1 (MQ=255)
                                                                             ataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa                 <  1:72965/100‑1 (MQ=255)
                                                                              taGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAt                >  1:929663/1‑100 (MQ=255)
                                                                                 tccgatctgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCCAAATAGTCGCCCGTGTTTCATTGCCCATTTc             <  1:1063016/93‑1 (MQ=255)
                                                                                      aGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTc        <  2:452379/100‑1 (MQ=255)
                                                                                         gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAGATAGTGGCCCGTGTGTCATTGCCCATTTcagatcgga     >  2:1063010/1‑92 (MQ=255)
                                                                                           cACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGcat  >  1:40512/1‑100 (MQ=255)
                                                                                               |                                                                                                
GTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCAT  >  NC_000913/1665094‑1665284

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