Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 307,588 C→T 15.0% A45T (GCT→ACT)  ylaC ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE307,5880CT15.0% 180.5 / 20.9 73A45T (GCT→ACT) ylaCpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (36/26);  new base T (4/7);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CAACCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTTTGTTGT  >  minE/307519‑307656
                                                                     |                                                                    
cAACCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGc                                                                      <  1:395759/70‑1 (MQ=255)
    cAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACg                                                                  >  1:2970981/1‑70 (MQ=255)
    cAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACg                                                                  >  1:3931382/1‑70 (MQ=255)
         gCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACAt                                                             <  1:1116513/70‑1 (MQ=255)
           cagacagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTa                                                                <  1:3688745/65‑1 (MQ=255)
           cagacagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATAcc                                                          <  1:1359800/71‑1 (MQ=255)
            agacagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATAcc                                                          >  1:1340771/1‑70 (MQ=255)
            agacagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCg                                                         <  1:384722/71‑1 (MQ=255)
              acagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGAt                                                       <  1:669740/71‑1 (MQ=255)
               cagaGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGAt                                                       <  1:3270070/70‑1 (MQ=255)
                agagCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGc                                                                      >  1:3159999/1‑54 (MQ=255)
                agagCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATaa                                                     <  1:3462431/71‑1 (MQ=255)
                 gagCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATaaa                                                    <  1:943252/71‑1 (MQ=255)
                 gagCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATaaa                                                    <  1:3356538/71‑1 (MQ=255)
                  agCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATaaa                                                    >  1:3812928/1‑70 (MQ=255)
                  agCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATaa                                                     >  1:2840167/1‑69 (MQ=255)
                  agCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAAc                                                   <  1:2194591/71‑1 (MQ=255)
                    ccTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAg                                                 >  1:3579794/1‑71 (MQ=255)
                    ccTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAg                                                 >  1:1284315/1‑71 (MQ=255)
                      tGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACATcc                                               >  1:2868856/1‑71 (MQ=255)
                      tGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACATcc                                               >  1:3623364/1‑71 (MQ=255)
                      tGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACATcc                                               >  1:3408109/1‑71 (MQ=255)
                        aCAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCg                                                         <  1:2921892/59‑1 (MQ=255)
                         cAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCg                                            >  1:1322705/1‑71 (MQ=255)
                         cAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCg                                            >  1:3045978/1‑71 (MQ=255)
                            cGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATaaca                                                   <  1:3624818/61‑3 (MQ=255)
                            cGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAAAGTAc                                                               >  1:1580372/1‑49 (MQ=255)
                               ttCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGtt                                      >  1:1816115/1‑71 (MQ=255)
                                tCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAAc                                                                   >  1:3071636/1‑41 (MQ=255)
                                tCGGACTGCAACATCACCGCCAGGGTGCCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGtt                                      >  1:399518/1‑70 (MQ=255)
                                  ggACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTa                                    <  1:1801467/70‑1 (MQ=255)
                                  ggACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTa                                    <  1:161208/70‑1 (MQ=255)
                                  ggACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTAc                                   <  1:1144863/71‑1 (MQ=255)
                                      tGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATa                                                      <  1:3757522/48‑1 (MQ=255)
                                      tGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGa                                >  1:2948881/1‑70 (MQ=255)
                                      tGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGa                                >  1:1074508/1‑70 (MQ=255)
                                        cAACATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGAt                               <  1:1549585/69‑1 (MQ=255)
                                        cAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACa                                                              >  1:994191/1‑38 (MQ=255)
                                        cAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACa                                                              >  1:533358/1‑38 (MQ=255)
                                        cAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACa                                                              >  1:110464/1‑38 (MQ=255)
                                        cAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTAc                                   >  1:2040987/1‑65 (MQ=255)
                                         aaCATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaa                            <  1:397853/71‑1 (MQ=255)
                                         aaCATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaa                            <  1:1981205/71‑1 (MQ=255)
                                         aaCATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaa                            <  1:2028888/71‑1 (MQ=255)
                                          aCATCACCGCCAGGGTGGCAAAAAAAGCAACGTACa                                                              >  1:3871587/1‑36 (MQ=255)
                                          aCATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaaa                           >  1:2155838/1‑71 (MQ=255)
                                           cATCACCGCCAGGGTGGCAAAAAAAGTAACGTACATACAGATAAACAGCCCCGGATGTTTa                                    >  1:1683429/1‑61 (MQ=255)
                                            aTCACCGCCAGGGTGGCAAAAAAAGCAACGTACATAcc                                                          >  1:1790503/1‑38 (MQ=255)
                                             tCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTg                        <  1:774775/71‑1 (MQ=255)
                                                  gCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGAt                                         <  1:804799/49‑1 (MQ=255)
                                                  gCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGa                       <  1:2080596/67‑1 (MQ=255)
                                                   ccAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTa                  <  1:219259/71‑1 (MQ=255)
                                                     aGGGTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGAt                               >  1:2584152/1‑56 (MQ=255)
                                                      gggTGGCAAAAAAAGCAACGTACATACCGATAAACAg                                                 >  1:586469/1‑37 (MQ=255)
                                                       ggTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaaa                           >  1:1945818/1‑58 (MQ=255)
                                                       ggTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATaaaa                           >  1:678352/1‑58 (MQ=255)
                                                        gTGGCAAAAAAAGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAg               <  1:1657613/69‑1 (MQ=255)
                                                            cAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGt         >  1:3007526/1‑71 (MQ=255)
                                                             aaaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGtt        >  1:3141174/1‑71 (MQ=255)
                                                             aaaaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGtt        <  1:1549581/71‑1 (MQ=255)
                                                               aaaaaGCAACGTACATACCGATAAACAGCCCCGGAt                                         >  1:1327092/1‑36 (MQ=255)
                                                               aaaaaGCAACGTACATACCGATAAACAGCCCCGGAt                                         >  1:467522/1‑36 (MQ=255)
                                                               aaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGtt        >  1:933874/1‑69 (MQ=255)
                                                               aaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttg      <  1:768360/71‑1 (MQ=255)
                                                               aaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttg      >  1:3419272/1‑71 (MQ=255)
                                                               aaaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttg      >  1:510832/1‑71 (MQ=255)
                                                                aaaaGCAACGTACATACCGGTAAACAGCCCCGGATGTTTACGGATaaaa                           >  1:1208044/1‑49 (MQ=255)
                                                                aaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttgt     >  1:3290283/1‑71 (MQ=255)
                                                                aaaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttg      <  1:3581406/70‑1 (MQ=255)
                                                                 aaaGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttgtt    <  1:2046725/71‑1 (MQ=255)
                                                                 aaaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttgtt    >  1:1577626/1‑71 (MQ=255)
                                                                  aaGTAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttgttg   <  1:184796/71‑1 (MQ=255)
                                                                  aaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGTATAAAACTGATACTAAAGCg             <  1:871368/61‑1 (MQ=255)
                                                                  aaGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTttgttg   <  1:1357389/71‑1 (MQ=255)
                                                                   aGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTTtgttgt  <  1:2938493/71‑1 (MQ=255)
                                                                     cAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGt         <  1:2537558/62‑1 (MQ=255)
                                                                     |                                                                    
CAACCAGTAGCCAGACAGAGCCTGACAGCGTTTCGGACTGCAACATCACCGCCAGGGTGGCAAAAAAAGCAACGTACATACCGATAAACAGCCCCGGATGTTTACGGATAAAACTGATACTAAAGCGGGGTTTGTTGT  >  minE/307519‑307656

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: