Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,480 A→C 16.5% A72A (GCT→GCG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4800AC16.5% 159.7 / 35.9 85A72A (GCT→GCGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base C (7/7);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAG  >  minE/2919410‑2919548
                                                                      |                                                                    
cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                                      <  1:1212140/71‑1 (MQ=255)
 gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                                                     >  1:2676850/1‑71 (MQ=255)
 gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                                                     >  1:3674753/1‑71 (MQ=255)
      cAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCAccc                                                                >  1:1171111/1‑71 (MQ=255)
        gtgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGccc                                                                <  1:1796656/69‑1 (MQ=255)
         tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                                             >  1:305225/1‑71 (MQ=255)
         tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                                             >  1:3772954/1‑71 (MQ=255)
          gtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                                             <  1:1451349/70‑1 (MQ=255)
          gtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGt                                                            >  1:2564401/1‑71 (MQ=255)
           tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                                      >  1:3227485/1‑60 (MQ=255)
           tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                                      >  1:140120/1‑60 (MQ=255)
           tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                                      >  1:135462/1‑60 (MQ=255)
             tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                                                     >  1:456162/1‑59 (MQ=255)
             tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCa                                                          >  1:3480390/1‑70 (MQ=255)
             tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCa                                                          >  1:910340/1‑70 (MQ=255)
               tAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCaaat                                                       >  1:3334071/1‑71 (MQ=255)
               tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCatat                                                       >  1:1032572/1‑71 (MQ=255)
               tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAAGCGAGCACCCTGGTCatat                                                       >  1:3624586/1‑71 (MQ=255)
                aaTTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCt                                                               >  1:2810014/1‑62 (MQ=255)
                  tCGACTCTTTATGACCAATCTGGCTGACGCTTGGGCCGAGATAAGTAACGCGAGCACCCt                                                               >  1:1648351/1‑60 (MQ=255)
                     aCTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCa                                                 <  1:1679702/71‑1 (MQ=255)
                      ctctTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCaa                                                >  1:3308221/1‑71 (MQ=255)
                       tctTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGc                                                  >  1:1662501/1‑68 (MQ=255)
                       tctTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCaa                                                >  1:2704415/1‑70 (MQ=255)
                       tctTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAAc                                               >  1:788008/1‑71 (MQ=255)
                         tttATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAAc                                               <  1:962910/69‑1 (MQ=255)
                         tttATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCATGGTCATATGCGGCAAc                                               <  1:1954740/69‑1 (MQ=255)
                          ttATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGc                                                  <  1:2668362/65‑1 (MQ=255)
                           tATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTc                                            <  1:339524/70‑1 (MQ=255)
                           tATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTc                                            >  1:575732/1‑70 (MQ=255)
                            aTGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCg                                           <  1:417330/70‑1 (MQ=255)
                            aTGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCg                                           >  1:3808433/1‑70 (MQ=255)
                             tGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                         <  1:1746299/71‑1 (MQ=255)
                             tGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGa                                          <  1:3506041/70‑1 (MQ=255)
                             tGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGa                                          <  1:773765/70‑1 (MQ=255)
                              gACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGc                                                     <  1:503415/58‑1 (MQ=255)
                              gACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGAACCCTTGACATATGc                                                     <  1:1341441/58‑1 (MQ=255)
                               cccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGaaa                                        <  1:3064838/69‑1 (MQ=255)
                               cccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAag                                       <  1:3312819/70‑1 (MQ=255)
                               aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                                      >  1:965550/1‑40 (MQ=255)
                               aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaaa                                        <  1:244577/70‑1 (MQ=255)
                               aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAag                                       <  1:1603931/71‑1 (MQ=255)
                               aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAag                                       <  1:2414506/71‑1 (MQ=255)
                                ccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGaaa                                        >  1:1550738/1‑69 (MQ=255)
                                ccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGaaa                                        >  1:3400372/1‑69 (MQ=255)
                                 cAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCCCCCTGGTCATATGCGGCAACTTCGAAagag                                     >  1:2090553/1‑71 (MQ=255)
                                 cAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAagag                                     >  1:2074159/1‑71 (MQ=255)
                                   aTCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCa                                   <  1:3100397/71‑1 (MQ=255)
                                    tCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCAt                                  >  1:1266996/1‑71 (MQ=255)
                                     cTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCa                                                                   >  1:2226734/1‑37 (MQ=255)
                                      tGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGc                                                     >  1:2709100/1‑50 (MQ=255)
                                      tGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCg                                >  1:1136939/1‑71 (MQ=255)
                                        gCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATc                                 <  1:988330/68‑1 (MQ=255)
                                         cTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGt                             <  1:1150435/71‑1 (MQ=255)
                                            ccGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                           >  1:1021955/1‑70 (MQ=255)
                                            ccGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTAc                           >  1:1352634/1‑70 (MQ=255)
                                             cGCTTGGGCCGAGATAAGTAACGCGAGCACGCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                           <  1:1755589/69‑1 (MQ=255)
                                             cGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcga                         <  1:47219/71‑1 (MQ=255)
                                             cGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAaga                                      >  1:1256399/1‑58 (MQ=255)
                                              gCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgag                        <  1:3292804/71‑1 (MQ=255)
                                               cTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATg                                                      <  1:2390389/40‑1 (MQ=255)
                                               cTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATc                                 <  1:3930069/61‑1 (MQ=255)
                                               cTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagt                       >  1:175560/1‑71 (MQ=255)
                                                 tGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                      >  1:764412/1‑70 (MQ=255)
                                                  gggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcg                     <  1:2916072/70‑1 (MQ=255)
                                                  gggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                      <  1:733628/69‑1 (MQ=255)
                                                  gggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                      <  1:3227791/69‑1 (MQ=255)
                                                   ggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcgag                   <  1:509753/71‑1 (MQ=255)
                                                     ccGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACtt                                             >  1:3289058/1‑43 (MQ=255)
                                                      cGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                         <  1:2694615/46‑1 (MQ=255)
                                                      cGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                         <  1:519303/46‑1 (MQ=255)
                                                      cGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCt                >  1:676457/1‑71 (MQ=255)
                                                       gagaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGt                             >  1:3365593/1‑57 (MQ=255)
                                                         gaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCtttt             >  1:2066507/1‑71 (MQ=255)
                                                          aTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACt                                              <  1:1050373/37‑1 (MQ=255)
                                                          aTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCttttt            <  1:1716666/71‑1 (MQ=255)
                                                            aaGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCa                                   <  1:2993297/46‑1 (MQ=255)
                                                              gTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa        <  1:2915227/71‑1 (MQ=255)
                                                               tAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                         <  1:2984119/37‑1 (MQ=255)
                                                                aaCGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTAcgagtcg                     <  1:1874776/56‑1 (MQ=255)
                                                                 aCGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGAt     <  1:3790563/71‑1 (MQ=255)
                                                                 aCGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                         <  1:3311641/35‑1 (MQ=255)
                                                                 aCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTc           <  1:635119/65‑1 (MQ=255)
                                                                 aCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGAt     <  1:902519/71‑1 (MQ=255)
                                                                 aCGCGAGCACCCTGGACATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGAt     <  1:3738397/71‑1 (MQ=255)
                                                                  cgcgcgcgCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCGTTTTCGAAGAt     >  1:1109394/1‑70 (MQ=255)
                                                                  cgcgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGt                             <  1:278775/46‑1 (MQ=255)
                                                                   gcgcgcgcCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGAt     <  1:590861/69‑1 (MQ=255)
                                                                   gcgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGAt     >  1:1052411/1‑69 (MQ=255)
                                                                   gcgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGa   <  1:360040/71‑1 (MQ=255)
                                                                    cgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAg  <  1:2008898/71‑1 (MQ=255)
                                                                     gAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAg       >  1:1141367/1‑65 (MQ=255)
                                                                      |                                                                    
CGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAG  >  minE/2919410‑2919548

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: