Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,648,042 T→C intergenic (‑254/‑486) dinQ ← / → arsR UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,648,0420TC97.3% 252.6 / ‑2.2 75intergenic (‑254/‑486)dinQ/arsRUV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (41/32);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

TATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG  >  NC_000913/3647945‑3648137
                                                                                                 |                                                                                               
tATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCAt                                                                                               >  1:1676083/1‑100 (MQ=255)
  tGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATgg                                                                                             <  2:1349379/100‑1 (MQ=255)
   gcgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGt                                                                                            >  2:273194/1‑100 (MQ=255)
    cgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTa                                                                                           >  2:1336818/1‑100 (MQ=255)
     gcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTAt                                                                                          >  1:1248928/1‑100 (MQ=255)
     gcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTAt                                                                                          >  1:1467183/1‑100 (MQ=255)
        tctccAGAGTGAAAATGCCTCAGATCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAg                                                                                       <  2:1359908/97‑1 (MQ=255)
          gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                     >  2:773308/1‑100 (MQ=255)
           ccAGAGTGAAAATGCCTCAGATCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAgaga                                                                                    >  1:1359911/1‑97 (MQ=255)
              gagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTg                                                                                  >  1:1554206/1‑99 (MQ=255)
              gagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGa                                                                                 >  1:1114537/1‑100 (MQ=255)
              gagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGa                                                                                 <  1:1816298/100‑1 (MQ=255)
                gTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATg                                                                               <  2:440048/100‑1 (MQ=255)
                 tGAAAATGCCTCAGAGCTAACTGGAGAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGTAATAGTTTCAGCATGGTATTAGCATTGATGc                                                                              <  2:1767536/100‑1 (MQ=255)
                  gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                             >  1:1703507/1‑100 (MQ=255)
                  gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                             <  2:1245581/100‑1 (MQ=255)
                  gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                             >  1:452927/1‑100 (MQ=255)
                    aaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGtt                                                                           <  1:1107483/100‑1 (MQ=255)
                    aaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGtt                                                                           <  1:490981/100‑1 (MQ=255)
                      aTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAg                                                                         <  2:14495/100‑1 (MQ=255)
                      aTGCCTCAGAGCTAACTGGATAATCATAAAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAg                                                                         <  2:1011293/100‑1 (MQ=255)
                       tGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAga                                                                        <  2:1281597/100‑1 (MQ=255)
                        gCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgat                                                                       <  2:511775/100‑1 (MQ=255)
                         ccTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatg                                                                      <  2:524059/100‑1 (MQ=255)
                         ccTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatg                                                                      <  2:2079504/100‑1 (MQ=255)
                           tCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatgat                                                                    >  2:1115499/1‑100 (MQ=255)
                           tCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatgat                                                                    >  1:1369704/1‑100 (MQ=255)
                             agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                  >  2:894148/1‑100 (MQ=255)
                             agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                  <  1:547761/100‑1 (MQ=255)
                             agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                  <  2:1914938/100‑1 (MQ=255)
                                gCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTa                                                               >  1:76256/1‑100 (MQ=255)
                                    aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtctc                                                           >  1:494700/1‑100 (MQ=255)
                                     cTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCa                                                          >  2:1824065/1‑100 (MQ=255)
                                       ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                        >  1:2215764/1‑100 (MQ=255)
                                       ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                        >  2:2081054/1‑100 (MQ=255)
                                       ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                        >  2:1664860/1‑100 (MQ=255)
                                           aaTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAg                                                    <  1:1559205/100‑1 (MQ=255)
                                            ntCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGt                                                   >  2:1235431/2‑100 (MQ=255)
                                            aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGt                                                   >  2:1705779/1‑100 (MQ=255)
                                             tCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTc                                                  >  1:1339981/1‑100 (MQ=255)
                                             tCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTc                                                  <  2:1554202/100‑1 (MQ=255)
                                                 aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag                                              >  1:380229/1‑100 (MQ=255)
                                                 aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag                                              <  2:76254/100‑1 (MQ=255)
                                                 aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag                                              <  2:432597/100‑1 (MQ=255)
                                                 aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag                                              <  1:1212452/100‑1 (MQ=255)
                                                   aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcc                                            >  2:126539/1‑100 (MQ=255)
                                                    gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcca                                           <  1:1137786/100‑1 (MQ=255)
                                                    gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcca                                           >  2:1860219/1‑100 (MQ=255)
                                                         tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                      >  2:169607/1‑100 (MQ=255)
                                                         tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                      >  2:1596386/1‑100 (MQ=255)
                                                           cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTc                                    >  2:1938466/1‑100 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                  >  2:234901/1‑100 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                  >  1:1292055/1‑100 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                  >  2:431367/1‑100 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                  >  2:2229865/1‑100 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                  >  1:777919/1‑100 (MQ=255)
                                                                  taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                             >  1:1613466/1‑100 (MQ=255)
                                                                  taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                             >  2:1764399/1‑100 (MQ=255)
                                                                   aaaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTgg                            >  1:792381/1‑100 (MQ=255)
                                                                    aaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGTTGGCTATCTGACTGCAGTCAGAGCCACCAACTCAGGGCTGGa                           >  2:640625/1‑100 (MQ=39)
                                                                       ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                        <  2:2016042/100‑1 (MQ=255)
                                                                       ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                        <  1:1764404/100‑1 (MQ=255)
                                                                       ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                        <  1:773309/100‑1 (MQ=255)
                                                                           cACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaaa                    <  2:872784/100‑1 (MQ=255)
                                                                             cGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaaaaa                  >  2:2070623/1‑100 (MQ=255)
                                                                                ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc               >  2:640699/1‑100 (MQ=255)
                                                                                ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc               >  2:1720913/1‑100 (MQ=255)
                                                                                ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc               <  2:158549/100‑1 (MQ=255)
                                                                                 cTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCg              <  2:1114535/100‑1 (MQ=255)
                                                                                   tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc            <  1:126541/100‑1 (MQ=255)
                                                                                    gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg           <  2:380227/100‑1 (MQ=255)
                                                                                    gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg           <  1:1336821/100‑1 (MQ=255)
                                                                                          gTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGt     >  1:589968/1‑100 (MQ=255)
                                                                                           tttCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACGCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGActttgcgtc    <  1:1720918/100‑9 (MQ=25)
                                                                                             tCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCgg  >  2:193939/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG  >  NC_000913/3647945‑3648137

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: