Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,139,470 G→A 58.1% D546D (GAC→GAT 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,4700GA58.1% 21.5 / 54.5 74D546D (GAC→GATinfBfused protein chain initiation factor 2, IF2
Reads supporting (aligned to +/- strand):  ref base G (6/25);  new base A (13/30);  total (19/55)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTC  >  minE/2139416‑2139529
                                                      |                                                           
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATc                                               <  1:2405090/69‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATc                                               <  1:642440/69‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATc                                               <  1:1500168/69‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATc                                               <  1:3086056/69‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATc                                               <  1:3027149/69‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:545384/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:24358/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:2340034/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:2238373/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:2146302/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:2025024/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:2949073/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:1892968/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:1819630/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:1729242/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:3083853/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:492093/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:3418183/68‑1 (MQ=255)
cATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGAATCCAGCAGTTCAt                                                <  1:615467/68‑1 (MQ=255)
cATACCTTTACGTACAGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATc                                               <  1:512725/69‑1 (MQ=255)
    ccTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGAt                                            >  1:3047362/1‑68 (MQ=255)
    ccTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATa                                           >  1:3305435/1‑69 (MQ=255)
          cGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:2876640/58‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGTACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:707039/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:1049474/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:2866393/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:313428/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:2458758/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:380722/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:40694/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:3011418/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:572305/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:1312670/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:1599542/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:791962/68‑1 (MQ=255)
           gTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:174762/68‑1 (MQ=255)
            tACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCATCAATACCGGTa                                     <  1:1634696/67‑1 (MQ=255)
             aCCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATc                                               <  1:1934783/56‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:1567664/62‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:120192/62‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:3083444/62‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:3248910/62‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:3007097/62‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:194449/62‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:155941/62‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:2099723/62‑1 (MQ=255)
               cGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCgg                                       <  1:2568221/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)
                  tttCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCATCCAGCAGTTCAt                                                <  1:180394/50‑1 (MQ=255)
                   ttCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCttt                             >  1:706850/1‑68 (MQ=255)
                   ttCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCttt                             >  1:856666/1‑68 (MQ=255)
                      aGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCg                           <  1:950598/67‑1 (MQ=255)
                      aGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCg                           <  1:1568392/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:792330/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:2887950/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:571249/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:3179133/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:2383937/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:2348868/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:2196697/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:1822065/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:1740267/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:1536157/1‑69 (MQ=255)
                                          cagcagGATAGCATCCAGCAGTTCATCAATACCGGTACCCGCTTTTGCAGATACGTGTACGAACTGGCt     >  1:1529497/1‑69 (MQ=255)
                                             cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc  <  1:3248211/69‑1 (MQ=255)
                                             cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc  <  1:873639/69‑1 (MQ=255)
                                             cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc  <  1:942385/69‑1 (MQ=255)
                                             cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc  <  1:1216469/69‑1 (MQ=255)
                                             cagGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTc  <  1:970476/69‑1 (MQ=255)
                                                      |                                                           
CATACCTTTACGTACCGCTTTCAGCTCCAGAACTTCCGCCTGCAGCAGGATAGCGTCCAGCAGTTCATCGATACCGGTACCCGCTTTCGCAGATACGTGTACGAACTGGCTTTC  >  minE/2139416‑2139529

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: