Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 244,476 A→G 14.0% Y316H (TAT→CAT)  sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE244,4760AG14.0% 258.6 / 31.5 100Y316H (TAT→CAT) sbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (47/39);  new base G (6/8);  total (53/47)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 108 total aligned reads displayed.

TACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGC  >  minE/244408‑244546
                                                                    |                                                                      
tACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACt                                                                      <  1:1075896/71‑1 (MQ=255)
tACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACt                                                                      <  1:843047/71‑1 (MQ=255)
  ccAGCAATACTTCGACAGTCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtca                                                                    <  1:2859039/71‑1 (MQ=255)
  ccAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtca                                                                    >  1:876967/1‑71 (MQ=255)
  ccAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtca                                                                    >  1:2710620/1‑71 (MQ=255)
  ccAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtc                                                                     <  1:2052306/70‑1 (MQ=255)
   cAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtcat                                                                   >  1:3792589/1‑71 (MQ=255)
   cAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtcat                                                                   >  1:3374143/1‑71 (MQ=255)
     gCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATa                                                                        >  1:2486053/1‑64 (MQ=255)
      gaaTACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCag                                                                <  1:1591959/70‑1 (MQ=255)
      cAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCag                                                                <  1:1149630/71‑1 (MQ=255)
       aaTACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCag                                                                >  1:2826043/1‑70 (MQ=255)
         tACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtcat                                                                   <  1:1549083/65‑1 (MQ=255)
          aCTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATGCTCATCagta                                                              >  1:1960767/1‑69 (MQ=255)
           cTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATa                                                                        >  1:1134239/1‑58 (MQ=255)
           cTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATa                                                                        >  1:1129435/1‑58 (MQ=255)
           cTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATGCTCATCagtagt                                                            >  1:3333949/1‑70 (MQ=255)
           cTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtcat                                                                   >  1:3234852/1‑63 (MQ=255)
           cTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagtagt                                                            <  1:3178716/70‑1 (MQ=255)
           cTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagtagt                                                            >  1:2598560/1‑70 (MQ=255)
            ttCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTAATCAGTAGTGa                                                          >  1:4146100/1‑71 (MQ=255)
             tCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGAt                                                         >  1:1779755/1‑71 (MQ=255)
               gACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGAtt                                                        <  1:3832427/70‑1 (MQ=255)
                aCAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTc                                                      <  1:2728953/71‑1 (MQ=255)
                  aGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagta                                                              >  1:48746/1‑61 (MQ=255)
                   ggCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATGCTCATCAGTAGTGAttt                                                       >  1:3024507/1‑67 (MQ=255)
                   ggCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTg                                                           >  1:320866/1‑63 (MQ=255)
                      aaTGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtcatca                                                                 >  1:3727745/1‑54 (MQ=255)
                      aaTGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagtagt                                                            <  1:2995156/59‑1 (MQ=255)
                      aaTGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagta                                                              <  1:2788718/57‑1 (MQ=255)
                      aaTGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATc                                                >  1:3870722/1‑71 (MQ=255)
                      aaTGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATc                                                >  1:888780/1‑71 (MQ=255)
                        tGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGa                                                    >  1:2666647/1‑65 (MQ=255)
                        tGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGa                                                    >  1:3482946/1‑65 (MQ=255)
                         gATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccag                                             >  1:1841284/1‑71 (MQ=255)
                          aTTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATGCTCATCAGTAGTGATTTCGATATccagc                                            <  1:387882/71‑1 (MQ=255)
                           ttCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATGCtcatca                                                                 >  1:166240/1‑49 (MQ=255)
                           ttCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATGCtcatc                                                                  >  1:3896111/1‑48 (MQ=255)
                            tCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATa                                                                        >  1:46251/1‑41 (MQ=255)
                             cGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccagcca                                          <  1:1904187/70‑1 (MQ=255)
                             cGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccagcca                                          <  1:2965016/70‑1 (MQ=255)
                             cGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccagcca                                          <  1:904560/70‑1 (MQ=255)
                             cGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccagc                                            >  1:1052797/1‑68 (MQ=255)
                             cGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccagc                                            >  1:3882146/1‑68 (MQ=255)
                             cGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATccagc                                            >  1:2199723/1‑68 (MQ=255)
                              ggTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtca                                                                    >  1:2053598/1‑43 (MQ=255)
                              ggTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGa                                        <  1:2175034/71‑1 (MQ=255)
                               gTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATa                                                                        >  1:1029405/1‑38 (MQ=255)
                               gTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagtagt                                                            <  1:2513375/50‑1 (MQ=255)
                                ttAATGCCTGGATTTTGCGCTGAATATCATGCAGATACtc                                                                     <  1:1389640/40‑1 (MQ=255)
                                ttAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagt                                                               <  1:3746459/46‑1 (MQ=255)
                                ttAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCagt                                                               <  1:1597337/46‑1 (MQ=255)
                                ttAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGAtt                                                        <  1:771266/53‑1 (MQ=255)
                                ttAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGAtt                                                        <  1:162491/53‑1 (MQ=255)
                                ttAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCTGACa                                      <  1:880702/71‑1 (MQ=255)
                                ttAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACa                                      <  1:1491003/71‑1 (MQ=255)
                                 tAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCATCCAGACAg                                     >  1:3530232/1‑71 (MQ=255)
                                 tAATGCCTGGATTTTGCGCTGAATATCATGAAGATACTCATCAGTAGTGa                                                          >  1:919456/1‑50 (MQ=255)
                                   aTGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGa                                                          >  1:3428054/1‑48 (MQ=255)
                                   aTGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGAt                                                         <  1:2133884/49‑1 (MQ=255)
                                    tGCCTGGATTTTGCGCTGAGTATCATGCAGATACtcat                                                                   <  1:1128998/38‑1 (MQ=255)
                                    tGCCTGGATTTTGCGCTGAATATCATGCAGATGCTCATCAGTAGTGATTTc                                                      <  1:1761861/51‑1 (MQ=255)
                                    tGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGa                                        >  1:927479/1‑65 (MQ=255)
                                    tGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGa                                        >  1:2405699/1‑65 (MQ=255)
                                    tGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTg                                  >  1:1399627/1‑71 (MQ=255)
                                      ccTGGATTTTGCGCTGAATATCATGCAGATACTCATCagta                                                              >  1:188849/1‑41 (MQ=255)
                                      ccTGGATTTTGCGCTGAATATCATGCAGATACTCATCagta                                                              >  1:35335/1‑41 (MQ=255)
                                       cTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCt                               >  1:3560318/1‑71 (MQ=255)
                                       cTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCt                               <  1:2664654/71‑1 (MQ=255)
                                        tGGATTTTGCGCTGAATATCATGCAGATGCTCATCAGTAGTGATTTCGATATCCAGCCAGAc                                       <  1:2480836/62‑1 (MQ=255)
                                        tGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATcca                                              >  1:3082909/1‑55 (MQ=255)
                                         ggATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGa                                                          >  1:3462340/1‑42 (MQ=255)
                                          gATTTTGCGCTGAATATCATGCAGATGCTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGc                                >  1:1650136/1‑67 (MQ=255)
                                          gATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCt                            >  1:3771451/1‑71 (MQ=255)
                                           aTTTTGCGCTGAATATCATGCAGATACTCATCagta                                                              <  1:1418905/36‑1 (MQ=255)
                                           aTTTTGCGCTGAATATCATGCAGATACTCATCagta                                                              <  1:1549274/36‑1 (MQ=255)
                                           aTTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCt                            <  1:1754765/70‑1 (MQ=255)
                                           aTTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTg                           >  1:952623/1‑71 (MQ=255)
                                            ttttGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGc                          <  1:287876/71‑1 (MQ=255)
                                               tGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGa                        >  1:1270186/1‑70 (MQ=255)
                                                gcgcTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGa                        <  1:3326516/69‑1 (MQ=255)
                                                gcgcTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGAt                       <  1:631749/70‑1 (MQ=255)
                                                 cgcTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATAc                     <  1:1066053/71‑1 (MQ=255)
                                                   cTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACAt                   <  1:2449510/71‑1 (MQ=255)
                                                    tGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGa                        <  1:2040610/65‑1 (MQ=255)
                                                    tGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACAt                   >  1:2816433/1‑70 (MQ=255)
                                                    tGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATc                  >  1:1993783/1‑71 (MQ=255)
                                                     gAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATc                  >  1:1670402/1‑70 (MQ=255)
                                                      aaTATCATGCAGATGCTCATCAGTAGTGATTTCGATATccagcca                                          <  1:3015022/45‑1 (MQ=255)
                                                      aaTATCATGCAGATGCTCATCAGTAGTGATTTCGATATccagcca                                          <  1:1537936/45‑1 (MQ=255)
                                                      aaTATCATGCAGATGCTCATCAGTAGTGATTTCGATATccagcca                                          <  1:1542828/45‑1 (MQ=255)
                                                      aaTATCATGCAGATGCTCATCAGTAGTGATTTCGATATccagcca                                          <  1:2746914/45‑1 (MQ=255)
                                                        tatCATGCAGATGCTCATCAGTAGTGATTTCGATATCCAGCCAGAc                                       <  1:210444/46‑1 (MQ=255)
                                                           cATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAgg                                    >  1:912566/1‑46 (MQ=255)
                                                            aTGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACt          >  1:639869/1‑71 (MQ=255)
                                                             tGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGAc                                       >  1:531244/1‑41 (MQ=255)
                                                                 gATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTcc     <  1:3009518/71‑1 (MQ=255)
                                                                   tACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAg   >  1:3529391/1‑71 (MQ=255)
                                                                    aCTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGa                        >  1:3065529/1‑49 (MQ=255)
                                                                    aCTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGc  >  1:2463646/1‑71 (MQ=255)
                                                                    |                                                                      
TACCAGCAATACTTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGC  >  minE/244408‑244546

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: