Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,809,581 G→A 55.8% 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,5810GA55.8% 23.9 / 120.9 86G52G (GGG→GGAmprADNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base A (24/24);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTT  >  W3110S.gb/2809514‑2809648
                                                                   |                                                                   
cTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGAtt                                                                  <  1:981777/71‑1 (MQ=255)
    tCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACg                                                              <  1:2263159/71‑1 (MQ=255)
      tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATAAAc                                                               >  1:555719/1‑68 (MQ=255)
      tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAAc                                                               >  1:1144491/1‑68 (MQ=255)
       ttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACgaga                                                           >  1:553149/1‑71 (MQ=255)
         gcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                                         >  1:1108752/1‑71 (MQ=255)
         gcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                                         >  1:628918/1‑71 (MQ=255)
             gcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgt                                                     >  1:1871685/1‑71 (MQ=255)
             gcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgt                                                     >  1:1597927/1‑71 (MQ=255)
             gcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgt                                                     >  1:2221692/1‑71 (MQ=255)
              cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                                    <  1:172931/71‑1 (MQ=255)
               acatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                                    >  1:1323928/1‑70 (MQ=255)
               acatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGttt                                                   <  1:1218010/71‑1 (MQ=255)
                   gcaAAGCAAGCTGCTGTAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTa                                                  >  1:906099/1‑68 (MQ=255)
                   gcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACgaga                                                           >  1:730540/1‑59 (MQ=255)
                   gcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                                      >  1:2146348/1‑64 (MQ=255)
                   gcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATgg                                               >  1:2399935/1‑71 (MQ=255)
                   gcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTa                                                  >  1:1651167/1‑68 (MQ=255)
                   gcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATg                                                >  1:397538/1‑70 (MQ=255)
                    caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTAt                                                 >  1:2424891/1‑68 (MQ=255)
                    caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGc                                              >  1:1510491/1‑71 (MQ=255)
                     aaaGCAAGCTTCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTAt                                                 >  1:505145/1‑67 (MQ=255)
                        gcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTa                                                                 >  1:2569813/1‑48 (MQ=255)
                        gcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgt                                                        >  1:2477165/1‑57 (MQ=255)
                        gcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGtt                                           <  1:1610675/70‑1 (MQ=255)
                         caagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAAtt                                                                  <  1:1412858/46‑1 (MQ=255)
                         caagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACg                                                              >  1:1923088/1‑50 (MQ=255)
                           agcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACTAGACgttgtt                                                    <  1:465684/58‑1 (MQ=255)
                            gctgctGGAGAACCGCAATAAAATGCTGAAGGCTCATGGAAtt                                                                  <  1:4217/43‑1 (MQ=38)
                              tgctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACg                                    >  1:2314301/1‑71 (MQ=255)
                                ctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                                               >  1:636159/1‑42 (MQ=255)
                                ctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACg                                    <  1:798937/69‑1 (MQ=255)
                                 tgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTAc                                     <  1:2024243/67‑1 (MQ=255)
                                   gagaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg                                 <  1:2254296/69‑1 (MQ=255)
                                   gagaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg                                 <  1:1356967/69‑1 (MQ=255)
                                   gagaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                               <  1:1533531/71‑1 (MQ=255)
                                   gagaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                               <  1:267926/71‑1 (MQ=255)
                                     gaACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGt                             >  1:140166/1‑71 (MQ=255)
                                     gaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGAc                                                          >  1:225469/1‑42 (MQ=255)
                                      aaCCGCAATAAAATGCTGAAGGCTCAGGGAATTaa                                                                <  1:1077402/35‑1 (MQ=255)
                                      aaCCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAt                                        <  1:2339465/59‑1 (MQ=255)
                                      aaCCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtc                            >  1:1251717/1‑71 (MQ=255)
                                       aCCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtct                           <  1:1376265/71‑1 (MQ=255)
                                          gCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                                        <  1:219134/55‑1 (MQ=255)
                                           cAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGttt                                                   >  1:1020368/1‑43 (MQ=255)
                                           cAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                               >  1:1011927/1‑63 (MQ=255)
                                            aaTAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGa                      <  1:2260652/71‑1 (MQ=255)
                                             aTAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGa                                         >  1:2590057/1‑51 (MQ=255)
                                             aTAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGa                      >  1:2342284/1‑70 (MQ=255)
                                             aTAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGa                      >  1:1618071/1‑70 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                           <  1:262015/48‑1 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaa                    <  1:983611/71‑1 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaa                    <  1:1761228/71‑1 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGAATTAACGAGACgt                                                        <  1:1007046/35‑1 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAt                                        <  1:843864/51‑1 (MQ=255)
                                              tAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTa                                      <  1:737841/53‑1 (MQ=255)
                                               aaaaTGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                                        >  1:856892/1‑50 (MQ=255)
                                               aaaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                         >  1:857661/1‑49 (MQ=255)
                                               aaaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTTCGCTGGAGTCTCAGGaaaa                   >  1:2253412/1‑71 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                         >  1:261605/1‑48 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                         >  1:2624132/1‑48 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACg                                    >  1:1900778/1‑53 (MQ=255)
                                                aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                     >  1:814545/1‑68 (MQ=255)
                                                 aaTGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                                        <  1:1227064/48‑1 (MQ=255)
                                                   tGCTGAAGGCTCAGGGGATTAACGAGACGTTGttt                                                   <  1:2037732/35‑1 (MQ=255)
                                                     cTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACcaca              <  1:931131/70‑1 (MQ=255)
                                                     cTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAg             >  1:182890/1‑71 (MQ=255)
                                                     cTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAg             >  1:734609/1‑71 (MQ=255)
                                                     cTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAg             <  1:1205582/71‑1 (MQ=255)
                                                     cTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAg             >  1:1281599/1‑71 (MQ=255)
                                                       gAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAt                                        <  1:435483/42‑1 (MQ=255)
                                                       gAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAt                                        <  1:334928/42‑1 (MQ=255)
                                                        aaGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGa                                         >  1:1595071/1‑40 (MQ=255)
                                                        aaGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTAt          <  1:1902688/71‑1 (MQ=255)
                                                        aaGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                         >  1:666124/1‑40 (MQ=255)
                                                        aaGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg                       <  1:118035/58‑1 (MQ=255)
                                                         aGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTAtt         <  1:2258024/71‑1 (MQ=255)
                                                           gCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACg                                    >  1:1106965/1‑42 (MQ=255)
                                                            cTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACcaca              <  1:620867/63‑1 (MQ=255)
                                                            cTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAg      <  1:516416/71‑1 (MQ=255)
                                                            cTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAg      <  1:1398415/71‑1 (MQ=255)
                                                             tCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                     <  1:1085161/55‑1 (MQ=255)
                                                             tCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                     <  1:1063022/55‑1 (MQ=255)
                                                             tCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                         <  1:1674386/35‑1 (MQ=255)
                                                             tCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                         <  1:2638093/35‑1 (MQ=255)
                                                             tCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                         <  1:1088907/35‑1 (MQ=255)
                                                                gggtATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaa                    <  1:1905318/53‑1 (MQ=255)
                                                                ggggATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg                                 <  1:36579/40‑1 (MQ=255)
                                                                gggaATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCtt  <  1:2597364/71‑1 (MQ=255)
                                                                 gggATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                     <  1:736908/51‑1 (MQ=255)
                                                                   |                                                                   
CTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTT  >  W3110S.gb/2809514‑2809648

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: