Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,665,195:1 (T)7→8 100% 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.T89.4% 139.1 / 11.8 47intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base T (19/23);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281
                                                                                                     |                                                                                                
ttCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGGTTACGCGCAATACGGCAATAATAGGTTAAAGTGTCACg                                                                                                   >  1:2366344/1‑100 (MQ=255)
  ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑tt                                                                                             >  2:264880/1‑100 (MQ=255)
     gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑ttttt                                                                                          >  1:1330021/1‑100 (MQ=255)
      tCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑tttttt                                                                                         <  2:1500892/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑ttttttt                                                                                        <  1:1036426/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑ttttttt                                                                                        >  1:462145/1‑100 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑ttttttt                                                                                        <  2:338569/100‑1 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑ttttttt                                                                                        <  1:1974232/99‑1 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑tttttttta                                                                                      <  2:748621/100‑3 (MQ=255)
          ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑ttttttttat                                                                                     <  1:141041/100‑4 (MQ=255)
            ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑ttttttttatct                                                                                   <  2:1127938/100‑6 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑TTTTTTTATCTCTTaaa                                                                              <  2:2065755/100‑1 (MQ=255)
                     tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCTCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAgca                                                                           <  2:1168198/100‑1 (MQ=255)
                     tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAgca                                                                           <  2:33704/100‑1 (MQ=255)
                      ntttCGATAGGATCACCGTAATATTGCCGGAAGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGAGCCACGT‑‑TTTTTTTATCTCTTAAAgca                                                                           >  2:2208132/2‑99 (MQ=255)
                      gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAgcac                                                                          >  1:394868/1‑100 (MQ=255)
                       tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAgcacg                                                                         <  1:1246962/100‑1 (MQ=255)
                       tttCGATAGGATCACCGTAATATTGCCGGAGGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAgcacg                                                                         <  2:1333412/100‑1 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAgcacgc                                                                        >  1:288755/1‑100 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTTATCTCTTAAAgcacgc                                                                        >  1:821895/1‑99 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:1697237/1‑100 (MQ=255)
                               gcaTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:1624358/3‑100 (MQ=255)
                                 aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACTCACTGCtttt                                                               >  1:865884/1‑100 (MQ=255)
                                  tCACCGTAATTTTGCCGGATGGTTACGATTTCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCttt                                                                <  1:163375/98‑1 (MQ=255)
                                    aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            <  1:2382925/100‑1 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  2:762420/1‑100 (MQ=255)
                                      cGTAATATTGCCGGAGGGTGTTGATTACGCGCACTCTTTCAATAATAGGTTACATTGTCTCTT‑‑‑TTTTTTTATCTCTTTAATCACGCAGTGCTTTTg                                                              <  2:1036092/96‑1 (MQ=255)
                                           tcttGCCTGCTGTTTACGATTACGCGCAATCCTTCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACTTACTTCTTTTGCGGCTGGcc                                                     <  2:884786/98‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     <  2:621358/100‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc                                                      <  1:400636/99‑1 (MQ=255)
                                             ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                   <  2:311680/100‑1 (MQ=255)
                                              tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGAACGCACTGCTTTTGCGGCTGGCctct                                                  >  1:1838646/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTTTTTCTTTTATCACGCACTGCTTTTGCGGCt                                                         >  2:1437973/1‑92 (MQ=255)
                                                 cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                <  2:1236185/99‑1 (MQ=255)
                                                    aTGGTGACGACTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAATCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc                                            <  2:1573989/100‑1 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTACGGCTGGCCTCTTTTGACGCa                                         >  1:370572/1‑100 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGATACAGTGGCACGG‑‑‑TTTTTGTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         >  1:2281700/1‑100 (MQ=255)
                                                       gTGACGATTACGCGCAATCAGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         <  2:650361/100‑1 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaaa                                      >  1:1213753/1‑99 (MQ=255)
                                                             aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     <  1:1143039/98‑1 (MQ=255)
                                                               tACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTc                                 <  1:264883/100‑1 (MQ=255)
                                                                aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                <  1:1555146/100‑1 (MQ=255)
                                                                      aaTCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATATCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTAGCCCgtg                          >  2:224061/1‑100 (MQ=255)
                                                                       aTCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                         >  2:1724079/1‑100 (MQ=255)
                                                                       aTCCGCCAATAATAGGTTACAGTCCCACG‑‑‑TCTATCTTAACTCTTAAAGCACCCACTCCTCATGCGGCTAGCCTCTTTTGCCGCCACATatccgc                               >  2:2223249/1‑89 (MQ=255)
                                                                        tCCGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                        <  1:413462/100‑1 (MQ=255)
                                                                          cGGCAATAATAGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                      >  1:1246367/1‑100 (MQ=255)
                                                                               ataataGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                     <  1:221997/96‑1 (MQ=255)
                                                                                taataGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                <  2:288752/100‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCACG‑‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa              >  2:401270/1‑100 (MQ=255)
                                                                                    aGGTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGTCTCTTTTGCTGCAAAATAGTCGCCCGTGTTTCATTGCCCAtt            >  2:1686521/1‑100 (MQ=255)
                                                                                     ggTTACAGTGTCACGT‑‑TTTTTTTATCTTTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt           >  2:1890481/1‑100 (MQ=255)
                                                                                      gTTACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          >  2:806650/1‑100 (MQ=255)
                                                                                       ttACTGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTGTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCt         <  1:35568/100‑1 (MQ=255)
                                                                                       ttACAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCt         >  1:1909645/1‑100 (MQ=255)
                                                                                         aCAGTGTCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc       <  1:289537/100‑1 (MQ=255)
                                                                                              gtCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  1:1681953/100‑1 (MQ=255)
                                                                                               tCACGT‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  1:617683/99‑1 (MQ=255)
                                                                                                     |                                                                                                
TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281

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