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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,0210AG 5.2% 266.2 / 5.5 96G225G (GGT→GGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (45/45);  new base G (3/2);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

CAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGG  >  minE/2918951‑2919090
                                                                      |                                                                     
caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCa                                                                       >  1:1749111/1‑71 (MQ=255)
caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCa                                                                       >  1:788911/1‑71 (MQ=255)
 aatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAc                                                                      <  1:1655502/71‑1 (MQ=255)
  atCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAcc                                                                     >  1:1118060/1‑71 (MQ=255)
   tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCt                                                                    >  1:1191696/1‑71 (MQ=255)
    ccGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGccct                                                                   >  1:3953663/1‑69 (MQ=255)
     cGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTc                                                                  <  1:3473949/71‑1 (MQ=255)
       ttCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCaa                                                                <  1:491313/71‑1 (MQ=255)
       ttCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCaa                                                                <  1:2543313/71‑1 (MQ=255)
             ccATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                                                          >  1:4072501/1‑71 (MQ=255)
             ccATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                                                          >  1:4035780/1‑71 (MQ=255)
              cATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCtt                                                         >  1:207836/1‑71 (MQ=255)
              cATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCtt                                                         >  1:3879067/1‑71 (MQ=255)
               aTTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTc                                                        >  1:1376554/1‑71 (MQ=255)
                tttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCa                                                                 <  1:2886486/61‑1 (MQ=255)
                 ttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc                                                      >  1:55714/1‑71 (MQ=255)
                 ttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc                                                      >  1:1179593/1‑71 (MQ=255)
                 ttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc                                                      >  1:3351464/1‑71 (MQ=255)
                 ttCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAAcg                                                              <  1:3108798/63‑2 (MQ=255)
                  tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc                                                      <  1:389013/70‑1 (MQ=255)
                  tctcTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGccct                                                                   >  1:1783061/1‑55 (MQ=255)
                   ctctTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgcg                                                     <  1:3346098/70‑1 (MQ=255)
                   ctctTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgcg                                                     <  1:2465857/70‑1 (MQ=255)
                    tctTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGa                                                    >  1:912239/1‑70 (MQ=255)
                      ttttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGa                                                    >  1:1207318/1‑68 (MQ=255)
                        ttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCt                                               >  1:1708280/1‑71 (MQ=255)
                        ttGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCt                                               >  1:906231/1‑71 (MQ=255)
                         tGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCtt                                              >  1:2821745/1‑71 (MQ=255)
                              cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCAcc                                                                     >  1:2565414/1‑43 (MQ=255)
                              cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                          >  1:187321/1‑70 (MQ=255)
                              cccccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                          >  1:232742/1‑70 (MQ=255)
                                 ccATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTcc                                                            >  1:892692/1‑49 (MQ=255)
                                   aTCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGt                                    <  1:3874861/71‑1 (MQ=255)
                                     cGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  >  1:1279946/1‑71 (MQ=255)
                                     cGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  >  1:1888978/1‑71 (MQ=255)
                                      gACACCGACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  >  1:533363/1‑70 (MQ=255)
                                      gACACCCACACATCGGTATAGATATAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  <  1:2831453/70‑1 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGa                                                    >  1:2628101/1‑52 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGAc                                                   <  1:950714/53‑1 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGAc                                                   <  1:1431769/53‑1 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCg                                     >  1:1866518/1‑67 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  >  1:567329/1‑70 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  >  1:1843052/1‑70 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  <  1:2329144/70‑1 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  <  1:25592/70‑1 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  >  1:1419252/1‑70 (MQ=255)
                                      gACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAat                                 <  1:3871625/71‑1 (MQ=255)
                                      gACACACACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGa                                                    >  1:4143229/1‑52 (MQ=255)
                                        cacCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCTACATCTTCAGTCAGCGTAat                                 >  1:1353509/1‑69 (MQ=255)
                                        cacCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGa                                                    >  1:2165919/1‑50 (MQ=255)
                                        cacCCACACATCGGTATAGATAAAGTCAGCACATTc                                                                  >  1:525436/1‑36 (MQ=255)
                                         acCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                              <  1:3445834/71‑1 (MQ=255)
                                         acCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATAtt                              <  1:2641345/71‑1 (MQ=255)
                                          cccACACATCGGTATAGATAAAGTCAGCACCTTCa                                                                 <  1:819463/35‑1 (MQ=255)
                                          cccACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                                                          <  1:1630634/42‑1 (MQ=255)
                                          cccACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGa                                                    <  1:2442358/48‑1 (MQ=255)
                                            cacacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTc                                                        >  1:1366151/1‑42 (MQ=255)
                                            cacacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGAc                                                   >  1:1417169/1‑47 (MQ=255)
                                            cacacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGt                                    <  1:3474849/62‑1 (MQ=255)
                                            cacacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTcc                            >  1:814860/1‑70 (MQ=255)
                                             acacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTa                                   <  1:3881004/62‑1 (MQ=255)
                                              cacgTCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAATTTTcc                            >  1:2494269/1‑68 (MQ=255)
                                              cacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtc                                             <  1:625418/51‑1 (MQ=255)
                                              cacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccac                         >  1:2873968/1‑71 (MQ=255)
                                              cacaTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccac                         >  1:2303700/1‑71 (MQ=255)
                                               acaTCGGTATAGATAAAGTCAGCCCCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCcca                          >  1:3146804/1‑69 (MQ=255)
                                               acaTCGGTATAGATAAAGTCAGCACCTTCAACTccc                                                           <  1:277066/36‑1 (MQ=255)
                                                caTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAata                                >  1:1440053/1‑62 (MQ=255)
                                                caTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCccacca                       <  1:3201278/71‑1 (MQ=255)
                                                 gtCGGTATAGATAAAGTCCGCGCCCTTAACGCCTGCCGCCACGTCTTCCGTCAGAGTAATTTTCccgccg                       <  1:2594638/69‑2 (MQ=255)
                                                 aTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAt                      <  1:117986/71‑1 (MQ=255)
                                                 aTCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAt                      <  1:2323555/71‑1 (MQ=255)
                                                  tCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgcg                                                     <  1:435819/39‑1 (MQ=255)
                                                   cGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagtc                                         >  1:1639524/1‑50 (MQ=255)
                                                    ggTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCAtttt                   <  1:3512839/71‑1 (MQ=255)
                                                      tataGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCg                                     <  1:3975790/51‑1 (MQ=255)
                                                       ataGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagtca                                        <  1:2388900/47‑1 (MQ=255)
                                                       ataGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCt                >  1:3518149/1‑71 (MQ=255)
                                                          gATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTcc                            >  1:1268009/1‑56 (MQ=255)
                                                          gATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGACAGCGTAATATTCCCACCATTTTGCtgt              <  1:286310/70‑1 (MQ=255)
                                                           aTAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  <  1:1023977/49‑1 (MQ=255)
                                                           aTAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  <  1:851402/49‑1 (MQ=255)
                                                            tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTtcagt                                          >  1:1902287/1‑40 (MQ=255)
                                                            tAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGcc           <  1:3676623/71‑1 (MQ=255)
                                                             aaaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGc            >  1:1057399/1‑69 (MQ=255)
                                                              aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGt                                    <  1:2252295/44‑1 (MQ=255)
                                                              aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCa          <  1:601629/70‑1 (MQ=255)
                                                              aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCa          <  1:691638/70‑1 (MQ=255)
                                                              aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAg         >  1:4007706/1‑71 (MQ=255)
                                                              aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAg         >  1:1264572/1‑71 (MQ=255)
                                                              aaGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAg         <  1:1350344/71‑1 (MQ=255)
                                                                 tCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTaa                                  <  1:3160720/43‑1 (MQ=255)
                                                                 tCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGcc           <  1:1180054/66‑1 (MQ=255)
                                                                  cAGCACCTTCAACTCCCTTCGCGACATGTTCAGTCAGCGTAATATTCCCACCATTTTGCtg               >  1:481170/1‑61 (MQ=255)
                                                                     cACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCgg  <  1:540574/71‑1 (MQ=255)
                                                                     cACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCgg  <  1:1036299/71‑1 (MQ=255)
                                                                     cACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCgg  <  1:3839254/71‑1 (MQ=255)
                                                                     cACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTCCTGTGCCAGGGCGCgg  <  1:3144181/71‑1 (MQ=255)
                                                                      |                                                                     
CAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGG  >  minE/2918951‑2919090

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: