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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE507,4250CA16.3% 188.8 / 14.2 67A71D (GCT→GAT) ybhHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/20);  new base A (3/8);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 4.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

ACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTG  >  minE/507356‑507491
                                                                     |                                                                  
aCGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                   <  1:594563/71‑1 (MQ=255)
   gTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGa                                                                 >  1:1444749/1‑70 (MQ=255)
      ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTc                                                             >  1:1484277/1‑71 (MQ=255)
      ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTc                                                             >  1:1094162/1‑71 (MQ=255)
        ggcggcGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                           >  1:1643407/1‑71 (MQ=255)
        ggcggcGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGa                                                           >  1:668829/1‑71 (MQ=255)
              ggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCt                                                     <  1:530867/71‑1 (MQ=255)
               gTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTg                                                    <  1:904493/71‑1 (MQ=255)
                tAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGt                                                   <  1:479779/71‑1 (MQ=255)
                tAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGt                                                   <  1:600602/71‑1 (MQ=255)
                tAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGt                                                   <  1:1187740/71‑1 (MQ=255)
                tAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGt                                                   <  1:2025606/71‑1 (MQ=255)
                tAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGt                                                   <  1:2146481/71‑1 (MQ=255)
                  aTCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGa                                                                 >  1:1681516/1‑55 (MQ=255)
                  aTCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGa                                                                 >  1:22913/1‑55 (MQ=255)
                   tCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTg                                                <  1:684670/71‑1 (MQ=255)
                   tCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTg                                                <  1:529354/71‑1 (MQ=255)
                     cACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGtt                                                  >  1:1015053/1‑67 (MQ=255)
                     cACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTgcg                                              >  1:12881/1‑71 (MQ=255)
                         gACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAg                                           >  1:403515/1‑70 (MQ=255)
                           ccAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCt                                                     >  1:1866028/1‑58 (MQ=255)
                            cAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCt                                                                   <  1:213584/43‑1 (MQ=255)
                                aaaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                   >  1:1419817/1‑71 (MQ=255)
                                aaaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGt                                   >  1:309958/1‑71 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATc                                                      >  1:420205/1‑51 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATc                                                      >  1:1643796/1‑51 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATc                                                      >  1:1936139/1‑51 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATc                                                      >  1:771432/1‑51 (MQ=255)
                                 aaGTCGCCATTATTAGCCGTTACAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCg                                    >  1:63483/1‑69 (MQ=255)
                                    tCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                               <  1:1628770/71‑1 (MQ=255)
                                     cGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                               >  1:1362622/1‑70 (MQ=255)
                                      gCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAgg                                          >  1:198633/1‑58 (MQ=255)
                                      gCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                               >  1:520425/1‑69 (MQ=255)
                                      gCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCAt                               >  1:848173/1‑69 (MQ=255)
                                      gCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                             <  1:39850/71‑1 (MQ=255)
                                         attattAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAg                            >  1:313720/1‑69 (MQ=255)
                                         attattAGCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGc                           >  1:301351/1‑70 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                         <  1:2022894/71‑1 (MQ=255)
                                          ttattaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCaa                         >  1:477524/1‑71 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGttt                                                 >  1:1223329/1‑45 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGttt                                                 >  1:1243218/1‑45 (MQ=255)
                                            attaGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAAc                        <  1:1616689/70‑1 (MQ=255)
                                               aGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                             >  1:1250102/1‑62 (MQ=255)
                                                gCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTc                   <  1:1084232/71‑1 (MQ=255)
                                                gCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTc                   <  1:1031246/71‑1 (MQ=255)
                                                gCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACgtgt                    <  1:469591/70‑1 (MQ=255)
                                                gCCGTTCCAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACgtgt                    <  1:997419/70‑1 (MQ=255)
                                                 ccGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAg                            <  1:362347/61‑1 (MQ=255)
                                                  cGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGa                 <  1:1171179/71‑1 (MQ=255)
                                                   gTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGAt                <  1:2118205/71‑1 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGa                             >  1:1012831/1‑55 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATa               >  1:380399/1‑69 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc              >  1:1425501/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc              >  1:813307/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc              >  1:395138/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc              >  1:1700392/1‑70 (MQ=255)
                                                      ccAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc              <  1:388201/70‑1 (MQ=255)
                                                      ccAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc              <  1:1331863/70‑1 (MQ=255)
                                                      ccAGCGATCCGCGTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATAc              <  1:1198880/70‑1 (MQ=255)
                                                       cAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAAc                        <  1:755013/59‑1 (MQ=255)
                                                        aGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAg                            >  1:2110944/1‑54 (MQ=255)
                                                         gCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACg          >  1:1997233/1‑71 (MQ=255)
                                                         gCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACg          <  1:480219/71‑1 (MQ=255)
                                                         gCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACg          >  1:1115571/1‑71 (MQ=255)
                                                         gCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACg          >  1:1106843/1‑71 (MQ=255)
                                                             tCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTa      <  1:680770/71‑1 (MQ=255)
                                                                 cgTGATGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTg  >  1:1720346/1‑71 (MQ=255)
                                                                     |                                                                  
ACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTG  >  minE/507356‑507491

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: