Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 507,425 C→A 32.2% A71D (GCT→GAT)  ybhH → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE507,4250CA32.2% 126.0 / 39.6 62A71D (GCT→GAT) ybhHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base A (14/6);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

ACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGG  >  minE/507356‑507494
                                                                     |                                                                     
aCGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                      <  1:3811456/71‑1 (MQ=255)
   gTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGAt                                                                   <  1:3162020/71‑1 (MQ=255)
   gTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGAt                                                                   <  1:1451332/71‑1 (MQ=255)
      ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGAt                                                                   <  1:684528/68‑1 (MQ=255)
        ggcggcGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                              >  1:2951130/1‑71 (MQ=255)
        ggcggcGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                              >  1:1318991/1‑71 (MQ=255)
        ggcggcGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                              >  1:1654438/1‑71 (MQ=255)
           ggcggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGtgatga                                                                    >  1:1344417/1‑62 (MQ=255)
                tAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                      <  1:522630/55‑1 (MQ=255)
                 aaTCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGtt                                                     >  1:1824407/1‑71 (MQ=255)
                       cTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                      <  1:3436037/48‑1 (MQ=255)
                       cTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                      <  1:1378931/48‑1 (MQ=255)
                       cTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTgcg                                                 >  1:205680/1‑69 (MQ=255)
                         gACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTgc                                                  >  1:396687/1‑66 (MQ=255)
                         gACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTgc                                                  >  1:4046476/1‑66 (MQ=255)
                         gACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTgc                                                  >  1:1307399/1‑66 (MQ=255)
                         gACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTgc                                                  >  1:1197983/1‑66 (MQ=255)
                           ccAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTa                                           <  1:3043223/71‑1 (MQ=255)
                               tAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCg                                       >  1:854041/1‑71 (MQ=255)
                               tAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCg                                       <  1:606333/71‑1 (MQ=255)
                                aaaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATc                                        >  1:4227718/1‑69 (MQ=255)
                                aaaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                      >  1:107680/1‑71 (MQ=255)
                                aaaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                      >  1:2387103/1‑71 (MQ=255)
                                aaaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                      >  1:3117712/1‑71 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTgatgat                                                                   <  1:1054289/41‑1 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTgatgat                                                                   <  1:3735250/41‑1 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTgatgat                                                                   <  1:3553378/41‑1 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTgatgat                                                                   <  1:3384362/41‑1 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTgatgat                                                                   <  1:1255331/41‑1 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGtt                                     >  1:3909295/1‑71 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGtt                                     >  1:3052838/1‑71 (MQ=255)
                                   gTCGCCATTATTAGCCGTTCCAGCGATCCGTGTgatgat                                                                   <  1:3653091/39‑1 (MQ=38)
                                    tCGCCATTATTTGCCGTTCCAGCGATCGGCGTGCTGATGTCGATTATc                                                         <  1:1486169/48‑1 (MQ=255)
                                    tCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                                  <  1:1963732/71‑1 (MQ=255)
                                      gCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                                  >  1:1304063/1‑69 (MQ=255)
                                      gCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                                  >  1:3734907/1‑69 (MQ=255)
                                        cATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                              <  1:3582269/71‑1 (MQ=255)
                                         attattAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCt                                                        <  1:2747593/44‑1 (MQ=255)
                                          ttattaGTCGTTCCATCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                            >  1:2010327/1‑71 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                            >  1:2462343/1‑71 (MQ=255)
                                           tattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGt                                                      >  1:4017429/1‑44 (MQ=255)
                                           tattaGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAAc                           >  1:68854/1‑71 (MQ=255)
                                           tattaGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAAc                           >  1:3931475/1‑71 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACg                          >  1:1297519/1‑71 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGATGATGTCGAt                                                             >  1:685385/1‑36 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGATGATGTCGAt                                                             >  1:2500025/1‑36 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTc                                    >  1:2592257/1‑61 (MQ=255)
                                            attaGCAGTTCCAGCGATCCGCGTGATGATGTCGAt                                                             >  1:247069/1‑36 (MQ=37)
                                                gCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTc                      <  1:2994921/71‑1 (MQ=255)
                                                gCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTc                      <  1:1658756/71‑1 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                 >  1:3762596/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc                 >  1:3223221/1‑70 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa               <  1:3186730/71‑1 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa               <  1:3430225/71‑1 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa               <  1:3425305/71‑1 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa               <  1:374307/71‑1 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCa               <  1:4179364/71‑1 (MQ=255)
                                                         gCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACg             >  1:3269116/1‑71 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                            >  1:1329040/1‑54 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                            >  1:2951345/1‑54 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                            >  1:1013734/1‑54 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                            >  1:1645447/1‑54 (MQ=255)
                                                           gATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                            >  1:1798406/1‑54 (MQ=255)
                                                             tCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTa         <  1:3767065/71‑1 (MQ=255)
                                                             tCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTa         <  1:3196849/71‑1 (MQ=255)
                                                             tCCGCGTGCTGATGTCGATTATCTGTTTGCGCAAGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTa         <  1:1564407/71‑1 (MQ=255)
                                                                 cgTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTg     <  1:3103279/71‑1 (MQ=255)
                                                                 cgTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTg     <  1:1941655/71‑1 (MQ=255)
                                                                 cgTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTg     <  1:3330119/71‑1 (MQ=255)
                                                                     cTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCgg  >  1:3248158/1‑70 (MQ=255)
                                                                     |                                                                     
ACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGG  >  minE/507356‑507494

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: