Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,216,701 T→C 19.3% K638K (AAA→AAG yhdA ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,216,7010TC19.3% 126.1 / 22.9 57K638K (AAA→AAGyhdAconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base C (5/6);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGC  >  minE/2216631‑2216765
                                                                      |                                                                
cGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATAtttt                                                                  <  1:589400/71‑1 (MQ=255)
  cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATAtttttt                                                                >  1:2483599/1‑71 (MQ=255)
  cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATAtttctt                                                                <  1:2318608/71‑1 (MQ=255)
   aGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTc                                                               <  1:160478/71‑1 (MQ=255)
      gcgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTCTTCACg                                                            >  1:2937501/1‑71 (MQ=255)
       cgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGt                                                           <  1:182253/71‑1 (MQ=255)
       cgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGt                                                           <  1:14258/71‑1 (MQ=255)
         cATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATAttttt                                                                 >  1:2116820/1‑63 (MQ=255)
           ttattCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGt                                                       <  1:499667/71‑1 (MQ=255)
           ttattCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGt                                                       <  1:1416848/71‑1 (MQ=255)
            tattCTACGTGAAAACGGATTAAACGTCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTa                                                      >  1:3217715/1‑71 (MQ=255)
            tattCTACGTGAAAACGGATTAAACGGCAGGTTACACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTa                                                      >  1:1192243/1‑71 (MQ=255)
             attCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTAt                                                     >  1:1292017/1‑71 (MQ=255)
                 tACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTa                                                         >  1:7359/1‑63 (MQ=255)
                  aCGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGt                                                >  1:955775/1‑71 (MQ=255)
                     tGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCaa                                                  <  1:1873634/66‑1 (MQ=255)
                        aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          >  1:2565097/1‑71 (MQ=255)
                        aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtg                                           <  1:2132173/70‑1 (MQ=255)
                        aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtg                                           <  1:3707893/70‑1 (MQ=255)
                        aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCt                                            >  1:1647853/1‑69 (MQ=255)
                         aaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          >  1:393569/1‑70 (MQ=255)
                         aaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          >  1:1894603/1‑70 (MQ=255)
                         aaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          >  1:861221/1‑70 (MQ=255)
                          aCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATAttttt                                                                 <  1:1512080/46‑1 (MQ=255)
                          aCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATAttttt                                                                 <  1:2916745/46‑1 (MQ=255)
                          aCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGa                                        <  1:768462/71‑1 (MQ=255)
                          aCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGa                                        <  1:1324626/71‑1 (MQ=255)
                            ggATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          >  1:1844358/1‑67 (MQ=255)
                            ggATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAgg                                      <  1:381860/71‑1 (MQ=255)
                            ggATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAg                                       >  1:1112851/1‑70 (MQ=255)
                              aTTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCt                                            <  1:1090505/63‑1 (MQ=255)
                              aTTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCt                                            <  1:948790/63‑1 (MQ=255)
                                tAAACGGCAGGTTAAACCTAGTATCTTTGTGAATATTTTTTCACGtt                                                          <  1:2846912/47‑1 (MQ=255)
                                 aaaCGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATc                                                    >  1:3390180/1‑52 (MQ=255)
                                 aaaCGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGAcgc                                  >  1:436728/1‑70 (MQ=255)
                                    cGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTa                                                      >  1:2356592/1‑47 (MQ=255)
                                         ggTTAAACCGAGTATCTTTGTGAATATTTCTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGaaa                              <  1:505396/66‑1 (MQ=255)
                                         ggTTAAACCGAGTATCTTTGTGAATATTTCTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGaaa                              <  1:426917/66‑1 (MQ=255)
                                         ggTTAAACCGAGTATCTTTGTGAATATTTCTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGaaa                              <  1:2017481/66‑1 (MQ=255)
                                          gTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          <  1:154454/53‑1 (MQ=255)
                                           ttAAACCGAGTATCTTTGTGAATATTTTTTCACGt                                                           <  1:1773221/35‑1 (MQ=255)
                                           ttAAACCGAGTATCTTTGTGAATATTTCTTCACGTTAGTATc                                                    <  1:2347281/42‑1 (MQ=255)
                                           ttAAACCGAGTATCTTTGTGAATATTTCTTCACGTTAGTATc                                                    <  1:3714512/42‑1 (MQ=255)
                                            tAAACGGAGTATGTTTGTGAATATTTCTTCACGTTAGTATc                                                    <  1:2844600/41‑1 (MQ=37)
                                            tAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATccc                       >  1:220161/1‑70 (MQ=255)
                                            tAAACCGAGTATCTTTGTGAATATTTCTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcccc                      >  1:695361/1‑71 (MQ=255)
                                               aCCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGa                                        >  1:3274316/1‑50 (MQ=255)
                                               aCCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGa                                        >  1:3160564/1‑50 (MQ=255)
                                                 cGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATccc                       >  1:2110758/1‑65 (MQ=255)
                                                  gAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGAcgc                                  >  1:2775059/1‑53 (MQ=255)
                                                  gAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGaaa                              <  1:3730738/57‑1 (MQ=255)
                                                  gAGTATCTTTGTGAATATATCTTCACGTTAGTATCAAGTg                                               >  1:2643433/1‑40 (MQ=255)
                                                   aGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCg               <  1:283122/71‑1 (MQ=255)
                                                       tCTTTGTGAATATTTTTTCACGTTTGTATCAAGTGGCtgt                                          <  1:3453991/40‑1 (MQ=255)
                                                       tCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          <  1:993776/40‑1 (MQ=255)
                                                       tCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          <  1:583609/40‑1 (MQ=255)
                                                       tCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          <  1:2858593/40‑1 (MQ=255)
                                                       tCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCtgt                                          <  1:2666204/40‑1 (MQ=255)
                                                          ttGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTaa        >  1:1575047/1‑71 (MQ=255)
                                                           tgtgAATATTTCTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGc                  >  1:1030320/1‑60 (MQ=255)
                                                            gtgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGt          <  1:1926791/67‑1 (MQ=255)
                                                              gAATATTTCTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACAcc    >  1:893498/1‑71 (MQ=255)
                                                                atatTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcgc  >  1:83054/1‑71 (MQ=255)
                                                                   ttttttCACGTTAGTATCAAGTGGCTGTGAGGACGCGaa                               <  1:2122631/39‑1 (MQ=255)
                                                                   ttttttCACGTTAGTATCAAGTGGCTGTGAGGACGCGaa                               <  1:1553026/39‑1 (MQ=255)
                                                                   ttttttCACGTTAGTATCAAGTGGCTGTGAGGACGCGaa                               <  1:1187509/39‑1 (MQ=255)
                                                                      |                                                                
CGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGC  >  minE/2216631‑2216765

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: