Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 675,913 A→G 31.0% S280S (TCT→TCC 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,9130AG31.0% 76.1 / 42.7 87S280S (TCT→TCCompAouter membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base G (15/12);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAC  >  minE/675852‑675977
                                                             |                                                                
ccATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAg                                                           <  1:3146594/69‑1 (MQ=255)
ccATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAg                                                           <  1:1179102/69‑1 (MQ=255)
ccATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAg                                                           <  1:851438/69‑1 (MQ=255)
ccATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAg                                                           <  1:2642207/69‑1 (MQ=255)
     ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTTAGCACgg                                                      >  1:2568979/1‑69 (MQ=255)
     ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGa                                                            >  1:357714/1‑63 (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:2430053/1‑69 (MQ=255)
     ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                      >  1:3145640/1‑69 (MQ=255)
     ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                      >  1:1997802/1‑69 (MQ=255)
     ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACgg                                                      >  1:1898845/1‑69 (MQ=255)
     ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACg                                                       >  1:704882/1‑68 (MQ=255)
     ccACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACg                                                       >  1:1880882/1‑68 (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:2909050/66‑2 (MQ=255)
                             gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                 <  1:1672858/66‑2 (MQ=255)
                             gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                 <  1:1551535/66‑2 (MQ=255)
                             gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                 <  1:570502/66‑2 (MQ=255)
                             gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                 <  1:2821477/66‑2 (MQ=255)
                             gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                 <  1:1929386/66‑2 (MQ=255)
                             gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                 <  1:2262365/66‑2 (MQ=255)
                             gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                 <  1:1986338/66‑2 (MQ=255)
                             gggATACCTTTAGAGATCAGGTAATCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGtta                                 <  1:1261276/66‑2 (MQ=255)
                             gggATACCTTTAGAGATAAGGTAATCAACAACGGACTGAGCACGACGCTCAGACAGAGCCTGGtta                                 <  1:541442/66‑2 (MQ=255)
                                 tACCTTTGGAGATCAGGTAATCAAGAACAGACTGAg                                                           <  1:3356527/36‑1 (MQ=255)
                                   ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                         >  1:1288734/1‑68 (MQ=255)
                                   ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                         >  1:2048330/1‑68 (MQ=255)
                                   ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                         >  1:194314/1‑68 (MQ=255)
                                   ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                         >  1:2704605/1‑68 (MQ=255)
                                   ccTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTc                         >  1:2957978/1‑68 (MQ=255)
                                     tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                       <  1:3127624/68‑1 (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:50809/68‑1 (MQ=255)
                                     tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                       <  1:524226/68‑1 (MQ=255)
                                     tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                       <  1:1655509/68‑1 (MQ=255)
                                     tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                       <  1:1321393/68‑1 (MQ=255)
                                     tttGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAg                       <  1:900529/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: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:1060715/51‑1 (MQ=255)
                                              cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                               <  1:2502578/51‑1 (MQ=255)
                                              cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                               <  1:2158529/51‑1 (MQ=255)
                                              cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                               <  1:2062254/51‑1 (MQ=255)
                                              cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                               <  1:1321615/51‑1 (MQ=255)
                                              cAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa                               <  1:1340432/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:1433699/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:1310775/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:1197877/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:1130608/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:2692606/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:221610/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:2234867/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:1023313/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:3249370/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:2289707/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:2981108/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:3009181/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:3027224/1‑66 (MQ=255)
                                                    aTCAACAACGGACTGAGCACGACGCTCGGACAGAGCCTGGTTATAAGCGTCAGAACCGATGcggtcag        >  1:228046/1‑66 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:2481059/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:946280/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:849164/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:736396/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:2466465/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:1666385/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:1822060/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:2421040/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:2295642/1‑69 (MQ=255)
                                                       aacaacAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTa    >  1:462452/1‑69 (MQ=255)
                                                         caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc  <  1:1253775/69‑1 (MQ=255)
                                                         caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc  <  1:2981404/69‑1 (MQ=255)
                                                         caacaGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAc  <  1:2278233/69‑1 (MQ=255)
                                                             |                                                                
CCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTAAGCGTCAGAACCGATGCGGTCGGTGTAAC  >  minE/675852‑675977

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