Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,483 A→G 16.8% G71G (GGT→GGC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4830AG16.8% 183.9 / 29.5 84G71G (GGT→GGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (27/42);  new base G (7/7);  total (35/49)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

CAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCG  >  minE/2919416‑2919552
                                                                   |                                                                     
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)
       tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCa                                                              >  1:910340/1‑70 (MQ=255)
       tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCa                                                              >  1:3480390/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:3808433/1‑70 (MQ=255)
                      aTGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCg                                               <  1:417330/70‑1 (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)
                         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)
                           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)
                                                                  cACCCTGGTCATATGCTGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGcgcg  <  1:2405417/71‑1 (MQ=255)
                                                                  cACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                               <  1:764474/42‑1 (MQ=255)
                                                                  cACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAgc     <  1:2483171/68‑1 (MQ=255)
                                                                   |                                                                     
CAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCG  >  minE/2919416‑2919552

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: