Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 201,711 C→T 32.0% V18V (GTC→GTT gmhB → D,D‑heptose 1,7‑bisphosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE201,7110CT32.0% 76.1 / 40.4 84V18V (GTC→GTTgmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (21/6);  total (50/34)
Fisher's exact test for biased strand distribution p-value = 3.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGC  >  minE/201645‑201775
                                                                  |                                                                
ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCg                                                                 <  1:409917/68‑1 (MQ=255)
ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCg                                                                 <  1:3380534/68‑1 (MQ=255)
ggAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCg                                                                 <  1:2096753/68‑1 (MQ=255)
 gAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCg                                                                 <  1:758572/67‑1 (MQ=255)
 gAGCTATAAAAGGTGGCGAAGAGAGTACCCGCAATTTTTCTTGAGCGTGATGGCACCATTAATGTCg                                                                 <  1:193742/67‑1 (MQ=255)
      ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACg                                                           >  1:1747947/1‑68 (MQ=255)
      ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACg                                                           >  1:1243783/1‑68 (MQ=255)
      ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACg                                                           >  1:1593486/1‑68 (MQ=255)
      ataAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACg                                                           >  1:1288178/1‑68 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:1731457/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:2735871/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:1768720/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:262895/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:2853224/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:3164361/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:646348/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:862791/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:915095/67‑1 (MQ=255)
         aaaGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:967716/67‑1 (MQ=255)
            ggTTGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:106372/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTAt                                                      >  1:1469153/1‑67 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:303313/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:936707/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:1147991/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:864294/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:805308/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:726222/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:456696/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:3002658/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:2920028/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:2841483/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:1018849/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:2733261/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:2324645/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:219567/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:2003379/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:1800797/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:1996438/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATCACGGCTATg                                                     >  1:1744421/1‑68 (MQ=255)
            ggTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTTGATAACGGCTAt                                                      >  1:1218725/1‑67 (MQ=255)
                gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTcc                                                  >  1:1287318/1‑67 (MQ=255)
                gCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTcc                                                  >  1:3430468/1‑67 (MQ=255)
                   aaGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:3359509/57‑1 (MQ=255)
                         gTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGAc                                       <  1:2001256/69‑1 (MQ=255)
                            cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:525369/48‑1 (MQ=255)
                            cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:1266939/48‑1 (MQ=255)
                            cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:1981056/48‑1 (MQ=255)
                            cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:2665577/48‑1 (MQ=255)
                            cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGc                                                         <  1:2786565/48‑1 (MQ=255)
                            cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAAc                                    >  1:2073087/1‑69 (MQ=255)
                            cccGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAAc                                    >  1:1434083/1‑69 (MQ=255)
                              cGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATc                                                              >  1:2335693/1‑41 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:380135/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:3333710/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:169585/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:60441/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:314381/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:1231992/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:3047224/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:1976424/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAttt                           >  1:93757/1‑69 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAtt                            >  1:1646654/1‑68 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAAtt                            >  1:1414502/1‑68 (MQ=255)
                                     tttCTTGACCGTGATGGCACCATTAATGTCGATAACGGCTATGTCCATGAGATCGACAACTTTGAAtt                            >  1:837669/1‑68 (MQ=255)
                                         ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATc                        <  1:2929374/68‑1 (MQ=255)
                                         ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATc                        <  1:1314247/68‑1 (MQ=255)
                                         ttGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATc                        <  1:170017/68‑1 (MQ=255)
                                                tGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGgtgt                <  1:2103242/69‑1 (MQ=255)
                                                tGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGgtgt                <  1:4543/69‑1 (MQ=255)
                                                 gATGGCACCATTAATGTTGATCACGGCTATGTCCATTAGATCGACAACTTTGAATTTATCGACGGTGtt               <  1:1755677/69‑1 (MQ=255)
                                                 gATGGCACCATTAATGTTGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGtt               <  1:1399799/69‑1 (MQ=255)
                                                 gATGGCACCATTAATGTTGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGtt               <  1:3157754/69‑1 (MQ=255)
                                                 gATGGCACCATTAATGTTGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGtt               <  1:3121822/69‑1 (MQ=255)
                                                 gATGGCACCATTAATGTTGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGtt               <  1:845096/69‑1 (MQ=255)
                                                 gATGGCACCATTAATGTTGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGtt               <  1:3056942/69‑1 (MQ=255)
                                                 gATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGtt               <  1:3275772/69‑1 (MQ=255)
                                                  aTGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACgg                   >  1:1451538/1‑64 (MQ=255)
                                                  aTGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACgg                   >  1:763855/1‑64 (MQ=255)
                                                  aTGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACgg                   >  1:592056/1‑64 (MQ=255)
                                                  aTGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACgg                   >  1:1715807/1‑64 (MQ=255)
                                                  aTGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACgg                   >  1:2968946/1‑64 (MQ=255)
                                                  aTGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACg                    >  1:1884847/1‑63 (MQ=255)
                                                             aaTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATg   >  1:3313921/1‑69 (MQ=255)
                                                              aTGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATgc  >  1:298345/1‑69 (MQ=255)
                                                                  |                                                                
GGAGCTATAAAAGGTGGCGAAGAGCGTACCCGCAATTTTTCTTGACCGTGATGGCACCATTAATGTCGATCACGGCTATGTCCATGAGATCGACAACTTTGAATTTATCGACGGTGTTATTGACGCCATGC  >  minE/201645‑201775

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: