Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,760 A→G 12.6% I312I (ATT→ATC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,7600AG12.6% 204.3 / 18.3 79I312I (ATT→ATCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base G (7/3);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCG  >  minE/2918690‑2918827
                                                                      |                                                                   
attattTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACa                                                                     <  1:2201761/71‑1 (MQ=255)
attattTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACa                                                                     <  1:1452984/71‑1 (MQ=255)
   attTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAAt                                                                   >  1:1688317/1‑70 (MQ=255)
      tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCt                                                                <  1:3499375/70‑1 (MQ=255)
      tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                               >  1:1084860/1‑71 (MQ=255)
      tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                               >  1:795772/1‑71 (MQ=255)
      tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                               >  1:3032873/1‑71 (MQ=255)
       aCTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                               <  1:609998/70‑1 (MQ=255)
         tGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                                                            >  1:1903999/1‑71 (MQ=255)
         caaGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgcc                                                            >  1:527391/3‑70 (MQ=255)
          aaGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgccg                                                           >  1:1215440/2‑69 (MQ=255)
           aGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcgg                                                           >  1:1836579/1‑70 (MQ=255)
           aGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcgg                                                           <  1:3623585/70‑1 (MQ=255)
               tCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAt                                                      <  1:2444571/71‑1 (MQ=255)
                 gcgACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTc                                                    >  1:1460201/1‑71 (MQ=255)
                  cgACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTc                                                    >  1:1284918/1‑70 (MQ=255)
                    aCCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCg                                                   >  1:3350507/1‑69 (MQ=255)
                     tcATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGAt                                                                   >  1:2305279/2‑52 (MQ=255)
                      cATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgccgccg                                                        >  1:511659/1‑60 (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)
                                         tGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                                                            <  1:2991960/39‑1 (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)
                                             tACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAtt                                                     <  1:519165/42‑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:3720328/70‑1 (MQ=255)
                                                    ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                  <  1:3396146/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:3643457/71‑1 (MQ=255)
                                                         ccTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa            <  1:3791017/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:196137/1‑71 (MQ=255)
                                                                aaaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttctt     <  1:101615/71‑1 (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)
                                                                      |                                                                   
ATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCG  >  minE/2918690‑2918827

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: