Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 675,938 C→G 30.6% G272A (GGT→GCT)  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,9380CG30.6% 105.1 / 59.1 88G272A (GGT→GCT) ompAouter membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base C (25/36);  new base G (15/12);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGAACCGTCTTT  >  minE/675870‑676002
                                                                    |                                                                
tCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGAc                                                                  <  1:543566/69‑1 (MQ=255)
 cTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGAc                                                                  <  1:505182/68‑1 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:2909050/66‑2 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:1929386/66‑2 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:1672858/66‑2 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:1986338/66‑2 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:1551535/66‑2 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:2821477/66‑2 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:1261276/66‑2 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:570502/66‑2 (MQ=255)
           gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                                          <  1:2262365/66‑2 (MQ=255)
           gggATACCTTTAGAGATAAGGTAATCAACAACGGACTGAGCACGACGCTCAGACAGAGCCTGGtta                                                          <  1:541442/66‑2 (MQ=255)
                 ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                                  >  1:1288734/1‑68 (MQ=255)
                 ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                                  >  1:2048330/1‑68 (MQ=255)
                 ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                                  >  1:2957978/1‑68 (MQ=255)
                 ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                                  >  1:194314/1‑68 (MQ=255)
                 ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                                                  >  1:2704605/1‑68 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:983244/68‑1 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:3209340/68‑1 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:395887/68‑1 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:1655509/68‑1 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:3127624/68‑1 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:50809/68‑1 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:1321393/68‑1 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:900529/68‑1 (MQ=255)
                   tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                                                <  1:524226/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:501626/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:2502578/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:263568/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:2900348/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:2962861/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:3123330/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:1060715/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:2158529/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:2062254/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:1340432/51‑1 (MQ=255)
                            cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                                                        <  1:1321615/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:1310775/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:1197877/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:1433699/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:3249370/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:228046/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:3009181/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:2981108/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:1130608/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:3027224/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:221610/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:1023313/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:2234867/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:2289707/1‑66 (MQ=255)
                                  aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag                                 >  1:2692606/1‑66 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:2421040/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:1666385/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:462452/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:2295642/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:736396/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:1822060/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:849164/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:2481059/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:2466465/1‑69 (MQ=255)
                                     aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             >  1:946280/1‑69 (MQ=255)
                                       caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc                           <  1:2981404/69‑1 (MQ=255)
                                       caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc                           <  1:1253775/69‑1 (MQ=255)
                                       caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc                           <  1:2278233/69‑1 (MQ=255)
                                              cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCg                                           >  1:2968604/1‑46 (MQ=255)
                                              cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAcc                          >  1:3310487/1‑63 (MQ=255)
                                              cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaac                    >  1:2350235/1‑69 (MQ=255)
                                              cTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGaac                    >  1:3353048/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:3304523/1‑69 (MQ=255)
                                                       gcgcTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGa           >  1:3163070/1‑69 (MQ=255)
                                                               aCAGACCCTGTTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             <  1:3134190/43‑1 (MQ=255)
                                                               aCAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             <  1:1680294/43‑1 (MQ=255)
                                                               aCAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             <  1:1753768/43‑1 (MQ=255)
                                                               aCAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             <  1:3091800/43‑1 (MQ=255)
                                                               aCAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             <  1:2253610/43‑1 (MQ=255)
                                                               aCAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa                             <  1:2823816/43‑1 (MQ=255)
                                                                cAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGAACCGTCttt  >  1:1871502/1‑69 (MQ=255)
                                                                cAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGAACCGTCttt  >  1:2287221/1‑69 (MQ=255)
                                                                    |                                                                
TCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAACCCAGAACAACTACGGAACCGTCTTT  >  minE/675870‑676002

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: