Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,250,792 A→C 54.1% L167L (CTT→CTG ruvC ← component of RuvABC resolvasome, endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,250,7920AC54.1% 42.2 / 64.1 61L167L (CTT→CTGruvCcomponent of RuvABC resolvasome, endonuclease
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base C (13/20);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCA  >  minE/1250725‑1250862
                                                                   |                                                                      
gagggagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAg                                                                      <  1:132145/70‑1 (MQ=255)
  gggagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGtt                                                                   <  1:3376865/71‑1 (MQ=255)
   ggagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTc                                                                  <  1:1267935/71‑1 (MQ=255)
   ggagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTc                                                                  <  1:1546080/71‑1 (MQ=255)
    gagTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCa                                                                 <  1:3439476/71‑1 (MQ=255)
           aaCGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg                                                                     >  1:3366016/1‑60 (MQ=255)
           aaCGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGa                                                          >  1:897646/1‑71 (MQ=255)
            aCGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGAt                                                         >  1:704328/1‑71 (MQ=255)
             cGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGAtt                                                        <  1:3459006/71‑1 (MQ=255)
                cTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAg                                                                      >  1:1404709/1‑54 (MQ=255)
                 tCAGCCGGAAGTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCt                                                    >  1:1593720/1‑71 (MQ=255)
                 tCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGAtt                                                        <  1:2524925/67‑1 (MQ=255)
                 tCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCt                                                    >  1:949991/1‑71 (MQ=255)
                   aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCa                                                  >  1:303735/1‑71 (MQ=255)
                   aGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTc                                                   <  1:355047/70‑1 (MQ=255)
                    gCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCAt                                                 >  1:284504/1‑71 (MQ=255)
                    gCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCAt                                                 >  1:2718212/1‑71 (MQ=255)
                     ccGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCa                                                  <  1:3232482/69‑1 (MQ=255)
                         aaCTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCcgc                                                            <  1:271920/55‑1 (MQ=255)
                         aaCTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCa                                            >  1:2751237/1‑71 (MQ=255)
                         aaCTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCa                                            <  1:1568360/71‑1 (MQ=255)
                         aaCTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCa                                            >  1:2158845/1‑71 (MQ=255)
                              aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCg                                         <  1:2894932/69‑1 (MQ=255)
                              aCCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCg                                         <  1:236572/69‑1 (MQ=255)
                               ccGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCa                                       <  1:577593/70‑1 (MQ=255)
                                cGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCAtt                                     <  1:50235/71‑1 (MQ=255)
                                  aGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                   >  1:247452/1‑71 (MQ=255)
                                    gCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCa                                            <  1:135771/60‑1 (MQ=255)
                                       gTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGAtt                                                        <  1:1624996/45‑1 (MQ=255)
                                        tataACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGAtt                                                        <  1:844813/44‑1 (MQ=255)
                                        tataACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCa                                            <  1:304336/56‑1 (MQ=255)
                                          taACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGc                                             <  1:732766/53‑1 (MQ=255)
                                          taACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAAc                            <  1:2096687/70‑1 (MQ=255)
                                              ttAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTaa                             <  1:2111905/65‑1 (MQ=255)
                                              ttAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTaa                             <  1:137700/65‑1 (MQ=255)
                                              ttAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTaa                             <  1:866174/65‑1 (MQ=255)
                                              ttAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTaa                             <  1:1497819/65‑1 (MQ=255)
                                              ttAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGt                          <  1:717332/68‑1 (MQ=255)
                                               tAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                        >  1:1777447/1‑69 (MQ=255)
                                               tAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTgg                        >  1:3292912/1‑69 (MQ=255)
                                                aaCGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCt                                   >  1:1141166/1‑57 (MQ=255)
                                                aaCGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACg                           >  1:1758371/1‑65 (MQ=255)
                                                aaCGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTg                         >  1:1792170/1‑67 (MQ=255)
                                                aaCGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGc                       <  1:2696166/69‑1 (MQ=255)
                                                 aCGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGt                    <  1:3719878/71‑1 (MQ=255)
                                                  cGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTg                   >  1:2908209/1‑71 (MQ=255)
                                                  cGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTg                   >  1:3495484/1‑71 (MQ=255)
                                                  cGCAGTCGCCCTCTCGCAAGGTTAAGCCGCGATTCGCTCATCTGCATAGCATTCTGACTAACGTGGCAGTg                   >  1:1143060/1‑71 (MQ=255)
                                                       tCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCa                                            <  1:1266173/41‑1 (MQ=255)
                                                       tCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCa                                            <  1:1924988/41‑1 (MQ=255)
                                                       tCGCCCTCTCGCCAGCTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTTACTAACGTGGCAGTGGGTGa              <  1:1116716/71‑1 (MQ=255)
                                                       tCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAAc                            <  1:785167/57‑1 (MQ=255)
                                                       tCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGa              <  1:1402318/71‑1 (MQ=255)
                                                         gCCCTCTCGCCAGGTTCAGCCGCGATTCGCTAATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGa              >  1:1953268/1‑69 (MQ=255)
                                                           ccTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAg           >  1:1670295/1‑70 (MQ=255)
                                                           ccTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGc          <  1:1095557/71‑1 (MQ=255)
                                                            ctctCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGc          >  1:479333/1‑70 (MQ=255)
                                                               tcGCCAGGTTCAGCCGCGATTTGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCaa        <  1:2091331/69‑1 (MQ=255)
                                                               tcGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATa            >  1:1565399/1‑65 (MQ=255)
                                                                 gCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATc      <  1:2645832/69‑1 (MQ=255)
                                                                  cAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGc    >  1:3544342/1‑70 (MQ=255)
                                                                   aaGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGc    >  1:1712060/1‑69 (MQ=255)
                                                                   aaGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGc    >  1:851149/1‑69 (MQ=255)
                                                                   aaGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCa  <  1:704963/71‑1 (MQ=255)
                                                                   |                                                                      
GAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCAAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCA  >  minE/1250725‑1250862

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: