Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,139,485 G→A 35.3% I541I (ATC→ATT infB ← fused protein chain initiation factor 2, IF2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,139,4850GA35.3% 55.7 / 44.3 79I541I (ATC→ATTinfBfused protein chain initiation factor 2, IF2
Reads supporting (aligned to +/- strand):  ref base G (19/32);  new base A (14/14);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGCAG  >  minE/2139420‑2139550
                                                                 |                                                                 
ccTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGAt                                                                 >  1:3047362/1‑68 (MQ=255)
ccTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATa                                                                >  1:3305435/1‑69 (MQ=255)
       gTACCGCTTTCAGCTCCAGTACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:707039/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:791962/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:2458758/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:572305/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:40694/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:2866393/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:3011418/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:174762/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:313428/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:1599542/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:380722/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:1312670/68‑1 (MQ=255)
       gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:1049474/68‑1 (MQ=255)
        tACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                                          <  1:1634696/67‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:1567664/62‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:155941/62‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:3083444/62‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:3248910/62‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:194449/62‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:2568221/62‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:120192/62‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:3007097/62‑1 (MQ=255)
           cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:2099723/62‑1 (MQ=255)
            gCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:670976/61‑1 (MQ=255)
            gCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCt                                                    >  1:1471742/1‑69 (MQ=255)
            gCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCt                                                    >  1:1980930/1‑69 (MQ=255)
               ttCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCttt                                                  >  1:856666/1‑68 (MQ=255)
               ttCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCttt                                                  >  1:706850/1‑68 (MQ=255)
                  aGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCg                                                <  1:1568392/67‑1 (MQ=255)
                  aGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCg                                                <  1:950598/67‑1 (MQ=255)
                             ttCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                                            <  1:637686/44‑1 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTgg                            >  1:580478/1‑67 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGc                           >  1:3318041/1‑68 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:2383937/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:1529497/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:1536157/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:3179133/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:1740267/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:1822065/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:792330/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:2887950/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:571249/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:2196697/1‑69 (MQ=255)
                                      cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt                          >  1:2348868/1‑69 (MQ=255)
                                         cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       <  1:942385/69‑1 (MQ=255)
                                         cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       <  1:970476/69‑1 (MQ=255)
                                         cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       <  1:3248211/69‑1 (MQ=255)
                                         cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       <  1:873639/69‑1 (MQ=255)
                                         cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       <  1:1216469/69‑1 (MQ=255)
                                                   tCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCTCCCCAc             <  1:7927/69‑1 (MQ=255)
                                                   tCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCAc             <  1:621142/69‑1 (MQ=255)
                                                   tCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCAc             <  1:105333/69‑1 (MQ=255)
                                                   tCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCAc             <  1:2864560/69‑1 (MQ=255)
                                                   tCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCAc             <  1:79391/69‑1 (MQ=255)
                                                   tCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCAc             <  1:2128405/69‑1 (MQ=255)
                                                   tCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCAc             <  1:962867/69‑1 (MQ=255)
                                                   tCCAGCAGTTCATCGATACAGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCAc             <  1:752179/69‑1 (MQ=255)
                                                     cagcagTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:369260/1‑57 (MQ=255)
                                                     cagcagTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:2570785/1‑57 (MQ=255)
                                                     cagcagTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:2081503/1‑57 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:702544/1‑49 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:978482/1‑49 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:376783/1‑49 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:3388946/1‑49 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:3190124/1‑49 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:2946090/1‑49 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:2475791/1‑49 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:198754/1‑49 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGc    >  1:27755/1‑68 (MQ=255)
                                                             cATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGc    >  1:180734/1‑68 (MQ=255)
                                                             cATCAATACCGGTACCAGCTTTCGCAGATACGTGTACGAACTGGCTTTc                       >  1:1632884/1‑49 (MQ=39)
                                                               tCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGCAg  <  1:580714/68‑1 (MQ=255)
                                                               tCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGCAg  <  1:2907632/68‑1 (MQ=255)
                                                               tCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGCAg  <  1:2578427/68‑1 (MQ=255)
                                                               tCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGCAg  <  1:1985149/68‑1 (MQ=255)
                                                               tCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGCAg  <  1:1222482/68‑1 (MQ=255)
                                                               tCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCATCCGGCAg  <  1:3111268/68‑1 (MQ=255)
                                                                 |                                                                 
CCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTCACCGCCCCACTCTTCCGGCAG  >  minE/2139420‑2139550

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: