Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,982 C→G 8.4% G238G (GGG→GGC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,9820CG 8.4% 246.4 / 9.8 83G238G (GGG→GGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (46/30);  new base G (5/2);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GCATCATCTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGT  >  minE/2918912‑2919050
                                                                      |                                                                    
gCATCATCTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTcc                                                                      >  1:2768666/1‑71 (MQ=255)
 catcatCTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTccc                                                                     >  1:585629/1‑71 (MQ=255)
     atcTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTccccatc                                                                 <  1:3901372/71‑4 (MQ=255)
     atcTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCa                                                                  >  1:2545675/1‑70 (MQ=255)
        ttGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGa                                                              >  1:889573/1‑71 (MQ=255)
        ttGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGa                                                              >  1:2551267/1‑71 (MQ=255)
        ttGCTGTTCACCTGATATTCACGCAGCAATGCAATACGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGa                                                              >  1:73977/1‑71 (MQ=255)
            tGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACAcc                                                          >  1:249702/1‑71 (MQ=255)
            cgTTCACCTGATACCCGCGCAGCAGCGCAATCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACAcc                                                          <  1:2222881/70‑1 (MQ=255)
              ttCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCca                                                        <  1:2847787/71‑1 (MQ=255)
               tCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCAt                                                                 <  1:2018299/61‑1 (MQ=255)
                cACCTGATACCCGCGCAGCAGCGCAATCCGCTCTGCCCACTTCTCTTTGGCCTCGCCCATCGACACCcaca                                                      >  1:605667/1‑71 (MQ=255)
                 aCCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCcaca                                                      >  1:3039307/1‑70 (MQ=255)
                  ccTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATc                                                                <  1:1095336/59‑1 (MQ=255)
                   cTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTccc                                                                     >  1:474834/1‑53 (MQ=255)
                    tGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCa                                                                  >  1:3618576/1‑55 (MQ=255)
                    tGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCa                                                                  >  1:1720744/1‑55 (MQ=255)
                    tGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGa                                                              <  1:3209557/59‑1 (MQ=255)
                      atatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCAt                                                                 >  1:2483606/1‑54 (MQ=255)
                      atatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCAt                                                                 >  1:1160846/1‑54 (MQ=255)
                      atatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGacac                                                           <  1:2452746/60‑1 (MQ=255)
                      atatTAACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATc                                                  <  1:2266520/69‑1 (MQ=255)
                       tatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTcccc                                                                    <  1:3254467/50‑1 (MQ=255)
                       tatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCgg                                                <  1:1678059/70‑1 (MQ=255)
                       tatTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCCGt                                               <  1:1926732/71‑1 (MQ=255)
                        atTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGTCACCCACACATCGGta                                              >  1:795956/1‑71 (MQ=255)
                        atTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCgg                                                <  1:1665479/69‑1 (MQ=255)
                        atTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              >  1:558292/1‑71 (MQ=255)
                           cACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTcccc                                                                    >  1:1790900/1‑46 (MQ=255)
                           cACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              >  1:1378896/1‑68 (MQ=255)
                            aCGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAg                                           <  1:2325358/70‑1 (MQ=255)
                            aCGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAg                                           <  1:1465726/70‑1 (MQ=255)
                              gcagcaATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGAt                                         >  1:3417693/1‑70 (MQ=255)
                              gcagcaATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATa                                        <  1:2929801/71‑1 (MQ=255)
                                agcaATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATaaa                                      <  1:3002880/71‑1 (MQ=255)
                                  caatgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGt                                    >  1:2884755/1‑71 (MQ=255)
                                    atgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGa                                                              <  1:2208449/43‑1 (MQ=255)
                                    atgcaatCCGTTCCGCCCATATCTCTTTTGCTTCCCCCATCGACACCCACACATCgg                                                <  1:3095397/57‑1 (MQ=255)
                                     tgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGacac                                                           >  1:4089998/1‑45 (MQ=255)
                                     tgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCcac                                                       <  1:2757391/49‑1 (MQ=255)
                                     tgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              >  1:3246022/1‑58 (MQ=255)
                                     tgcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAg                                 <  1:2892110/71‑1 (MQ=255)
                                      gcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGaca                                                            >  1:1106913/1‑43 (MQ=255)
                                      gcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGaca                                                            >  1:1116336/1‑43 (MQ=255)
                                      gcaatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAg                                 >  1:3278866/1‑70 (MQ=255)
                                       caatCCGTTCCGCCCATTTCTCTTTTGCTTccccc                                                                   >  1:712698/1‑35 (MQ=255)
                                       caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              >  1:442242/1‑56 (MQ=255)
                                       caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGt                                               >  1:3187149/1‑55 (MQ=255)
                                       caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGc                                <  1:35147/70‑1 (MQ=255)
                                       caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCa                               >  1:788911/1‑71 (MQ=255)
                                       caatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCa                               >  1:1749111/1‑71 (MQ=255)
                                        aatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGc                                <  1:2004618/69‑1 (MQ=255)
                                        aatCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAc                              <  1:1655502/71‑1 (MQ=255)
                                         atCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCAcc                             >  1:1118060/1‑71 (MQ=255)
                                          tCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCt                            >  1:1191696/1‑71 (MQ=255)
                                           ccGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCa                                  >  1:3883369/1‑64 (MQ=255)
                                           ccGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCa                                  >  1:1900714/1‑64 (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)
                                                 cGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCgg                                                <  1:1340670/44‑1 (MQ=255)
                                                    ccATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                  >  1:4035780/1‑71 (MQ=255)
                                                    ccATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCt                  >  1:4072501/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)
                                                       tttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATa                                        <  1:1956282/46‑1 (MQ=255)
                                                       tttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCa                         <  1:2886486/61‑1 (MQ=255)
                                                        ttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc              >  1:55714/1‑71 (MQ=255)
                                                        ttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc              >  1:3351464/1‑71 (MQ=255)
                                                        ttCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc              >  1:1179593/1‑71 (MQ=255)
                                                        ttCTCTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGCCCTTAAcg                      <  1:3108798/63‑2 (MQ=255)
                                                         tctcTTTTGCTTCCCCCATCGACACCCACACATCGGta                                              >  1:1988104/1‑38 (MQ=255)
                                                         tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATaaa                                      >  1:3149849/1‑46 (MQ=255)
                                                         tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGt                                    <  1:2301328/48‑1 (MQ=255)
                                                         tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTc                                   >  1:167274/1‑49 (MQ=255)
                                                         tctcTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgc              <  1:389013/70‑1 (MQ=255)
                                                         tctcTTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATAAAGTCCGCGccct                           >  1:1783061/1‑55 (MQ=255)
                                                          ctctTTTGCTTCCCCCATCGACACCCACACATCGGtata                                            >  1:3681552/1‑39 (MQ=255)
                                                          ctctTTTGCTTCCCCCATCGACACCCACACATCGGtata                                            >  1:4041049/1‑39 (MQ=255)
                                                          ctctTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgcg             <  1:3346098/70‑1 (MQ=255)
                                                          ctctTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTcgcg             <  1:2465857/70‑1 (MQ=255)
                                                          ctctTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGAt                                         >  1:4030659/1‑42 (MQ=25)
                                                          ctctTTGGCCTCGCCCATCGACACCCACACGTCGGTATAGATaaa                                      >  1:3599960/1‑45 (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)
                                                                      |                                                                    
GCATCATCTTGCTGTTCACCTGATATTCACGCAGCAATGCAATCCGTTCCGCCCATTTCTCTTTTGCTTCCCCCATCGACACCCACACATCGGTATAGATAAAGTCAGCACCTTCAACTCCCTTCGCGACATCTTCAGT  >  minE/2918912‑2919050

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: