Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 675,925 G→A 32.4% R276R (CGC→CGT ompA ← outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE675,9250GA32.4% 71.8 / 40.8 83R276R (CGC→CGTompAouter membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base G (22/34);  new base A (15/12);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGA  >  minE/675857‑675993
                                                                    |                                                                    
ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTTAGCACgg                                                                      >  1:2568979/1‑69 (MQ=255)
ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                                      >  1:904126/1‑69 (MQ=255)
ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                                      >  1:2736552/1‑69 (MQ=255)
ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                                      >  1:1228900/1‑69 (MQ=255)
ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                                      >  1:1898845/1‑69 (MQ=255)
ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                                      >  1:2430053/1‑69 (MQ=255)
ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                                      >  1:1997802/1‑69 (MQ=255)
ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                                      >  1:3145640/1‑69 (MQ=255)
         gagaTCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCACCAACAGACTGAGCACGGCGCTCg                                                                >  1:3159652/1‑66 (MQ=255)
            aTCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGa                                                          <  1:1154544/69‑1 (MQ=255)
            aTCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGa                                                          <  1:728629/69‑1 (MQ=255)
             tCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGAc                                                         <  1:543566/69‑1 (MQ=255)
              cTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGAc                                                         <  1:505182/68‑1 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:2262365/66‑2 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:2909050/66‑2 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:1986338/66‑2 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:2821477/66‑2 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:1929386/66‑2 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:1672858/66‑2 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:570502/66‑2 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:1261276/66‑2 (MQ=255)
                        gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:1551535/66‑2 (MQ=255)
                        gggATACCTTTAGAGATAAGGTAATCAACAACGGACTGAGCACGACGCTCAGACAGAGCCTGGtta                                                 <  1:541442/66‑2 (MQ=255)
                              ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                         >  1:2704605/1‑68 (MQ=255)
                              ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                         >  1:194314/1‑68 (MQ=255)
                              ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                         >  1:1288734/1‑68 (MQ=255)
                              ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                         >  1:2957978/1‑68 (MQ=255)
                              ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                         >  1:2048330/1‑68 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:983244/68‑1 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:900529/68‑1 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:1321393/68‑1 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:3127624/68‑1 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:524226/68‑1 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:3209340/68‑1 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:395887/68‑1 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:50809/68‑1 (MQ=255)
                                tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:1655509/68‑1 (MQ=255)
                                 ttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:2898122/67‑1 (MQ=255)
                                   ggAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                       <  1:1910114/65‑1 (MQ=255)
                                      gaTCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGAt                                <  1:1929403/69‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:3123330/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:2962861/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:2900348/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:263568/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:501626/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:1060715/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:1321615/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:1340432/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:2062254/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:2502578/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                               <  1:2158529/51‑1 (MQ=255)
                                         cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGTTGTa                                               <  1:2442357/50‑1 (MQ=255)
                                            gTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:3178384/46‑2 (MQ=255)
                                             tAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                 <  1:2754855/45‑2 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTg                                                     >  1:354932/1‑39 (MQ=37)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:1023313/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:3249370/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:1433699/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:2692606/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:1310775/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:221610/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:2234867/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:3027224/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:3009181/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:2981108/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:1197877/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:228046/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:2289707/1‑66 (MQ=255)
                                               aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                        >  1:1130608/1‑66 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:736396/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:849164/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:946280/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:1666385/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:1822060/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:462452/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:2421040/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:2481059/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:2466465/1‑69 (MQ=255)
                                                  aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                    >  1:2295642/1‑69 (MQ=255)
                                                    caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc                  <  1:1253775/69‑1 (MQ=255)
                                                    caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc                  <  1:2278233/69‑1 (MQ=255)
                                                    caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc                  <  1:2981404/69‑1 (MQ=255)
                                                           cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCg                                  >  1:2968604/1‑46 (MQ=255)
                                                           cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAcc                 >  1:3310487/1‑63 (MQ=255)
                                                           cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaac           >  1:3353048/1‑69 (MQ=255)
                                                           cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaac           >  1:2350235/1‑69 (MQ=255)
                                                           cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaac           >  1:2643324/1‑69 (MQ=255)
                                                           cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaac           >  1:2036517/1‑69 (MQ=255)
                                                                cACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaacaac        <  1:1256011/67‑1 (MQ=255)
                                                                cACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaacaac        <  1:917664/67‑1 (MQ=255)
                                                                 aCGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaacaac        <  1:2888098/66‑1 (MQ=255)
                                                                    gcgcTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGa  >  1:3163070/1‑69 (MQ=255)
                                                                    gcgcTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGa  >  1:3304523/1‑69 (MQ=255)
                                                                    |                                                                    
CCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGA  >  minE/675857‑675993

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: