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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,3870AG 6.6% 260.6 / 8.3 91G103G (GGT→GGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (50/35);  new base G (4/2);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

ATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGC  >  minE/2919320‑2919454
                                                                   |                                                                   
aTACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCg                                                                  <  1:547756/71‑1 (MQ=255)
 tACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATAc                                                                    <  1:3470373/68‑1 (MQ=255)
 tACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATAGCGt                                                                 <  1:1016432/71‑1 (MQ=255)
 caCCGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGt                                                                 <  1:3084572/70‑1 (MQ=255)
  aCCGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGt                                                                 >  1:257623/1‑70 (MQ=255)
   ccGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                               >  1:2365229/1‑71 (MQ=255)
   ccGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                               >  1:2075023/1‑71 (MQ=255)
   ccGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                               >  1:3326762/1‑71 (MQ=255)
    cGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCa                                                               >  1:2466837/1‑70 (MQ=255)
    cGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCa                                                               >  1:651248/1‑70 (MQ=255)
      gCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATAc                                                            <  1:2119899/71‑1 (MQ=255)
      gCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATAc                                                            <  1:3357700/71‑1 (MQ=255)
          gCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                               <  1:4011572/64‑1 (MQ=255)
           cTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCg                                                       <  1:1613966/71‑1 (MQ=255)
              gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAc                                                     >  1:3695835/1‑70 (MQ=255)
               cGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                    >  1:1433522/1‑70 (MQ=255)
               cGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                    <  1:2830463/70‑1 (MQ=255)
                 tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGt                                                                 <  1:2525052/55‑1 (MQ=255)
                 tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                    <  1:168086/68‑1 (MQ=255)
                 tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAg                                                 <  1:1179099/71‑1 (MQ=255)
                  aTTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAg                                                 >  1:1365459/1‑70 (MQ=255)
                   ttCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa                                               >  1:449873/1‑71 (MQ=255)
                   ttCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa                                               >  1:1606350/1‑71 (MQ=255)
                     ccGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                             >  1:2587752/1‑71 (MQ=255)
                      cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                    <  1:1292680/63‑1 (MQ=255)
                      cGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAcc                                             <  1:1021139/70‑2 (MQ=255)
                        ccAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCAt                                                              <  1:1134790/51‑1 (MQ=255)
                          aGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCg                                                       >  1:2856499/1‑56 (MQ=255)
                          aGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                      >  1:3363395/1‑57 (MQ=255)
                           gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCa                                       <  1:253905/71‑1 (MQ=255)
                            tgtTTCGATAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAg                                      >  1:2243959/1‑71 (MQ=255)
                            tgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAg                                      >  1:2286494/1‑71 (MQ=255)
                              tttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATa                                                             <  1:2094733/46‑1 (MQ=255)
                              tttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATa                                                             >  1:1494379/1‑46 (MQ=255)
                               ttCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGgtgt                                   >  1:1877756/1‑71 (MQ=255)
                                tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGt                                                                 >  1:578407/1‑40 (MQ=255)
                                tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGt                                                                 >  1:811116/1‑40 (MQ=255)
                                tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                               >  1:2151629/1‑42 (MQ=255)
                                tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                               >  1:3610495/1‑42 (MQ=255)
                                tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                    >  1:2718988/1‑53 (MQ=255)
                                tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtgt                                   <  1:1504363/70‑1 (MQ=255)
                                  gACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATa                                                             <  1:3222743/42‑1 (MQ=255)
                                   aCAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATa                                                             <  1:796291/41‑1 (MQ=255)
                                    cAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACa                                                           >  1:1724808/1‑42 (MQ=255)
                                    cAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCttt                               <  1:2495971/70‑1 (MQ=255)
                                     aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCaa                                                  >  1:495885/1‑50 (MQ=255)
                                     aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCaa                                                  >  1:3081558/1‑50 (MQ=255)
                                     aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCaa                                                  >  1:1678301/1‑50 (MQ=255)
                                     aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                             >  1:1184423/1‑71 (MQ=255)
                                     aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                             <  1:3214295/71‑1 (MQ=255)
                                      aTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAAt                            >  1:2553387/1‑71 (MQ=255)
                                        ctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTa                              >  1:3095037/1‑67 (MQ=255)
                                        ctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                           >  1:597122/1‑70 (MQ=255)
                                        ctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                           >  1:4052687/1‑70 (MQ=255)
                                         tcCTGACCATAGCCGCGATACTGAATACCGTCATACa                                                           >  1:373458/1‑37 (MQ=255)
                                         tcCTGACCATAGCCGCGATACTGAATACCGTCATACa                                                           >  1:240459/1‑37 (MQ=255)
                                           cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                           <  1:685844/67‑1 (MQ=255)
                                           cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                         >  1:2786790/1‑69 (MQ=255)
                                           cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                         <  1:1358783/69‑1 (MQ=255)
                                             gACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa                                               >  1:2257563/1‑45 (MQ=255)
                                             gACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtgt                                   >  1:947401/1‑57 (MQ=255)
                                             gACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtgt                                   >  1:309011/1‑57 (MQ=255)
                                             gACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActct                     >  1:1004185/1‑71 (MQ=255)
                                             gACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActct                     >  1:732057/1‑71 (MQ=255)
                                              aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActc                      <  1:209558/69‑1 (MQ=255)
                                                cATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActct                     <  1:1241767/68‑1 (MQ=255)
                                                 aTAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGGGGGCAGTGTCTTTAATCGACTCTTTAt                 <  1:3481903/71‑1 (MQ=255)
                                                 aTAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTa                              <  1:2810847/58‑1 (MQ=255)
                                                  tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActc                      <  1:3852490/65‑1 (MQ=255)
                                                  tAGCCGCGATACTGAATACCGTCATACATGCGACAAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                 >  1:355954/1‑70 (MQ=255)
                                                   aGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa               >  1:1128770/1‑71 (MQ=255)
                                                    gCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc              >  1:2593185/1‑71 (MQ=255)
                                                     ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCttt                               <  1:1833507/53‑1 (MQ=255)
                                                      cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCt                                 >  1:3722784/1‑50 (MQ=255)
                                                      cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                       >  1:3056088/1‑60 (MQ=255)
                                                      cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                       >  1:1069099/1‑60 (MQ=255)
                                                      cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                       >  1:1255088/1‑60 (MQ=255)
                                                      cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCa            <  1:1120459/71‑1 (MQ=255)
                                                          aTACTGAATACCGTCATACATGCGACCAAGCACGCGGGc                                        <  1:3246968/39‑1 (MQ=255)
                                                           tACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                    <  1:4109202/58‑1 (MQ=255)
                                                            aatGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                           >  1:2421525/3‑50 (MQ=255)
                                                            aCTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATc         >  1:1156875/1‑68 (MQ=255)
                                                             cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg      <  1:413946/70‑1 (MQ=255)
                                                              tGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTa                              <  1:627498/45‑1 (MQ=255)
                                                               gAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                             >  1:1442746/1‑45 (MQ=255)
                                                               gAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg   >  1:1953984/1‑71 (MQ=255)
                                                                aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                         >  1:1581934/1‑48 (MQ=255)
                                                                aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGc  >  1:3935998/1‑71 (MQ=255)
                                                                 aTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                           <  1:4037096/45‑1 (MQ=255)
                                                                  tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                             <  1:927870/42‑1 (MQ=255)
                                                                  tACCGTCATACATGCGACCAAGAACGCGGGCAGTGTCTTTAATCGa                         >  1:1483838/1‑46 (MQ=255)
                                                                   |                                                                   
ATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGC  >  minE/2919320‑2919454

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: