Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,465 G→T 14.5% L77L (CTC→CTA) ‡ argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4650GT14.5% 278.7 / 22.4 96L77L (CTC→CTA) ‡argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (43/39);  new base T (8/6);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 102 total aligned reads displayed.

GCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCT  >  minE/2919399‑2919534
                                                                  |                                                                     
gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                                                                    >  1:1671216/1‑70 (MQ=255)
gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                                                                    >  1:2474233/1‑70 (MQ=255)
gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                                                                    >  1:286668/1‑70 (MQ=255)
 cGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGAt                                                                   >  1:628717/1‑70 (MQ=255)
  gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                                                                 <  1:3385590/71‑1 (MQ=255)
  gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                                                                 <  1:1545066/71‑1 (MQ=255)
    ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                                                >  1:2606645/1‑70 (MQ=255)
    ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                                                >  1:3253438/1‑70 (MQ=255)
    ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                                                >  1:3744153/1‑70 (MQ=255)
    ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                                                >  1:1435759/1‑70 (MQ=255)
     cAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTa                                                              <  1:1338284/71‑1 (MQ=255)
     cAAGCACGCGGGCAGTGTCTTTAATAGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt                                                               >  1:3808790/1‑70 (MQ=255)
       aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc                                                            >  1:2019871/1‑71 (MQ=255)
       aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc                                                            <  1:800302/71‑1 (MQ=255)
         cACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                                                          <  1:198903/71‑1 (MQ=255)
         cACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcg                                                           >  1:1024561/1‑70 (MQ=255)
           cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                        <  1:1212140/71‑1 (MQ=255)
            gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                                       >  1:3674753/1‑71 (MQ=255)
            gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                                       >  1:2676850/1‑71 (MQ=255)
                 cAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCAccc                                                  >  1:1171111/1‑71 (MQ=255)
                   gtgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGccc                                                  <  1:1796656/69‑1 (MQ=255)
                    tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa                                                             >  1:2262576/1‑57 (MQ=255)
                    tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                               >  1:3772954/1‑71 (MQ=255)
                    tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                               >  1:305225/1‑71 (MQ=255)
                    tgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAAc                                                            >  1:1536508/1‑58 (MQ=255)
                     gtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                               <  1:1451349/70‑1 (MQ=255)
                     gtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGt                                              >  1:2564401/1‑71 (MQ=255)
                     gtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAAcgc                                                          <  1:1223077/59‑1 (MQ=255)
                      tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCga                                                                      >  1:1111400/1‑46 (MQ=255)
                      tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                        >  1:135462/1‑60 (MQ=255)
                      tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                        >  1:3227485/1‑60 (MQ=255)
                      tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                        >  1:140120/1‑60 (MQ=255)
                        tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt                                                               <  1:3029960/51‑1 (MQ=255)
                        tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                                                          <  1:2679631/56‑1 (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)
                            aTCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc                                                            <  1:3024510/50‑1 (MQ=255)
                             tCGACTCTTTATGACCAATCTGGCTGACGCTTGGGCCGAGATAAGTAACGCGAGCACCCt                                                 >  1:1648351/1‑60 (MQ=255)
                               gACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgag                                                                     <  1:1462822/38‑1 (MQ=255)
                                aCTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGAt                                                                   <  1:889712/39‑1 (MQ=255)
                                aCTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa                                                             <  1:1961293/45‑1 (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)
                                       aTGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa                                                             >  1:2501674/1‑38 (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:773765/70‑1 (MQ=255)
                                        tGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGa                            <  1:3506041/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)
                                          aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                                                          <  1:1027469/38‑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:3400372/1‑69 (MQ=255)
                                           ccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGaaa                          >  1:1550738/1‑69 (MQ=255)
                                            cAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                                                          >  1:657330/1‑36 (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)
                                                    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:519303/46‑1 (MQ=255)
                                                                 cGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCt  >  1:676457/1‑71 (MQ=255)
                                                                  gagaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGt               >  1:3365593/1‑57 (MQ=255)
                                                                  |                                                                     
GCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCT  >  minE/2919399‑2919534

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: