Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,757 A→C 14.3% V313V (GTT→GTG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,7570AC14.3% 137.5 / 24.6 70V313V (GTT→GTGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base C (7/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCG  >  minE/2918687‑2918827
                                                                      |                                                                      
gCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaaaa                                                                        <  1:1594911/71‑1 (MQ=255)
   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:3032873/1‑71 (MQ=255)
         tACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                                               >  1:795772/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:3623585/70‑1 (MQ=255)
              aGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcgg                                                           >  1:1836579/1‑70 (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:2800896/1‑70 (MQ=255)
                                gCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcat                                         >  1:2259504/1‑70 (MQ=255)
                                   ttGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAg                                      >  1:2546262/1‑70 (MQ=255)
                                    tGATAGTATGCATACGGTTTTCCGCCTGATCAAAAAc                                                                      <  1:1141652/37‑1 (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:3396146/70‑1 (MQ=255)
                                                       ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                  <  1:3720328/70‑1 (MQ=255)
                                                        tCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa                 <  1:3683894/70‑1 (MQ=255)
                                                          cGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                               <  1:860264/54‑1 (MQ=255)
                                                          cGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTc                               <  1:2049051/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: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: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)
                                                                      |                                                                      
GCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCG  >  minE/2918687‑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: