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

TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281
                                                                                                     |                                                                                                
ttCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGGTTACGCGCAATACGGCAATAATAGGTTAAAGTGTCAC‑g                                                                                                   >  1:2366344/1‑100 (MQ=255)
  ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑t                                                                                             >  2:264880/1‑100 (MQ=255)
     gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttt                                                                                          >  1:1330021/1‑100 (MQ=255)
      tCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑ttttt                                                                                         <  2:1500892/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttt                                                                                        <  1:1036426/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttt                                                                                        >  1:462145/1‑100 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttt                                                                                        <  2:338569/100‑1 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttt                                                                                        <  1:1974232/99‑1 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑ttttttta                                                                                      <  2:748621/100‑3 (MQ=255)
          ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttttat                                                                                     <  1:141041/100‑4 (MQ=255)
            ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttttatct                                                                                   <  2:1127938/100‑6 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑T‑TTTTTTATCTCTTaaa                                                                              <  2:2065755/100‑1 (MQ=255)
                     tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCTCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgca                                                                           <  2:1168198/100‑1 (MQ=255)
                     tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgca                                                                           <  2:33704/100‑1 (MQ=255)
                      ntttCGATAGGATCACCGTAATATTGCCGGAAGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGAGCCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgca                                                                           >  2:2208132/2‑99 (MQ=255)
                      gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgcac                                                                          >  1:394868/1‑100 (MQ=255)
                       tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgcacg                                                                         <  1:1246962/100‑1 (MQ=255)
                       tttCGATAGGATCACCGTAATATTGCCGGAGGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgcacg                                                                         <  2:1333412/100‑1 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgcacgc                                                                        >  1:288755/1‑100 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTT‑TTTTTTATCTCTTAAAgcacgc                                                                        >  1:821895/1‑99 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:1697237/1‑100 (MQ=255)
                               gcaTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:1624358/3‑100 (MQ=255)
                                 aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACTCACTGCtttt                                                               >  1:865884/1‑100 (MQ=255)
                                  tCACCGTAATTTTGCCGGATGGTTACGATTTCGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCttt                                                                <  1:163375/98‑1 (MQ=255)
                                    aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            <  1:2382925/100‑1 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  2:762420/1‑100 (MQ=255)
                                      cGTAATATTGCCGGAGGGTGTTGATTACGCGCACTCTTTCAATAATAGGTTACATTGTCTCTT‑‑‑T‑TTTTTTATCTCTTTAATCACGCAGTGCTTTTg                                                              <  2:1036092/96‑1 (MQ=255)
                                           tcttGCCTGCTGTTTACGATTACGCGCAATCCTTCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACTTACTTCTTTTGCGGCTGGcc                                                     <  2:884786/98‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     <  2:621358/100‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc                                                      <  1:400636/99‑1 (MQ=255)
                                             ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                   <  2:311680/100‑1 (MQ=255)
                                              tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGAACGCACTGCTTTTGCGGCTGGCctct                                                  >  1:1838646/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTTTTTCTTTTATCACGCACTGCTTTTGCGGCt                                                         >  2:1437973/1‑92 (MQ=255)
                                                 cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                <  2:1236185/99‑1 (MQ=255)
                                                    aTGGTGACGACTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAATCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc                                            <  2:1573989/100‑1 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTACGGCTGGCCTCTTTTGACGCa                                         >  1:370572/1‑100 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGATACAGTGGCACGG‑‑‑T‑TTTTGTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         >  1:2281700/1‑100 (MQ=255)
                                                       gTGACGATTACGCGCAATCAGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         <  2:650361/100‑1 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaaa                                      >  1:1213753/1‑99 (MQ=255)
                                                             aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     <  1:1143039/98‑1 (MQ=255)
                                                               tACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTc                                 <  1:264883/100‑1 (MQ=255)
                                                                aCGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                <  1:1555146/100‑1 (MQ=255)
                                                                      aaTCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATATCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTAGCCCgtg                          >  2:224061/1‑100 (MQ=255)
                                                                       aTCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                         >  2:1724079/1‑100 (MQ=255)
                                                                       aTCCGCCAATAATAGGTTACAGTCCCAC‑G‑‑‑TCTATCTTAACTCTTAAAGCACCCACTCCTCATGCGGCTAGCCTCTTTTGCCGCCACATatccgc                               >  2:2223249/1‑89 (MQ=255)
                                                                        tCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                        <  1:413462/100‑1 (MQ=255)
                                                                          cGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                      >  1:1246367/1‑100 (MQ=255)
                                                                               ataataGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                     <  1:221997/96‑1 (MQ=255)
                                                                                taataGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                <  2:288752/100‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCAC‑G‑‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa              >  2:401270/1‑100 (MQ=255)
                                                                                    aGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGTCTCTTTTGCTGCAAAATAGTCGCCCGTGTTTCATTGCCCAtt            >  2:1686521/1‑100 (MQ=255)
                                                                                     ggTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTTTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt           >  2:1890481/1‑100 (MQ=255)
                                                                                      gTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          >  2:806650/1‑100 (MQ=255)
                                                                                       ttACTGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTGTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCt         <  1:35568/100‑1 (MQ=255)
                                                                                       ttACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCt         >  1:1909645/1‑100 (MQ=255)
                                                                                         aCAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc       <  1:289537/100‑1 (MQ=255)
                                                                                              gtCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  1:1681953/100‑1 (MQ=255)
                                                                                               tCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  1:617683/99‑1 (MQ=255)
                                                                                                     |                                                                                                
TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  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: ‑