Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,423 T→C 8.9% K91K (AAA→AAG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4230TC 8.9% 252.4 / 12.5 90K91K (AAA→AAGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (45/37);  new base C (4/4);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 101 total aligned reads displayed.

ACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCA  >  minE/2919355‑2919492
                                                                    |                                                                     
aCAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                                                                      >  1:3666389/1‑70 (MQ=255)
  aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                                   >  1:3046045/1‑71 (MQ=255)
  aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                                   <  1:312579/71‑1 (MQ=255)
  aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                                                                   <  1:492581/71‑1 (MQ=255)
    tctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                                 >  1:2754572/1‑71 (MQ=255)
      tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                                               >  1:3655131/1‑71 (MQ=255)
      tcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                                               >  1:1711836/1‑71 (MQ=255)
       ccTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                                              >  1:1750873/1‑71 (MQ=255)
        cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                                 <  1:2676158/67‑1 (MQ=255)
        cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCg                                                                <  1:3462241/68‑1 (MQ=255)
        cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                                                             >  1:2337930/1‑71 (MQ=255)
        cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                                              >  1:1923432/1‑70 (MQ=255)
         tGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCt                                                                       <  1:1139500/60‑1 (MQ=255)
           aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCt                                                                       >  1:1550077/1‑58 (MQ=255)
           aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActc                                                            <  1:342752/69‑1 (MQ=255)
           aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                                          >  1:1105342/1‑71 (MQ=255)
           aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                                          >  1:3646189/1‑71 (MQ=255)
            ccATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCttt                                                         <  1:1678526/71‑1 (MQ=255)
             cATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                                                                      <  1:282461/57‑1 (MQ=255)
              aTAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                                       <  1:1988056/71‑1 (MQ=255)
               tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                                       >  1:1231583/1‑70 (MQ=255)
               tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                                       <  1:3021845/70‑1 (MQ=255)
               tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                                      <  1:1195531/71‑1 (MQ=255)
               tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                                      <  1:3527365/71‑1 (MQ=255)
               tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTCTCTTTAATCGACTCTTTAt                                                       >  1:2264886/1‑70 (MQ=255)
                aGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTa                                                                    <  1:2481641/56‑1 (MQ=255)
                  ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa                                                     >  1:2516378/1‑69 (MQ=255)
                  ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                                                    >  1:3877993/1‑70 (MQ=255)
                  ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                                                    >  1:287660/1‑70 (MQ=255)
                   cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAcc                                                   >  1:3140073/1‑70 (MQ=255)
                     cgATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAAt                                                >  1:3343308/1‑71 (MQ=255)
                     cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                                                <  1:1501430/71‑1 (MQ=255)
                     cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                                                <  1:1072865/71‑1 (MQ=255)
                       aTACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt                                              <  1:3140774/71‑1 (MQ=255)
                        tACTTAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt                                              <  1:3182861/70‑1 (MQ=255)
                        tACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTg                                             >  1:1151001/1‑71 (MQ=255)
                            gAAAACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAc                                                              <  1:44526/50‑1 (MQ=255)
                             aaaaCCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                         <  1:2421608/67‑1 (MQ=255)
                             aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg                                            >  1:1196747/1‑67 (MQ=255)
                             aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                         <  1:3400262/70‑1 (MQ=255)
                             aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                         <  1:1588042/70‑1 (MQ=255)
                             aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                         <  1:3576235/70‑1 (MQ=255)
                              aTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                                       >  1:3646598/1‑55 (MQ=255)
                               tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATg                                                      <  1:1969164/55‑1 (MQ=255)
                               tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                                      >  1:2875305/1‑71 (MQ=255)
                               tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                                      <  1:1014431/71‑1 (MQ=255)
                                 ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt                                    <  1:530444/71‑1 (MQ=255)
                                  cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                                   <  1:3216134/71‑1 (MQ=255)
                                  cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                                   >  1:390611/1‑71 (MQ=255)
                                    tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTgg                                 <  1:3916074/71‑1 (MQ=255)
                                    tCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCgg                                 >  1:3047201/1‑71 (MQ=255)
                                       tACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                              <  1:3711769/71‑1 (MQ=255)
                                         cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                                          <  1:1304227/41‑1 (MQ=255)
                                         cATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGcccta                            <  1:1101366/71‑3 (MQ=255)
                                          aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                              >  1:3325441/1‑68 (MQ=255)
                                          aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                              >  1:2753508/1‑68 (MQ=255)
                                            gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                                        >  1:64676/1‑40 (MQ=255)
                                            gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                          >  1:286668/1‑70 (MQ=255)
                                            gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                          >  1:2474233/1‑70 (MQ=255)
                                            gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                          >  1:1671216/1‑70 (MQ=255)
                                             cGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGAt                         >  1:628717/1‑70 (MQ=255)
                                              gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                                       <  1:3747826/39‑1 (MQ=255)
                                              gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                       <  1:3385590/71‑1 (MQ=255)
                                              gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                       <  1:1545066/71‑1 (MQ=255)
                                                ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGcc                                       <  1:1916688/53‑1 (MQ=255)
                                                ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                      >  1:3744153/1‑70 (MQ=255)
                                                ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                      >  1:1435759/1‑70 (MQ=255)
                                                ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                      >  1:2606645/1‑70 (MQ=255)
                                                ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                      >  1:3253438/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)
                                                       cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                                               <  1:3360299/38‑1 (MQ=255)
                                                       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)
                                                              ggtgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCt                                    <  1:1502732/41‑1 (MQ=37)
                                                               gtgtCTTTAATCGACTCTTTATGACCAATCTGGCTGcc                                       >  1:80866/1‑38 (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:140120/1‑60 (MQ=255)
                                                                  tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa              >  1:3227485/1‑60 (MQ=255)
                                                                  tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa              >  1:135462/1‑60 (MQ=255)
                                                                    tttAATCGACTCTTTATGACCAATCTGGCTGCCGCt                                    <  1:1706149/36‑1 (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)
                                                                    |                                                                     
ACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCA  >  minE/2919355‑2919492

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: