Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,778 G→A 6.3% F306F (TTC→TTT argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,7780GA 6.3% 244.2 / 4.6 79F306F (TTC→TTTargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (40/34);  new base A (3/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

ACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGT  >  minE/2918710‑2918846
                                                                    |                                                                    
aCCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCg                                                                      >  1:3350507/1‑69 (MQ=255)
   aTCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCaaaca                                                                  <  1:2276540/70‑3 (MQ=255)
    tCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGa                                                                     <  1:2097934/66‑1 (MQ=255)
       ccGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtca                                                             <  1:538164/71‑1 (MQ=255)
        cGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa                                                                  >  1:1843292/1‑65 (MQ=255)
         gCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcat                                                            >  1:2259504/1‑70 (MQ=255)
         gCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcat                                                            >  1:2800896/1‑70 (MQ=255)
            ttGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAg                                                         >  1:2546262/1‑70 (MQ=255)
             tGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCACATCAGt                                                        >  1:2475180/1‑70 (MQ=255)
             tGATAGTATGCATACGGTTTTCCACCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGt                                                        <  1:1433903/70‑1 (MQ=255)
                  gTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                                                  <  1:2289895/71‑1 (MQ=255)
                   tATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTcc                                                 <  1:3665075/71‑1 (MQ=255)
                   tATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                                                  >  1:2370965/1‑70 (MQ=255)
                     tGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCAt                                               <  1:2537598/71‑1 (MQ=255)
                      gCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGaa                                                                    >  1:3323714/1‑49 (MQ=255)
                      gCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATa                                              >  1:2586964/1‑71 (MQ=255)
                       cATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCAtcc                                             >  1:3815328/1‑69 (MQ=255)
                       cATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGa                                                                     <  1:231261/47‑1 (MQ=255)
                         tACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATAc                                             <  1:3482601/69‑1 (MQ=255)
                          ccGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCAcct                                                   >  1:456762/2‑60 (MQ=255)
                          ccGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGc                                          >  1:1610999/2‑71 (MQ=255)
                          aCGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa                                                                  >  1:627506/1‑47 (MQ=255)
                          aCGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc                                          >  1:3433599/1‑71 (MQ=255)
                           cGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGcc                                         <  1:2703955/71‑1 (MQ=255)
                           cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGAcc                                                                >  1:1417919/1‑48 (MQ=255)
                           cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcatca                                                          >  1:880587/1‑54 (MQ=255)
                           cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATa                                              >  1:483246/1‑66 (MQ=255)
                           cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGcc                                         <  1:985325/71‑1 (MQ=255)
                            ggTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCa                                        >  1:1931823/1‑71 (MQ=255)
                             gTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCAt                                       <  1:2494147/71‑1 (MQ=255)
                              ttttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATg                                      <  1:3403002/71‑1 (MQ=255)
                                ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGaa                                                                    >  1:1991256/1‑39 (MQ=255)
                                ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATg                                      >  1:883096/1‑69 (MQ=255)
                                ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                                     <  1:3396146/70‑1 (MQ=255)
                                ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                                     <  1:3720328/70‑1 (MQ=255)
                                 tCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa                                    <  1:3683894/70‑1 (MQ=255)
                                   cGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                                                  <  1:2049051/54‑1 (MQ=255)
                                   cGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                                                  <  1:860264/54‑1 (MQ=255)
                                     ccTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGc                                 <  1:3160204/69‑1 (MQ=255)
                                     ccTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa                               <  1:3791017/71‑1 (MQ=255)
                                     ccTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa                               <  1:3643457/71‑1 (MQ=255)
                                       tGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACtt                                                   <  1:3791737/49‑1 (MQ=255)
                                         aTCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCaaa                             <  1:3074562/69‑1 (MQ=255)
                                         aTCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAtt                           >  1:3867180/1‑71 (MQ=255)
                                          tCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCAt                                       >  1:190238/1‑58 (MQ=255)
                                           cAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAgg                                  >  1:631627/1‑62 (MQ=255)
                                           cAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGTCCAAAtt                           >  1:1350505/1‑69 (MQ=255)
                                           cAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa                               <  1:2587262/65‑1 (MQ=255)
                                           cAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttct                         >  1:1741476/1‑71 (MQ=255)
                                            aaaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa                                    >  1:1667542/1‑59 (MQ=255)
                                            aaaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttctt                        <  1:101615/71‑1 (MQ=255)
                                            aaaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttctt                        >  1:196137/1‑71 (MQ=255)
                                            aaaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttct                         <  1:1540500/70‑1 (MQ=255)
                                             aaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAtt                           >  1:1944187/1‑67 (MQ=255)
                                              aaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTcc                      >  1:74187/1‑71 (MQ=255)
                                              aaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTcc                      >  1:1925411/1‑71 (MQ=255)
                                               aaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCaa                              <  1:1472181/62‑1 (MQ=255)
                                               aaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTcc                      >  1:1214096/1‑70 (MQ=255)
                                               aaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCg                     <  1:2405858/71‑1 (MQ=255)
                                                aCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAgg                                  >  1:487055/1‑57 (MQ=255)
                                                 cAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAgg                                  >  1:746293/1‑56 (MQ=255)
                                                   aTGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCAt                 <  1:288285/71‑1 (MQ=255)
                                                    tGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCAtt                <  1:2109237/71‑1 (MQ=255)
                                                       tGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTc             >  1:2847188/1‑71 (MQ=255)
                                                         gcggcAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTc             <  1:2541658/69‑1 (MQ=255)
                                                         gcggcAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTc             <  1:3884315/69‑1 (MQ=255)
                                                         gcggcAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCtt           <  1:1210628/71‑1 (MQ=255)
                                                          cggcAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTg          >  1:756714/1‑71 (MQ=255)
                                                          cggcAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTg          >  1:3817298/1‑71 (MQ=255)
                                                          cggcAGATTCGAAGACCTCATAAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCtt           >  1:2976809/1‑70 (MQ=255)
                                                             cAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa                                    >  1:1807754/1‑42 (MQ=255)
                                                             cAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGcc        >  1:3310994/1‑70 (MQ=255)
                                                              aGATTCGAAGACCTCATCAGTGACTTCCATACCGCCAt                                       >  1:1284664/1‑38 (MQ=255)
                                                                aTTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCg                     <  1:2086230/54‑1 (MQ=255)
                                                                aTTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCg                     <  1:964462/54‑1 (MQ=255)
                                                                 ttCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa                               <  1:417495/43‑1 (MQ=255)
                                                                 ttCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa                               <  1:1876398/43‑1 (MQ=255)
                                                                  tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCa       >  1:2414540/1‑66 (MQ=255)
                                                                  tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGcgt  >  1:2419196/1‑71 (MQ=255)
                                                                  tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGcgt  >  1:98349/1‑71 (MQ=255)
                                                                   cGAAGAACTCATCAGTGACTTCCATACCGCCATGTAGGCCaaa                             >  1:2632021/1‑43 (MQ=255)
                                                                    gAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttcttc                       <  1:3394476/48‑1 (MQ=255)
                                                                    |                                                                    
ACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGT  >  minE/2918710‑2918846

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: