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.T95.3% 142.1 / ‑0.1 43intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (22/19);  total (23/20)
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.00e+00

TGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATC  >  NC_000913/1665093‑1665285
                                                                                                 |                                                                                                 
tGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑tttt                                                                                               <  2:920741/100‑1 (MQ=255)
   cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATATGTTACAGTGTCAC‑G‑ttttttt                                                                                            <  2:1209978/100‑1 (MQ=255)
   cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑ttttttt                                                                                            >  2:215701/1‑100 (MQ=255)
   cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑ttttttt                                                                                            <  2:327479/100‑1 (MQ=255)
   cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑ttttttt                                                                                            <  2:104920/100‑1 (MQ=255)
      ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑ttttttttat                                                                                         >  1:1121984/1‑97 (MQ=255)
       gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑ttttttttatc                                                                                        >  1:1081556/1‑96 (MQ=255)
       gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑ttttttttatc                                                                                        >  1:297322/1‑96 (MQ=255)
       gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑ttttttttatc                                                                                        <  2:734773/100‑5 (MQ=255)
         tAAACCGCTGTTTCGTTAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTAtctc                                                                                       <  1:1250282/100‑1 (MQ=255)
         tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTAtctc                                                                                       <  1:1111697/100‑1 (MQ=255)
         tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑tttttttttatct                                                                                      >  2:849878/1‑94 (MQ=255)
          aaaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATctct                                                                                      <  1:1111875/100‑1 (MQ=255)
             ccGCTGTTTCGATAGGATCACCGTATTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCtt                                                                                     <  2:553688/98‑1 (MQ=255)
              cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTaaa                                                                                  >  2:435768/1‑100 (MQ=255)
                 tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAgca                                                                               <  1:1086646/100‑1 (MQ=255)
                  gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAgcac                                                                              >  1:883227/1‑100 (MQ=255)
                   tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAgcacg                                                                             >  2:1396535/1‑100 (MQ=255)
                           ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                     >  1:329255/1‑100 (MQ=255)
                                aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                                <  1:1185331/100‑1 (MQ=255)
                                    tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                            >  1:1428968/1‑100 (MQ=255)
                                       tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                         <  2:1183265/100‑1 (MQ=255)
                                        atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                        >  2:91552/1‑100 (MQ=255)
                                        atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                        >  2:1465126/1‑100 (MQ=255)
                                         ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                       <  1:435770/100‑1 (MQ=255)
                                             cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtttt                                                   >  2:1519591/1‑100 (MQ=255)
                                             cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                    <  2:607461/99‑1 (MQ=255)
                                              ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                  >  2:470131/1‑100 (MQ=255)
                                                  ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                              <  1:1314829/100‑1 (MQ=255)
                                                       cGATTACACGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                         >  1:827947/1‑100 (MQ=255)
                                                            aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGG‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                    <  2:1249602/100‑1 (MQ=255)
                                                             cgcgCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                                   >  2:596972/1‑100 (MQ=255)
                                                              gcgcAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                                  >  1:204833/1‑100 (MQ=255)
                                                                   aTCCGGCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                             <  1:560722/100‑1 (MQ=255)
                                                                        gCAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                        >  2:1444489/1‑100 (MQ=255)
                                                                         cAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                       >  2:236489/1‑100 (MQ=255)
                                                                         cAATAATAGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                       >  2:836496/1‑100 (MQ=255)
                                                                           ataataGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                     >  1:1093134/1‑100 (MQ=255)
                                                                           ataataGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                     >  1:908461/1‑100 (MQ=255)
                                                                             aataGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                   <  2:1011896/100‑1 (MQ=255)
                                                                             aataGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                   >  2:220580/1‑100 (MQ=255)
                                                                              ataGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa                  <  1:1396545/100‑1 (MQ=255)
                                                                               taGGTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTTTTTCATTGCCCAt                 <  2:307436/100‑1 (MQ=255)
                                                                                 ggTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt               >  2:690883/1‑100 (MQ=255)
                                                                                 ggTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt               >  2:762188/1‑100 (MQ=255)
                                                                                  gTTACAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc              <  2:888358/100‑1 (MQ=255)
                                                                                     aCAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg            <  1:470134/99‑1 (MQ=255)
                                                                                     aCAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc           <  1:91552/100‑1 (MQ=255)
                                                                                     aCAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc           <  2:1495453/100‑1 (MQ=255)
                                                                                      cAGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCt          >  2:1225020/1‑100 (MQ=255)
                                                                                       aGTGTCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTc         <  2:329253/100‑1 (MQ=255)
                                                                                          gtCAC‑GTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg      >  2:538447/1‑100 (MQ=255)
                                                                                             aC‑G‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGcttc  >  2:923346/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
TGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATC  >  NC_000913/1665093‑1665285

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