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,0420TC98.2% 194.1 / NA 56intergenic (‑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 (0/0);  major base C (26/29);  minor base A (1/0);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGT  >  NC_000913/3647946‑3648138
                                                                                                |                                                                                                
aTGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATg                                                                                               <  2:592044/100‑1 (MQ=255)
aTGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATg                                                                                               >  1:798314/1‑100 (MQ=255)
 tGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATgg                                                                                              >  1:1070725/1‑100 (MQ=255)
     cATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTAtt                                                                                          <  2:181270/100‑1 (MQ=255)
        tGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGc                                                                                       >  1:948116/1‑100 (MQ=255)
           cAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAtt                                                                                    <  1:811602/100‑1 (MQ=255)
            agagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTg                                                                                   <  2:869942/100‑1 (MQ=255)
              agTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGAt                                                                                 >  1:315021/1‑100 (MQ=255)
               gTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATg                                                                                <  1:742234/100‑1 (MQ=255)
                 gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                              >  1:122889/1‑100 (MQ=255)
                 gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                              <  1:893941/100‑1 (MQ=255)
                  aaaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGt                                                                             >  2:326474/1‑100 (MQ=255)
                   aaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGtt                                                                            <  2:1119694/100‑1 (MQ=255)
                     aTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAAGCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGAGGCGTTAg                                                                          <  2:1254778/100‑1 (MQ=255)
                      tGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAga                                                                         >  2:507022/1‑100 (MQ=255)
                         cTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatga                                                                      >  1:997087/1‑100 (MQ=255)
                           cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                    <  2:199431/100‑1 (MQ=255)
                               gCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTGGCAATAGTTGCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTa                                                                <  2:1098868/100‑1 (MQ=255)
                                cTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAt                                                               >  2:1117838/1‑100 (MQ=255)
                                   aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtctc                                                            <  2:1204843/100‑1 (MQ=255)
                                    cTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCa                                                           <  2:228231/100‑1 (MQ=255)
                                       gATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTc                                                        <  1:973999/100‑1 (MQ=255)
                                       gATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTc                                                        >  2:539443/1‑100 (MQ=255)
                                       gATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTc                                                        >  1:391124/1‑100 (MQ=255)
                                        aTAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTcc                                                       <  1:913207/100‑1 (MQ=255)
                                         tAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCa                                                      >  2:623967/1‑100 (MQ=255)
                                           aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGt                                                    >  1:756481/1‑100 (MQ=255)
                                             cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa                                                  <  1:1019054/100‑1 (MQ=255)
                                             cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa                                                  >  2:341853/1‑100 (MQ=255)
                                              aTACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCag                                                 <  2:170810/100‑1 (MQ=255)
                                                aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag                                               <  1:850632/100‑1 (MQ=255)
                                                aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCaga                                                <  1:1145603/99‑1 (MQ=255)
                                                    tACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccac                                           <  1:1261341/100‑1 (MQ=255)
                                                     aCATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacc                                          <  1:268704/100‑1 (MQ=255)
                                                      cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                         <  1:1168713/100‑1 (MQ=255)
                                                      cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTAAGCATGGTATTAGCATTGATGCGTTAGATGATGGCGATCTCACTCCAGTCAGAGccacc                                          >  1:567721/1‑99 (MQ=255)
                                                        tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                       <  1:517777/100‑1 (MQ=255)
                                                        tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                       >  1:756908/1‑100 (MQ=255)
                                                        tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                       >  2:1073231/1‑100 (MQ=255)
                                                         gCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACt                                      >  1:531558/1‑100 (MQ=255)
                                                         gCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACt                                      >  1:535034/1‑100 (MQ=255)
                                                         gCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGAATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCATCAACt                                      >  2:779032/1‑100 (MQ=255)
                                                              ttATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACTAACTCAggg                                 >  1:1267178/1‑100 (MQ=255)
                                                               tataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                                >  2:446460/1‑100 (MQ=255)
                                                                 taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                              >  2:581897/1‑100 (MQ=255)
                                                                   aaaCCAGCACGTCCTTGCAATAGTTGCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                              <  2:723533/98‑1 (MQ=255)
                                                                      ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                         <  1:507026/100‑1 (MQ=255)
                                                                      ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                         <  2:249147/100‑1 (MQ=255)
                                                                        aGCACGTCCTTGCAATAGTTTCAGCATGGTAGTAGCATGGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGt                        <  2:1290298/99‑1 (MQ=255)
                                                                          cACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaaa                     <  2:463899/100‑1 (MQ=255)
                                                                               ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc                >  2:550410/1‑100 (MQ=255)
                                                                                atctcAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc                <  2:715615/95‑1 (MQ=255)
                                                                                    cAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACga           >  2:429832/1‑99 (MQ=255)
                                                                                    cAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGc           >  2:1011823/1‑100 (MQ=255)
                                                                                    cAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAAGc           >  1:715620/1‑100 (MQ=255)
                                                                                             cAGCATGGTATTAGCATTGATGCGTTAGATGCTGGCTATCTCACTCCAGTCAGAGCCACAAACTCAGGGCTGGAACGTAAAAAACCGACGCAAAGTCGGt  <  2:1088045/100‑1 (MQ=21)
                                                                                                |                                                                                                
ATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGT  >  NC_000913/3647946‑3648138

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: