Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,415,011 A→G L215L (CTT→CTC aslB ← predicted regulator of arylsulfatase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,415,0110AG98.2% 206.9 / NA 56L215L (CTT→CTCaslBpredicted regulator of arylsulfatase activity
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/35);  minor base T (0/1);  total (20/36)
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

CGAACTCCTCGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTAAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATAAACTGGATGAATTC  >  minE/2414944‑2415076
                                                                   |                                                                 
cGAACTCCTCGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGt                                                                >  1:1934181/1‑71 (MQ=255)
  aaCTCCTCGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTcc                                                              <  1:3512572/71‑1 (MQ=255)
   aCTCCTCGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCa                                                             <  1:109030/71‑1 (MQ=255)
   aCTCCTCGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCa                                                             <  1:337113/71‑1 (MQ=255)
        tcGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTcc                                                              <  1:320284/65‑1 (MQ=255)
         cGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAg                                                                 >  1:1175205/1‑61 (MQ=255)
         cGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGAATGTAACTTGAGTCCATCGCgg                                                       >  1:2319976/1‑71 (MQ=255)
             gcgcACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCa                                                             <  1:1976246/61‑1 (MQ=255)
               gcACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAgtt                                                 <  1:2905159/71‑1 (MQ=255)
                cACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAgttg                                                <  1:957465/71‑1 (MQ=255)
                  ccGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAg                                                                 >  1:981466/1‑52 (MQ=255)
                  ccGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAg                                                                 >  1:1195578/1‑52 (MQ=255)
                  ccGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAg                                                                 >  1:123291/1‑52 (MQ=255)
                  ccGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGt                                                                <  1:1312792/53‑1 (MQ=255)
                  ccGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAgttgtt                                              <  1:1555977/71‑1 (MQ=255)
                   cGACCATTCCGTTATCTCACCCTGAATATCACCAGGCGCATGTAACTTGAg                                                                 <  1:1773958/51‑1 (MQ=255)
                   cGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGttt                                             >  1:1184129/1‑71 (MQ=255)
                    gACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGttt                                             >  1:2912212/1‑70 (MQ=255)
                     aCCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCggg                                                      >  1:258618/1‑60 (MQ=255)
                     aCCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTtca                                           <  1:3549946/71‑1 (MQ=255)
                      ccATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTtcat                                          <  1:1805925/71‑1 (MQ=255)
                         ttCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAg                                       >  1:2679982/1‑71 (MQ=255)
                            cGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTtca                                           <  1:3182468/64‑1 (MQ=255)
                            cGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGcc                                     <  1:2150174/70‑1 (MQ=255)
                            cGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGcc                                     <  1:913180/70‑1 (MQ=255)
                            cGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGcc                                     <  1:855948/70‑1 (MQ=255)
                            cGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCa                                    <  1:2444435/71‑1 (MQ=255)
                             gTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCTTCAGCCAg                                   <  1:917803/71‑1 (MQ=255)
                              ttAGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAgg                                  <  1:2890889/71‑1 (MQ=255)
                                aGCTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGgc                                 <  1:3029048/70‑1 (MQ=255)
                                  cTCACCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAgttgt                                               <  1:1578039/54‑1 (MQ=255)
                                     aCCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTc                             <  1:917672/69‑1 (MQ=255)
                                     aCCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTc                             <  1:2105408/69‑1 (MQ=255)
                                     aCCCTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTc                             <  1:4140029/69‑1 (MQ=255)
                                      cccTGAATATCACCAGGCGCATGTAACTTGAGTcc                                                              <  1:4213/35‑1 (MQ=255)
                                      cccTGAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGac                          >  1:3452223/1‑71 (MQ=255)
                                          gAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGcc                                     >  1:132425/1‑56 (MQ=255)
                                          gAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGaccacc                      >  1:1688366/1‑71 (MQ=255)
                                          gAATATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGa                           >  1:1457153/1‑66 (MQ=255)
                                           aaTATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGaccac                       <  1:1304318/69‑1 (MQ=255)
                                           aaTATCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGCCCACCg                     >  1:434445/1‑71 (MQ=255)
                                            atatCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTTTTTCATCAGCCAGGCGCTCGACCACCgg                    <  1:2591342/71‑1 (MQ=255)
                                               tCACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGaa                  >  1:2252280/1‑70 (MQ=255)
                                                cACCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAAt                 >  1:3314602/1‑70 (MQ=255)
                                                 aCCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATaa               <  1:3054316/71‑1 (MQ=255)
                                                 aCCAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATa                >  1:4127104/1‑70 (MQ=255)
                                                  ccAGGCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGac                          >  1:3843034/1‑59 (MQ=255)
                                                     ggCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTtcat                                          >  1:1540493/1‑40 (MQ=255)
                                                     ggCGCATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATAAACTg           <  1:1044763/71‑1 (MQ=255)
                                                      gcgcATGTAACTTTAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGgcg                                <  1:2129018/49‑1 (MQ=255)
                                                      gcgcATGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATAAACTgg          <  1:3471846/71‑1 (MQ=255)
                                                          aTGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATAAACTGGATGa      <  1:833568/71‑1 (MQ=255)
                                                           tGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGc                                      <  1:3390977/38‑1 (MQ=255)
                                                           tGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATAAACTGGATGaa     <  1:779468/71‑1 (MQ=255)
                                                           tGTAACTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATAAACTGGATGaa     <  1:3249919/71‑1 (MQ=255)
                                                               aCTTGAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATAAACTGGATGAATTc  <  1:1170250/70‑1 (MQ=255)
                                                                   |                                                                 
CGAACTCCTCGGGGCGCACCGACCATTCCGTTAGCTCACCCTGAATATCACCAGGCGCATGTAACTTAAGTCCATCGCGGGCAGTTGTTTCATCAGCCAGGCGCTCGACCACCGGAATAAACTGGATGAATTC  >  minE/2414944‑2415076

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: