Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 244,761 C→T 39.1% G221R (GGG→AGG)  sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE244,7610CT39.1% 50.9 / 71.5 69G221R (GGG→AGG) sbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base T (13/14);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

ATGGGGGAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGT  >  minE/244697‑244829
                                                                |                                                                    
aTGGGGGAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAgcgc                                                                >  1:3979400/1‑71 (MQ=255)
aTGGGGGAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAgcgc                                                                <  1:3684698/71‑1 (MQ=255)
 tGGGGGAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAgcgc                                                                >  1:710494/1‑70 (MQ=255)
 tGGGGGAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGcgcg                                                               <  1:3381322/71‑1 (MQ=255)
    gggAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATg                                                            <  1:58207/71‑1 (MQ=255)
    gggAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATg                                                            <  1:2196045/71‑1 (MQ=255)
      gAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTa                                                          <  1:2556003/71‑1 (MQ=255)
      gAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTa                                                          <  1:3583324/71‑1 (MQ=255)
       aGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAg                                                         >  1:320207/1‑71 (MQ=255)
        gCCGCAATAGCGAACATGTTCCATGTCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGt                                                        >  1:4237853/1‑71 (MQ=255)
           gCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCgg                                                     >  1:1912817/1‑71 (MQ=255)
             aaTAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCt                                                   >  1:4232721/1‑71 (MQ=255)
                aGCGAACATGTTCCATGCCGCCATTAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCtggt                                                >  1:4614247/1‑71 (MQ=255)
                aGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCtggt                                                >  1:4608749/1‑71 (MQ=255)
                 gCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAg                                                         >  1:1621678/1‑61 (MQ=255)
                 gCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAg                                                         >  1:2118105/1‑61 (MQ=255)
                      cATGTTCCATGCCGCCAATAATCTGTGCGCGGTTAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAg                                          <  1:2305450/71‑1 (MQ=255)
                      cATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAg                                          <  1:4578208/71‑1 (MQ=255)
                      cATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAg                                          >  1:3057552/1‑71 (MQ=255)
                        tGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGtt                                        <  1:1636096/71‑1 (MQ=255)
                        tGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGt                                         <  1:3763664/70‑1 (MQ=255)
                        tGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAg                                                         <  1:2507124/54‑1 (MQ=255)
                        tGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGtt                                        <  1:955570/71‑1 (MQ=255)
                        tGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGtt                                        <  1:2049821/71‑1 (MQ=255)
                         gTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCtggtg                                               >  1:270619/1‑63 (MQ=255)
                         gTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGttt                                       >  1:276112/1‑71 (MQ=255)
                          ttCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGttt                                       <  1:4276158/70‑1 (MQ=255)
                          ttCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATg                                                            <  1:4287465/49‑1 (MQ=255)
                          ttCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAg                                                         <  1:2116602/52‑1 (MQ=255)
                             cATGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTtgt                                    <  1:3063905/70‑1 (MQ=255)
                              aTGCCGCCAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGc                                  >  1:1377725/1‑71 (MQ=255)
                               tGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATg                                                            <  1:1139409/44‑1 (MQ=255)
                                gccgccAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATg                                                            <  1:1133858/43‑1 (MQ=255)
                                 ccgccAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGcc                                 >  1:2350170/1‑69 (MQ=255)
                                 ccgccAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCgg                               <  1:3996868/71‑1 (MQ=255)
                                   gccAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGa                              >  1:1493106/1‑70 (MQ=255)
                                     cAATAATCTGTGCGCGGTGAATATGCCTGAGcgcg                                                               <  1:3229520/35‑1 (MQ=255)
                                     cAATAATCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAAc                           <  1:3507617/71‑1 (MQ=255)
                                      aataatCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTtgtg                                   >  1:2271962/1‑62 (MQ=255)
                                        taatCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGaa                                            <  1:3514358/51‑1 (MQ=255)
                                        taatCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAAcgcg                        <  1:3651241/71‑1 (MQ=255)
                                         aatCTGTGCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGtt                                        >  1:1424923/1‑54 (MQ=255)
                                          atCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCtg                                                  >  1:4497600/1‑43 (MQ=255)
                                               tgCGCGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAg                   <  1:3448480/69‑1 (MQ=255)
                                               tgCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCg                 >  1:375101/1‑71 (MQ=255)
                                               tgCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCg                 <  1:3077347/71‑1 (MQ=255)
                                                gcgcGGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGttt                                       <  1:2037544/48‑1 (MQ=255)
                                                 cgcgGTGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTg               >  1:4487220/1‑71 (MQ=255)
                                                  gcgGTGAATATGCCCGAGCGCGATGTGGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGc              <  1:2855936/71‑1 (MQ=255)
                                                  gcgGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGc              >  1:2505367/1‑71 (MQ=255)
                                                   cgGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGc              >  1:4372218/1‑70 (MQ=255)
                                                   cgGTGAATATGCCCGAGCGCGATGTAATCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGcc             <  1:1695029/71‑1 (MQ=255)
                                                      tGAATATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCaa           >  1:478308/1‑70 (MQ=255)
                                                      tGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGtt                                        <  1:2906398/41‑1 (MQ=255)
                                                        aaTATGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAatat        <  1:2788139/71‑1 (MQ=255)
                                                        aaTATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGtttt                                      >  1:4235935/1‑41 (MQ=255)
                                                         atatGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAatat        >  1:3660931/1‑70 (MQ=255)
                                                         atatGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAatat        >  1:85986/1‑70 (MQ=255)
                                                         atatGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAAtata       <  1:4157659/71‑1 (MQ=255)
                                                         atatGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAAtata       <  1:2197829/71‑1 (MQ=255)
                                                         atatGCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAAtata       <  1:537717/71‑1 (MQ=255)
                                                         atatGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAatat        >  1:161868/1‑70 (MQ=255)
                                                          tatGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGc                                  >  1:2128687/1‑43 (MQ=255)
                                                          tatGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCa            <  1:565045/65‑1 (MQ=255)
                                                          tatGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAAtata       <  1:941429/70‑1 (MQ=255)
                                                          tatGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATaa      <  1:166549/71‑1 (MQ=255)
                                                            tGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAAtata       <  1:4423294/68‑1 (MQ=255)
                                                             gCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATaa      >  1:3841255/1‑68 (MQ=255)
                                                             gCCTGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATaa      >  1:546742/1‑68 (MQ=255)
                                                              cccGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATg   <  1:3875163/70‑1 (MQ=255)
                                                              cccGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGt  >  1:329979/1‑71 (MQ=255)
                                                               ccGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGcc                                 >  1:1801636/1‑39 (MQ=255)
                                                                cGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGa                              <  1:2650914/41‑1 (MQ=255)
                                                                cGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCa            <  1:1424634/59‑1 (MQ=255)
                                                                cGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCa            <  1:1111497/59‑1 (MQ=255)
                                                                |                                                                    
ATGGGGGAGCCGCAATAGCGAACATGTTCCATGCCGCCAATAATCTGTGCGCGGTGAATATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGT  >  minE/244697‑244829

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: