Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,790 A→C 6.8% T302T (ACT→ACG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,7900AC 6.8% 188.5 / 7.6 74T302T (ACT→ACGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base C (3/2);  total (39/35)
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.98e-01

TTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCA  >  minE/2918722‑2918859
                                                                    |                                                                     
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)
          gCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATa                                                           >  1:2586964/1‑71 (MQ=255)
           cATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCAtcc                                                          >  1:3815328/1‑69 (MQ=255)
             tACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATAc                                                          <  1:3482601/69‑1 (MQ=255)
              ccGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCAcct                                                                >  1:456762/2‑60 (MQ=255)
              ccGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGc                                                       >  1:1610999/2‑71 (MQ=255)
              aCGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc                                                       >  1:3433599/1‑71 (MQ=255)
               cGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGcc                                                      <  1:2703955/71‑1 (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)
                    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:1925411/1‑71 (MQ=255)
                                  aaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTcc                                   >  1:74187/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:3884315/69‑1 (MQ=255)
                                             gcggcAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTc                          <  1:2541658/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:98349/1‑71 (MQ=255)
                                                      tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGcgt               >  1:2419196/1‑71 (MQ=255)
                                                       cGAAGAACTCATCAGTGACTTCCATACCGCCATGTAGGCCaaa                                          >  1:2632021/1‑43 (MQ=255)
                                                        gAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttcttc                                    <  1:3394476/48‑1 (MQ=255)
                                                         aaGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGcgtc              <  1:464896/69‑1 (MQ=255)
                                                          aGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCa                               >  1:807925/1‑51 (MQ=255)
                                                          aGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGtt           >  1:2142419/1‑71 (MQ=255)
                                                           gACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAATCTTCCGCCATTTTCTTGCCAAGCGTCGtt           >  1:1492204/1‑70 (MQ=255)
                                                             ccTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTgg        >  1:981076/1‑71 (MQ=255)
                                                                catcaGTGACTTCCATACCGCCATGTAGGCCAAATTCTTcc                                   >  1:1633500/1‑41 (MQ=255)
                                                                catcaGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGgtcg     >  1:2231690/1‑71 (MQ=255)
                                                                catcaGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGgtcg     >  1:1415679/1‑71 (MQ=255)
                                                                catcaGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGgtcg     >  1:1917534/1‑71 (MQ=255)
                                                                   caGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAg                  >  1:1645646/1‑55 (MQ=255)
                                                                   caGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGgtc      >  1:1883977/1‑67 (MQ=255)
                                                                   caGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCa  >  1:2177474/1‑71 (MQ=255)
                                                                    |                                                                     
TTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCA  >  minE/2918722‑2918859

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: