Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,736 A→C 15.5% R320R (CGT→CGG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,7360AC15.5% 112.0 / 21.8 58R320R (CGT→CGGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (27/22);  new base C (5/4);  total (32/26)
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.59e-01

CCGGCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGC  >  minE/2918667‑2918806
                                                                     |                                                                      
ccGGCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATa                                                                        <  1:3155333/70‑1 (MQ=255)
  ggCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACg                                                                      >  1:1256403/1‑70 (MQ=255)
     aTCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGttt                                                                  <  1:2928845/71‑1 (MQ=255)
     aTCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGttt                                                                  <  1:2762301/71‑1 (MQ=255)
       cTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTc                                                                >  1:623347/1‑71 (MQ=255)
        tttATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTcc                                                               >  1:1111846/1‑71 (MQ=255)
        tttATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTcc                                                               >  1:879681/1‑71 (MQ=255)
            taGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCt                                                           <  1:537347/71‑1 (MQ=255)
             aGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCt                                                           >  1:516079/1‑70 (MQ=255)
                   aGCAATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaa                                                     >  1:2296540/1‑70 (MQ=255)
                   aGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaa                                                     >  1:368211/1‑70 (MQ=255)
                   aGCAATTATTTACTGAGCGTCGAGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaa                                                     >  1:1151857/1‑70 (MQ=255)
                    gCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaaaa                                                   <  1:1594911/71‑1 (MQ=255)
                       attattTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACa                                                <  1:1452984/71‑1 (MQ=255)
                       attattTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACa                                                <  1:2201761/71‑1 (MQ=255)
                         tattTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGc                                                             >  1:210804/1‑56 (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)
                              aCTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCt                                                           >  1:2843328/1‑53 (MQ=255)
                              aCTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                                          <  1:609998/70‑1 (MQ=255)
                                tGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaa                                                     <  1:2005096/57‑1 (MQ=255)
                                tGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                                       >  1:1903999/1‑71 (MQ=255)
                                caaGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgcc                                       >  1:527391/3‑70 (MQ=255)
                                 aaGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgccg                                      >  1:1215440/2‑69 (MQ=255)
                                  aGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTcg                                                      <  1:1579608/54‑2 (MQ=255)
                                  aGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcgg                                      >  1:1836579/1‑70 (MQ=255)
                                  aGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcgg                                      <  1:3623585/70‑1 (MQ=255)
                                   gCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTg                                                          <  1:1589865/49‑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)
                                                   cGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTcg                                                      <  1:3110759/37‑2 (MQ=16)
                                                    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)
                                                                     aCGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa                          >  1:627506/1‑47 (MQ=255)
                                                                     aCGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc  >  1:3433599/1‑71 (MQ=255)
                                                                     |                                                                      
CCGGCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGC  >  minE/2918667‑2918806

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: