Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,695,852 G→A 38.9% A53V (GCC→GTC)  yfhQ ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,695,8520GA38.9% 43.0 / 56.2 54A53V (GCC→GTC) yfhQpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (21/12);  new base A (11/10);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGGCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAGCCACAGAT  >  minE/1695783‑1695915
                                                                     |                                                               
accaGGCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGc                                                                >  1:283306/1‑71 (MQ=255)
accaGGCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGc                                                                >  1:390615/1‑71 (MQ=255)
accaGGCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGc                                                                >  1:1915911/1‑71 (MQ=255)
 ccaGGCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGAcc                                                               <  1:2575458/71‑1 (MQ=255)
  caGGCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGAccc                                                              >  1:1026146/1‑71 (MQ=255)
  caGGCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGAccc                                                              >  1:2849674/1‑71 (MQ=255)
    ggCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCTATCACATCGCTGACCCct                                                            <  1:3757315/71‑1 (MQ=255)
    ggCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGcccc                                                             >  1:2501840/1‑70 (MQ=255)
     gCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCctg                                                           >  1:2576381/1‑71 (MQ=255)
     gCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCct                                                            <  1:701253/70‑1 (MQ=255)
     gCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGACCCct                                                            <  1:2711734/70‑1 (MQ=255)
         cAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCtgctgc                                                       >  1:2084675/1‑71 (MQ=255)
           aCCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGACCCctgct                                                         >  1:1610681/1‑67 (MQ=255)
               gCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGACCCctg                                                           <  1:3115453/61‑1 (MQ=255)
                      cTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAg                                                    >  1:1989147/1‑61 (MQ=255)
                      cTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGcc                                          >  1:3977198/1‑71 (MQ=255)
                      cTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGc                                           >  1:2200837/1‑70 (MQ=255)
                       ttCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCa                                                     <  1:3355094/59‑1 (MQ=255)
                         cgtcgAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAg                                                    >  1:1586953/1‑58 (MQ=255)
                          gtcgAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGACCCCtgctgc                                                       <  1:1764456/54‑1 (MQ=255)
                             gAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGt                                   >  1:1687568/1‑71 (MQ=255)
                              aGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGTc                                  <  1:2948568/71‑1 (MQ=255)
                               gCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCg                                 <  1:827858/71‑1 (MQ=255)
                               gCGTATCGACGATGTGAGCATTACCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGTCg                                 <  1:2181021/71‑1 (MQ=255)
                                  tATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAAt                                              >  1:1552463/1‑55 (MQ=255)
                                    tCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGttt                            <  1:3804908/71‑1 (MQ=255)
                                      gacgaTGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCa                          >  1:1769848/1‑71 (MQ=255)
                                        cgaTGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCAcc                        <  1:3410765/71‑1 (MQ=255)
                                        cgaTGTGAGCATTACCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCAcc                        <  1:870772/71‑1 (MQ=255)
                                         gaTGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCa                       <  1:2478714/71‑1 (MQ=255)
                                         gaTGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCa                       <  1:3215728/71‑1 (MQ=255)
                                           tgtgAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGt                     >  1:216531/1‑71 (MQ=255)
                                            gtgAGCATTACCAATCACATCGCTGACCCCTGCTGCCAGGGCaa                                               >  1:1637807/1‑44 (MQ=255)
                                            gtgAGCATTACCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTg                    >  1:3448824/1‑71 (MQ=255)
                                                gCATTACCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGAtt                <  1:649398/71‑1 (MQ=255)
                                                 cATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTc                                  <  1:1916477/52‑1 (MQ=255)
                                                 cATTACCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGAtt                >  1:3553896/1‑70 (MQ=255)
                                                   ttACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCAcc                        <  1:600597/60‑1 (MQ=255)
                                                   ttACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAAc             <  1:2716131/71‑1 (MQ=255)
                                                   ttACCAATCACATCGCTGACCCCTGCTGCCAGGGCaa                                               <  1:2170921/37‑1 (MQ=255)
                                                     aCCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAg                      >  1:2602539/1‑60 (MQ=255)
                                                     aCCAATCACATCGCTGACCCCTGCTGCCAGGGCAAt                                              >  1:3766156/1‑36 (MQ=255)
                                                     aCCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAAc             >  1:607069/1‑69 (MQ=255)
                                                     aCCAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAAc             >  1:2248869/1‑69 (MQ=255)
                                                      ccAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAg          >  1:3703034/1‑71 (MQ=255)
                                                      ccAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAg          >  1:3995261/1‑71 (MQ=255)
                                                       cAATCACATCGCTGACCCCTGCTGCCAGGGCAATCGc                                           >  1:117312/1‑37 (MQ=255)
                                                          tCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAGCca       >  1:751132/1‑70 (MQ=255)
                                                           cacaTCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCAc                         >  1:4137597/1‑51 (MQ=255)
                                                            acaTCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCaa                         >  1:324953/1‑49 (MQ=255)
                                                            acaTCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAGCca       >  1:3747210/1‑68 (MQ=255)
                                                            acaTCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAGCCACAg    <  1:243946/71‑1 (MQ=255)
                                                             caTCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAGCCACAGa   <  1:2228465/71‑1 (MQ=255)
                                                              aTCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAGCCACAGAt  >  1:248360/1‑71 (MQ=255)
                                                                     |                                                               
ACCAGGCTACAACCGGCTAACGCTTCGTCGAGCGTATCGACGATGTGAGCATTACCAATCACATCGCTGGCCCCTGCTGCCAGGGCAATCGCCTGGGAGTCGGGTTTCACCAGTGGATTAACCAGCCACAGAT  >  minE/1695783‑1695915

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: