Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,809,581 G→A 9.5% 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,5810GA 9.5% 124.2 / 3.9 42G52G (GGG→GGAmprADNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base A (2/2);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCT  >  W3110S.gb/2809522‑2809650
                                                           |                                                                     
tgcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACga                                                               <  1:889141/68‑1 (MQ=255)
tgcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGAc                                                            >  1:978156/1‑71 (MQ=255)
    tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                                                 <  1:1849198/62‑1 (MQ=255)
    tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                                        >  1:1508998/1‑71 (MQ=255)
      cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                                      >  1:2483012/1‑71 (MQ=255)
      cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                      <  1:1410262/71‑1 (MQ=255)
            caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTaa                                                                  >  1:1180625/1‑53 (MQ=255)
            caAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGtttt                                                    >  1:1071912/1‑66 (MQ=255)
             aaaGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                                      <  1:116335/64‑1 (MQ=255)
               agcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCg                                               <  1:1855787/69‑1 (MQ=255)
                gcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                             >  1:656892/1‑70 (MQ=255)
                gcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                             <  1:1640339/70‑1 (MQ=255)
                  aagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                                          <  1:2358806/71‑1 (MQ=255)
                     ctgctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTa                                                    >  1:677735/1‑58 (MQ=255)
                     ctgctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                                          <  1:424055/68‑1 (MQ=255)
                     ctgctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTAc                                       <  1:1208681/71‑1 (MQ=255)
                     ctgctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTAc                                       >  1:411911/1‑71 (MQ=255)
                      tgctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGAtt                                                                    <  1:1389831/41‑1 (MQ=255)
                         tgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACga                                                               <  1:1432046/43‑1 (MQ=255)
                           gagaACCTCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg                                  <  1:805857/70‑1 (MQ=255)
                           gagaACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg                                  <  1:905865/70‑1 (MQ=255)
                           gagaACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                                 <  1:443808/71‑1 (MQ=255)
                                ccGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                                        >  1:1717097/1‑43 (MQ=255)
                                ccGCAAAAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                                        >  1:359971/1‑43 (MQ=255)
                                  gCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGc                                     >  1:1965389/1‑60 (MQ=255)
                                  gCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGc                                     >  1:265895/1‑60 (MQ=255)
                                  gCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGc                                     >  1:292908/1‑60 (MQ=255)
                                     aTAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                                 <  1:675067/61‑1 (MQ=255)
                                     aTAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                       <  1:625038/71‑1 (MQ=255)
                                     aTAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                       <  1:621648/71‑1 (MQ=255)
                                      tAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTa                                                    <  1:1192477/41‑1 (MQ=255)
                                      tAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttgtt                                                      <  1:2772854/39‑1 (MQ=255)
                                        aaaTGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                    >  1:1352589/1‑71 (MQ=255)
                                        aaaTGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                    >  1:210438/1‑71 (MQ=255)
                                        aaaTGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                    >  1:846556/1‑71 (MQ=255)
                                        aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATgg                                                 >  1:1579008/1‑42 (MQ=255)
                                        aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATgg                                                 >  1:1612102/1‑42 (MQ=255)
                                            gCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTAt                                                   >  1:26279/1‑36 (MQ=255)
                                                aaGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTa             <  1:1856412/70‑1 (MQ=255)
                                                   gCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                             <  1:2726426/35‑1 (MQ=255)
                                                   gCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                             <  1:93765/35‑1 (MQ=255)
                                                   gCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                             <  1:1157653/35‑1 (MQ=255)
                                                        tgggATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg                                   <  1:2100964/39‑1 (MQ=255)
                                                          ggATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCt  >  1:1006609/1‑71 (MQ=255)
                                                           gATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGa                        >  1:989646/1‑48 (MQ=255)
                                                           |                                                                     
TGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCT  >  W3110S.gb/2809522‑2809650

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