Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,408 C→A 6.4% V96V (GTG→GTT argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4080CA 6.4% 265.4 / 2.8 78V96V (GTG→GTTargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (45/28);  new base A (2/3);  total (47/31)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAA  >  minE/2919339‑2919475
                                                                     |                                                                   
ttCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa                                                                    >  1:449873/1‑71 (MQ=255)
ttCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa                                                                    >  1:1606350/1‑71 (MQ=255)
  ccGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                                                  >  1:2587752/1‑71 (MQ=255)
   cGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAcc                                                                  <  1:1021139/70‑2 (MQ=255)
        gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCa                                                            <  1:253905/71‑1 (MQ=255)
         tgtTTCGATAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAg                                                           >  1:2243959/1‑71 (MQ=255)
         tgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAg                                                           >  1:2286494/1‑71 (MQ=255)
            ttCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGgtgt                                                        >  1:1877756/1‑71 (MQ=255)
             tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtgt                                                        <  1:1504363/70‑1 (MQ=255)
                 cAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCttt                                                    <  1:2495971/70‑1 (MQ=255)
                  aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                  <  1:3214295/71‑1 (MQ=255)
                  aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                  >  1:1184423/1‑71 (MQ=255)
                   aTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAAt                                                 >  1:2553387/1‑71 (MQ=255)
                     ctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTa                                                   >  1:3095037/1‑67 (MQ=255)
                     ctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                >  1:4052687/1‑70 (MQ=255)
                     ctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                >  1:597122/1‑70 (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:309011/1‑57 (MQ=255)
                          gACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtgt                                                        >  1:947401/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:1255088/1‑60 (MQ=255)
                                   cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                                            >  1:1069099/1‑60 (MQ=255)
                                   cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                                            >  1:3056088/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)
                                               tgccGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCg                     <  1:482992/69‑1 (MQ=255)
                                               tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                  <  1:927870/42‑1 (MQ=255)
                                               tACCGTCATACATGCGACCAAGAACGCGGGCAGTGTCTTTAATCGa                                              >  1:1483838/1‑46 (MQ=255)
                                                 ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGc                    <  1:3018203/70‑1 (MQ=255)
                                                 ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt                   <  1:1804524/71‑1 (MQ=255)
                                                 ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt                   <  1:2335345/71‑1 (MQ=255)
                                                   gTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCttt                                        <  1:254613/48‑1 (MQ=255)
                                                    tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                              <  1:1763417/57‑1 (MQ=255)
                                                    tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                              <  1:3181627/57‑1 (MQ=255)
                                                    tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                              <  1:414439/57‑1 (MQ=255)
                                                      aTACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                  >  1:1258748/1‑35 (MQ=255)
                                                      aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGcct                     >  1:2081716/1‑63 (MQ=255)
                                                      aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc              >  1:465663/1‑71 (MQ=255)
                                                      aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc              >  1:593064/1‑71 (MQ=255)
                                                      aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc              >  1:3858920/1‑71 (MQ=255)
                                                      aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc              >  1:2231304/1‑71 (MQ=255)
                                                      aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc              >  1:3359434/1‑71 (MQ=255)
                                                       tACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc             >  1:1527442/1‑71 (MQ=255)
                                                         cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                       >  1:2082932/1‑43 (MQ=255)
                                                         cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                       >  1:2087736/1‑43 (MQ=255)
                                                         cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa                                    >  1:3082744/1‑46 (MQ=255)
                                                         cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGc                          >  1:1975978/1‑56 (MQ=255)
                                                         cATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGccct            <  1:1830262/70‑2 (MQ=255)
                                                            gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgag          >  1:175127/1‑69 (MQ=255)
                                                             cGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTg                            >  1:2184810/1‑50 (MQ=255)
                                                               aCCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCa                                 >  1:3903239/1‑43 (MQ=255)
                                                                 cAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt    <  1:3181648/70‑1 (MQ=255)
                                                                  aaGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa  >  1:2889833/1‑71 (MQ=255)
                                                                   aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCga           <  1:1348478/61‑1 (MQ=255)
                                                                     |                                                                   
TTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAA  >  minE/2919339‑2919475

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: