Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 242,863 G→A 6.6% Q454* (CAA→TAA)  sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE242,8630GA 6.6% 296.5 / 7.5 91Q454* (CAA→TAA) sbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (47/38);  new base A (3/3);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

AGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTCCATGCAGC  >  minE/242793‑242933
                                                                      |                                                                      
aGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTg                                                                        <  1:3140184/71‑1 (MQ=255)
 gTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGc                                                                       >  1:224984/1‑71 (MQ=255)
    gctgcGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGt                                                                     <  1:895784/70‑1 (MQ=255)
    gctgcGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGt                                                                     <  1:3035096/70‑1 (MQ=255)
     ctgcGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTg                                                                    <  1:2171390/70‑1 (MQ=255)
        cGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACa                                                                 <  1:3832308/70‑1 (MQ=255)
        cGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACa                                                                 <  1:793247/70‑1 (MQ=255)
        cGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACa                                                                 <  1:3142750/70‑1 (MQ=255)
        cGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACAt                                                                <  1:3085542/71‑1 (MQ=255)
         gTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACAtt                                                               <  1:2729276/71‑1 (MQ=255)
          tcttttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACAtt                                                               <  1:441588/70‑1 (MQ=255)
           cttttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTACGTGa                                                                   <  1:1901237/65‑1 (MQ=255)
            ttttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTc                                                              >  1:2984996/1‑69 (MQ=255)
            ttttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTg                                                            >  1:3568630/1‑71 (MQ=255)
            ttttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTg                                                            >  1:3009743/1‑71 (MQ=255)
             tttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGAc                                                                  >  1:3161957/1‑64 (MQ=255)
             tttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGAc                                                                  >  1:3431698/1‑64 (MQ=255)
              ttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACAtt                                                               >  1:636827/1‑66 (MQ=255)
              ttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCt                                                             >  1:3253692/1‑68 (MQ=255)
                   tataACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGc                                                      <  1:677254/70‑1 (MQ=255)
                   tataACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTACGTGa                                                                   >  1:2286238/1‑57 (MQ=255)
                     taACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACa                                                                 <  1:1511812/57‑1 (MQ=255)
                     taACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCg                                                     <  1:3169097/69‑1 (MQ=255)
                     taACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGAc                                                   <  1:2085506/71‑1 (MQ=255)
                       aCGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGa                                                                   >  1:3540013/1‑53 (MQ=255)
                         gCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGAc                                                   >  1:645123/1‑67 (MQ=255)
                         gCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGt                                               <  1:781530/71‑1 (MQ=255)
                         gCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGt                                               <  1:1577584/71‑1 (MQ=255)
                          cTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTg                                                                        >  1:3068664/1‑45 (MQ=255)
                          cTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTg                                                                        >  1:1345122/1‑45 (MQ=255)
                            ggCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTg                                                            <  1:2938473/55‑1 (MQ=255)
                             gcgcATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAAc                                            <  1:1687072/70‑1 (MQ=255)
                              cgcATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTg                                                            >  1:4042846/1‑53 (MQ=255)
                              cgcATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTg                                                            >  1:2777830/1‑53 (MQ=255)
                                cATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAAc                                            >  1:153548/1‑67 (MQ=255)
                                cATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTgcg                                        >  1:3718571/1‑71 (MQ=255)
                                 aTTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGa                                                                   <  1:2134876/43‑1 (MQ=255)
                                  tttCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGt                                                                     <  1:3527934/40‑1 (MQ=255)
                                  tttCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGt                                                                     <  1:914090/40‑1 (MQ=255)
                                  tttCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACa                                                                 <  1:2905507/44‑1 (MQ=255)
                                    tCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTg                                                            <  1:3659018/47‑1 (MQ=255)
                                    tCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGa                                                    <  1:1084939/55‑1 (MQ=255)
                                    tCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTgcgc                                       >  1:2722128/1‑68 (MQ=255)
                                    tCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTgcgc                                       >  1:3350938/1‑68 (MQ=255)
                                    tCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTgcgc                                       >  1:2945153/1‑68 (MQ=255)
                                     cGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGa                                                          >  1:2335351/1‑48 (MQ=255)
                                         aaGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGAGATTCTGGa                                                          >  1:542564/1‑44 (MQ=255)
                                         aaGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACa                                                                 >  1:376879/1‑37 (MQ=255)
                                         aaGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGa                                                          >  1:1350999/1‑44 (MQ=255)
                                         aaGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGtt                               >  1:2173527/1‑71 (MQ=255)
                                            tGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGa                                                          >  1:2586987/1‑41 (MQ=255)
                                             gcggcgTTACGTTGCGTCTGTTTTTGCGTGACATTCTGGATAGc                                                      >  1:3488705/1‑44 (MQ=255)
                                             gcggcgTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGAt                                                         >  1:679556/1‑41 (MQ=255)
                                             gcggcgTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGAc                                                   >  1:1483328/1‑47 (MQ=255)
                                             gcggcgTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTaa                                             >  1:2355750/1‑53 (MQ=255)
                                               ggcgTTACGTTGCGTCTGTTCTTGCGTGACATTCTg                                                            <  1:1754292/36‑1 (MQ=255)
                                                   ttACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGttttt                            <  1:1383856/64‑1 (MQ=255)
                                                   ttACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTttgt                          <  1:780055/66‑1 (MQ=255)
                                                   ttACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTggg                      >  1:1608224/1‑70 (MQ=255)
                                                     aCGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGAc                                                   >  1:3641982/1‑39 (MQ=255)
                                                     aCGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAAc                                            >  1:1855310/1‑46 (MQ=255)
                                                     aCGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTttgttg                        >  1:2924546/1‑66 (MQ=255)
                                                     aCGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGaa                   >  1:3113624/1‑71 (MQ=255)
                                                     aCGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGaa                   >  1:2515588/1‑71 (MQ=255)
                                                     aCGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGaa                   >  1:1120377/1‑71 (MQ=255)
                                                     aCGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGaa                   >  1:2346815/1‑71 (MQ=255)
                                                     aCGTTGCGTCTGTTCTTACGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGaa                   >  1:855321/1‑71 (MQ=255)
                                                      cGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGaa                   >  1:1527150/1‑70 (MQ=255)
                                                       gTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTTCGCCAGACGTTTTTGTTGGGGAACa                 <  1:1642179/71‑1 (MQ=255)
                                                       gTTGCGTCTGTTCTTACGTGACATTCTGGATAGCGACCTGt                                               <  1:2603512/41‑1 (MQ=255)
                                                        ttGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGt                                               <  1:2897774/40‑1 (MQ=255)
                                                        ttGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGtt                               <  1:1809274/56‑1 (MQ=255)
                                                        ttGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGtt                               <  1:1457609/56‑1 (MQ=255)
                                                        ttGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACaa                >  1:3994234/1‑71 (MQ=255)
                                                        ttGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACaa                <  1:1056354/71‑1 (MQ=255)
                                                          gCGTCTGTTCTTGCGTGACATTCTGTATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAAt               <  1:2886894/70‑1 (MQ=255)
                                                          gCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTttgt                          >  1:2850586/1‑59 (MQ=255)
                                                          gCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAAt               <  1:1653859/70‑1 (MQ=255)
                                                            gTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTg            >  1:2281243/1‑71 (MQ=255)
                                                             tCTGTTCTTACGTGACATTCTGGATACCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGaa                   <  1:1664879/63‑1 (MQ=255)
                                                               tGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTc          >  1:2892480/1‑70 (MQ=255)
                                                               tGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTc          >  1:3525308/1‑70 (MQ=255)
                                                               tGTTCTTACGTGACATTCTGGATAGCGACCTGTACCTGCGCCAGACGTTTTTGTTggg                      >  1:2729346/1‑58 (MQ=38)
                                                                gTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTttgt                          >  1:1011985/1‑53 (MQ=255)
                                                                 ttCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACg                                 <  1:89279/45‑1 (MQ=255)
                                                                 ttCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGttt                              >  1:2751219/1‑48 (MQ=255)
                                                                 ttCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAAc                  <  1:414230/60‑1 (MQ=255)
                                                                 ttCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAAc                  <  1:1621359/60‑1 (MQ=255)
                                                                 ttCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACa                 >  1:4074965/1‑61 (MQ=255)
                                                                 ttCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACa                 <  1:1143621/61‑1 (MQ=255)
                                                                  tCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGtttt                             >  1:1997789/1‑48 (MQ=255)
                                                                  tCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGtttt                             >  1:1662459/1‑48 (MQ=255)
                                                                   cTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTCCATg      >  1:2494708/1‑70 (MQ=255)
                                                                   cTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTCCATg      <  1:2268888/70‑1 (MQ=255)
                                                                      gCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTCCATGc     >  1:618601/1‑68 (MQ=255)
                                                                      gCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTCCATGCa    <  1:282014/69‑1 (MQ=255)
                                                                      gCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTCCATGCa    <  1:3208877/69‑1 (MQ=255)
                                                                      gCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTCCATGCAgc  >  1:1251075/1‑71 (MQ=255)
                                                                      |                                                                      
AGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGGGAACAATCTGTCCATGCAGC  >  minE/242793‑242933

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: