Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,809,581 G→A 76.9% G52G (GGG→GGA mprA → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,809,5810GA76.9% 406.5 / 132.1 208G52G (GGG→GGAmprADNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base A (68/92);  total (90/118)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 209 total aligned reads displayed.

CTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCT  >  W3110S.gb/2809514‑2809647
                                                                   |                                                                  
cTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAAtt                                                                 <  1:773471/71‑1 (MQ=255)
cTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAAtt                                                                 <  1:283003/71‑1 (MQ=255)
 tCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAAtt                                                                 >  1:2073992/1‑70 (MQ=255)
   gTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                                              <  1:2109658/71‑1 (MQ=255)
   gTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTa                                                                <  1:2061444/69‑1 (MQ=255)
   gTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAAc                                                              <  1:617366/71‑1 (MQ=255)
   gTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAAc                                                              <  1:3414917/71‑1 (MQ=255)
   gTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAAc                                                              <  1:1974511/71‑1 (MQ=255)
   gTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAAc                                                              <  1:2344843/71‑1 (MQ=255)
     cTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                                              >  1:1743760/1‑69 (MQ=255)
      tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAAc                                                              >  1:3677932/1‑68 (MQ=255)
        tgcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACga                                                            <  1:2175088/68‑1 (MQ=255)
        tgcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACga                                                            <  1:2439270/68‑1 (MQ=255)
        tgcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGAc                                                         >  1:1628852/1‑71 (MQ=255)
         gcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGAc                                                         <  1:2976166/70‑1 (MQ=255)
         gcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGAc                                                         <  1:1400996/70‑1 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                                   >  1:2799555/1‑71 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                                   >  1:339755/1‑71 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTTACGAGACgttgtt                                                   >  1:1254579/1‑71 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACgaga                                                          >  1:1227161/1‑64 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:1839838/71‑1 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:874993/71‑1 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:2078688/71‑1 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:2035704/71‑1 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:217751/71‑1 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   >  1:2278742/1‑71 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   >  1:1505172/1‑71 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   >  1:1391418/1‑71 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:1053798/71‑1 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACAGCAATAAAATGGTGAAGGCTCAGGGAGTTAACGAGACgttgtt                                                   >  1:1344732/1‑71 (MQ=255)
                  tgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                                     <  1:1928055/65‑1 (MQ=255)
                  tgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                                     <  1:1587189/65‑1 (MQ=255)
                  tgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                                     <  1:1546950/65‑1 (MQ=255)
                  tgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAAtt                                                                 <  1:3199340/53‑1 (MQ=255)
                  tgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACga                                                            <  1:2607405/58‑1 (MQ=255)
                  tgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACga                                                            <  1:144959/58‑1 (MQ=255)
                   gcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATgg                                              >  1:3141316/1‑71 (MQ=255)
                   gcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATgg                                              >  1:1689617/1‑71 (MQ=255)
                    caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACgaga                                                          >  1:3127809/1‑58 (MQ=255)
                    caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACgaga                                                          >  1:2862455/1‑58 (MQ=255)
                    caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgt                                                    >  1:3333285/1‑64 (MQ=255)
                    caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTa                                                 >  1:1069971/1‑67 (MQ=255)
                    caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTa                                                 >  1:718790/1‑67 (MQ=255)
                    caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGc                                             >  1:2480100/1‑71 (MQ=255)
                      aagcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                                   <  1:2988514/63‑1 (MQ=255)
                       agcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                          >  1:1904284/1‑71 (MQ=255)
                       agcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGtt                                          <  1:358793/71‑1 (MQ=255)
                       agcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGtt                                          <  1:3494568/71‑1 (MQ=255)
                       agcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGtt                                          <  1:1274901/71‑1 (MQ=255)
                       agcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGtt                                          <  1:332056/71‑1 (MQ=255)
                       agcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGtt                                          <  1:292050/71‑1 (MQ=255)
                        gcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTa                                                 >  1:2117028/1‑63 (MQ=255)
                         caagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgtt                                                      >  1:2388613/1‑57 (MQ=255)
                         caagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTa                                                 >  1:427277/1‑62 (MQ=255)
                         caagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTg                                         <  1:3311840/70‑1 (MQ=255)
                          aagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                                       <  1:3544094/71‑1 (MQ=255)
                          aagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAt                                       <  1:2743381/71‑1 (MQ=255)
                          aagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAt                                       <  1:3514064/71‑1 (MQ=255)
                           agcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAt                                       <  1:3270831/70‑1 (MQ=255)
                                 tgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACgaga                                                          <  1:3228756/45‑1 (MQ=255)
                                 tgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTa                                     <  1:450000/66‑1 (MQ=255)
                                 tgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTa                                     <  1:679514/66‑1 (MQ=255)
                                 tgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTa                                     <  1:2761308/66‑1 (MQ=255)
                                  ggAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:3611738/51‑1 (MQ=255)
                                  ggAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:1098962/51‑1 (MQ=255)
                                  ggAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg                                >  1:2197247/1‑70 (MQ=255)
                                   gagaACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg                               <  1:1323486/70‑1 (MQ=255)
                                   gagaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg                               <  1:28926/70‑1 (MQ=255)
                                   gagaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg                               <  1:1755421/70‑1 (MQ=255)
                                   gagaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg                               <  1:211993/70‑1 (MQ=255)
                                      aaCCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtc                           <  1:2453213/71‑1 (MQ=255)
                                      aaCCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtc                           <  1:1825032/71‑1 (MQ=255)
                                       aCCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtc                           <  1:1236850/70‑1 (MQ=255)
                                             aTAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                   <  1:2676261/40‑1 (MQ=255)
                                             aTAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                    <  1:3645348/71‑1 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                        >  1:3442758/1‑50 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                        >  1:487521/1‑50 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                        >  1:3555496/1‑50 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                          <  1:1720230/46‑1 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                 >  1:654025/1‑71 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                 >  1:1419156/1‑71 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                 >  1:691927/1‑71 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                 >  1:2215955/1‑71 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                 >  1:533590/1‑71 (MQ=255)
                                                  aTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTa                                                 >  1:2728572/1‑37 (MQ=255)
                                                   tGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGtt                                          <  1:1931938/43‑1 (MQ=255)
                                                   tGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACca               >  1:2697551/1‑70 (MQ=255)
                                                    gCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaaa                  >  1:3018949/1‑66 (MQ=255)
                                                         aGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTAtt        >  1:3067957/1‑71 (MQ=255)
                                                         aGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTAtt        >  1:208621/1‑71 (MQ=255)
                                                         aGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTAt         >  1:1807168/1‑70 (MQ=255)
                                                         aGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTAt         >  1:3161100/1‑70 (MQ=255)
                                                            cTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAg     <  1:82175/71‑1 (MQ=255)
                                                            cTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAg     <  1:2136662/71‑1 (MQ=255)
                                                             tCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTa                                     >  1:575406/1‑38 (MQ=255)
                                                             tCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGc    >  1:1179264/1‑71 (MQ=255)
                                                               aGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGcc   <  1:2555287/70‑1 (MQ=255)
                                                               aGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGcc   <  1:2515530/70‑1 (MQ=255)
                                                               aGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCt  <  1:3706318/71‑1 (MQ=255)
                                                                gggaATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAggc                           <  1:3013495/45‑3 (MQ=38)
                                                                   |                                                                  
CTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCT  >  W3110S.gb/2809514‑2809647

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