Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,334,615 G→A 53.1% L72L (CTG→TTG)  nudD ← GDP‑mannose mannosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,334,6150GA53.1% ‑0.9 / 67.5 49L72L (CTG→TTG) nudDGDP‑mannose mannosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base G (12/11);  new base A (14/12);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CCAGAGAAGTTATCGTCATAAAAGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGTCTTTCTGCACGCGCCCTC  >  minE/1334545‑1334683
                                                                      |                                                                    
ccAGAGAAGTTATCGTCATAAAAGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAg                                                                      >  1:3005836/1‑71 (MQ=255)
  agagAATTTATCGTCATAAAAGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAAcc                                                                    <  1:1246343/71‑1 (MQ=255)
         ttATCGTCATAAAAGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAAccc                                                                   <  1:3206323/65‑1 (MQ=255)
           aTCGTCATAAAAGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCg                                                           >  1:1495666/1‑71 (MQ=255)
              gTCATAAAAGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCa                                                        >  1:376374/1‑71 (MQ=255)
              gTCATAAAAGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGcc                                                         <  1:3224984/70‑1 (MQ=255)
                      aGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGc                                                >  1:180522/1‑71 (MQ=255)
                             tctgACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCa                                          <  1:3761212/67‑1 (MQ=255)
                             ccAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCa                                                                 >  1:2342348/1‑47 (MQ=255)
                             ccAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCa                                          <  1:1345481/70‑1 (MQ=255)
                             ccAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCa                                          <  1:3455260/70‑1 (MQ=255)
                             ccAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCa                                          <  1:1888888/70‑1 (MQ=255)
                               aGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAAt                                       <  1:2219515/71‑1 (MQ=255)
                                 acacCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGc                                     <  1:604847/71‑1 (MQ=255)
                                 acacCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCa                                          <  1:1729528/66‑1 (MQ=255)
                                 acacCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATg                                      >  1:1940163/1‑70 (MQ=255)
                                    ccGTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATc                                                      <  1:3670225/51‑1 (MQ=255)
                                     cGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCt                                 <  1:1317938/71‑1 (MQ=255)
                                     cGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCt                                 <  1:3352234/71‑1 (MQ=255)
                                      gTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCAt                                                       >  1:3807535/1‑48 (MQ=255)
                                      gTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGc                                             <  1:2912221/58‑1 (MQ=255)
                                      gTAAAACTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCtt                                >  1:56535/1‑71 (MQ=255)
                                       tAAAACTGGCCTGCTGTTATCGGCAAACGCAGccc                                                                   >  1:501594/1‑35 (MQ=255)
                                       tAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATg                                      >  1:3926491/1‑64 (MQ=255)
                                       tAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATg                                      >  1:2860575/1‑64 (MQ=255)
                                        aaaaCTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTc                               >  1:1647896/1‑70 (MQ=255)
                                         aaaCTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTc                               >  1:1962450/1‑69 (MQ=255)
                                         aaaCTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCAGCCATCGTCAGCCGCTCAAATGCGGCTTc                               >  1:3293720/1‑69 (MQ=255)
                                         aaaCTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCa                                                  >  1:3098310/1‑50 (MQ=255)
                                         aaaCTGGCCTGCTGTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTc                               >  1:1694712/1‑69 (MQ=255)
                                            cTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCg                          >  1:1422185/1‑71 (MQ=255)
                                                 ctgctgTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCa                                          <  1:2775325/50‑1 (MQ=255)
                                                  tgctgTTATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCaa                                         <  1:605929/50‑1 (MQ=255)
                                                   gctgTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGtt                        <  1:2535670/66‑1 (MQ=255)
                                                    ctgTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGt                    <  1:1667902/69‑1 (MQ=255)
                                                       ttATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCg                                    <  1:4018211/50‑1 (MQ=255)
                                                       ttATCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGc                                  >  1:123472/1‑52 (MQ=255)
                                                         aTCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCaa                                         >  1:4083304/1‑43 (MQ=255)
                                                         aTCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCaa                                         >  1:1260950/1‑43 (MQ=255)
                                                         aTCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGttt                       >  1:3624114/1‑61 (MQ=255)
                                                         aTCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGTCTTTCTg             >  1:713287/1‑71 (MQ=255)
                                                         aTCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGTCTTTCTg             >  1:4019444/1‑71 (MQ=255)
                                                         aTCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGTCTTTCTg             >  1:1861523/1‑71 (MQ=255)
                                                         aTCGGCAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCCGCTTCCAg                            >  1:528576/1‑56 (MQ=39)
                                                          tCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGTCTTTCTg             <  1:2599361/70‑1 (MQ=255)
                                                            ggCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGTCTTTCTGCa           <  1:2512172/70‑1 (MQ=255)
                                                              cAAACGCAACCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGc                           >  1:1738966/1‑52 (MQ=255)
                                                               aaaCGCAACCCCAGTTCCGCCATCGTCAGCCGCTCa                                          <  1:2387532/36‑1 (MQ=255)
                                                                  cGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGc                           >  1:3927506/1‑48 (MQ=255)
                                                                    cAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGTCTTTCTGCACGCGCCCTc  <  1:3365385/71‑1 (MQ=255)
                                                                      |                                                                    
CCAGAGAAGTTATCGTCATAAAAGTGCTGCCAGACACCGTAAAACTGGCCTGCTGTTATCGGCAAACGCAGCCCCAGTTCCGCCATCGTCAGCCGCTCAAATGCGGCTTCCAGCGTTTCGTCTTTCTGCACGCGCCCTC  >  minE/1334545‑1334683

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: