Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,918,802 A→C 6.6% G298G (GGT→GGG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,918,8020AC 6.6% 237.4 / 9.2 91G298G (GGT→GGGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (43/42);  new base C (4/2);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

ATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAACGCGGGC  >  minE/2918734‑2918871
                                                                    |                                                                     
aTACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATAcc                                                                     >  1:4716868/1‑71 (MQ=255)
  ccGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGc                                                                   >  1:1221099/2‑71 (MQ=255)
   cGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGc                                                                   >  1:4455993/1‑70 (MQ=255)
   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc                                                                   >  1:2264626/1‑70 (MQ=255)
   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc                                                                   >  1:2323733/1‑70 (MQ=255)
   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc                                                                   <  1:1832383/70‑1 (MQ=255)
   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc                                                                   >  1:1584815/1‑70 (MQ=255)
   cGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGc                                                                   >  1:4756248/1‑70 (MQ=255)
      ttttCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGCCGTg                                                               <  1:818793/71‑1 (MQ=255)
      ttttCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTCACCTCCATCCCGCCGTg                                                               <  1:2992857/71‑1 (MQ=255)
      ttttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATg                                                               <  1:1080511/71‑1 (MQ=255)
        ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATg                                                               >  1:627066/1‑69 (MQ=255)
        ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                                                              <  1:5079649/70‑1 (MQ=255)
        ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                                                              <  1:1652947/70‑1 (MQ=255)
        ttCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                                                              <  1:132974/70‑1 (MQ=255)
          ccGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa                                                             <  1:2975957/69‑1 (MQ=255)
          ccGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAg                                                            <  1:3210685/70‑1 (MQ=255)
             ccTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa                                                        <  1:4915402/71‑1 (MQ=255)
              cTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa                                                        >  1:1270875/1‑70 (MQ=255)
              cTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCa                                                        >  1:4355450/1‑70 (MQ=255)
                 aTCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                                                              <  1:3617310/61‑1 (MQ=255)
                 aTCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa                                                             <  1:5396164/62‑1 (MQ=255)
                 aTCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAt                                                     <  1:4201521/70‑1 (MQ=255)
                 aTCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAt                                                     >  1:2346855/1‑70 (MQ=255)
                  tCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAt                                                     <  1:487271/69‑1 (MQ=255)
                   cAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAg                                                            >  1:5607975/1‑61 (MQ=255)
                    aaaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAgg                                                           >  1:2195743/1‑61 (MQ=255)
                    aaaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttctt                                                 >  1:1040042/1‑71 (MQ=255)
                    aaaaaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAAttctt                                                 >  1:1574564/1‑71 (MQ=255)
                       aaCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCg                                              >  1:2308722/1‑71 (MQ=255)
                        aCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATg                                                               >  1:2314143/1‑53 (MQ=255)
                        aCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGc                                             >  1:234154/1‑71 (MQ=255)
                        aCAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGc                                             <  1:2687385/71‑1 (MQ=255)
                         cAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTcc                                               >  1:2137888/1‑68 (MQ=255)
                          aaTGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCa                                           >  1:4400442/1‑71 (MQ=255)
                           aTGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCAt                                          >  1:1178884/1‑71 (MQ=255)
                             gCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATa                                                                       >  1:254083/1‑40 (MQ=255)
                                    gcAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGcc                                 <  1:3290718/71‑1 (MQ=255)
                                       gATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTc                                      >  1:5507081/1‑63 (MQ=255)
                                        aTTCGAAGACCTCATCAGTGACTTCCATACCGCCATGt                                                              <  1:1340324/38‑1 (MQ=255)
                                        aTTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCaaa                                                      <  1:1511454/46‑1 (MQ=255)
                                        aTTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCaaa                                                      <  1:3330796/46‑1 (MQ=255)
                                        aTTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCaaa                                                      <  1:2010187/46‑1 (MQ=255)
                                        aTTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGcc                                 <  1:3529722/67‑1 (MQ=255)
                                         ttCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGcg                            <  1:4182252/71‑1 (MQ=255)
                                          tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa                                                             <  1:3882795/37‑1 (MQ=255)
                                          tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTa                                                             <  1:1340242/37‑1 (MQ=255)
                                          tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGcc                                 >  1:4531460/1‑65 (MQ=255)
                                          tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGc                             >  1:434607/1‑69 (MQ=255)
                                          tCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGc                             >  1:3393762/1‑69 (MQ=255)
                                                 ccTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTgg                    <  1:2863999/71‑1 (MQ=255)
                                                    catcaGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGgtcg                 >  1:2329362/1‑71 (MQ=255)
                                                     atcaGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGc                                  >  1:336684/1‑53 (MQ=255)
                                                     atcaGTAACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGttt                      >  1:9642/1‑65 (MQ=255)
                                                     atcaGTAACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGttt                      >  1:3158/1‑65 (MQ=255)
                                                      tcaGTGACTTCCATACCGCCATGTAGGCCAAAttct                                                  <  1:2308195/36‑1 (MQ=255)
                                                      tcaGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGttt                      <  1:1716639/64‑1 (MQ=255)
                                                       caGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCa              >  1:3211959/1‑71 (MQ=255)
                                                        aGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCAt             >  1:3066760/1‑71 (MQ=255)
                                                        aGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCAt             >  1:3955391/1‑71 (MQ=255)
                                                        aGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCAt             <  1:4536450/71‑1 (MQ=255)
                                                         gTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCtt                                    <  1:262123/47‑1 (MQ=255)
                                                         gTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGgtcg                 >  1:5034903/1‑66 (MQ=255)
                                                          tGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCAttt                                        <  1:1504065/42‑1 (MQ=255)
                                                          tGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGtt                       <  1:1916432/59‑1 (MQ=255)
                                                          tGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGtt                       <  1:272027/59‑1 (MQ=255)
                                                          tGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTg                     <  1:2790097/61‑1 (MQ=255)
                                                          tGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGa           <  1:1019896/71‑1 (MQ=255)
                                                            aCTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGc                             >  1:1321065/1‑51 (MQ=255)
                                                             cTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAAc        <  1:4055129/71‑1 (MQ=255)
                                                             cTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGAGTCGTTTGGTCGTCATGAAAc        <  1:841406/71‑1 (MQ=255)
                                                              ttCCATACCGCCATGTAGGCCAAATTCTTCCGCCAtttt                                       <  1:2863700/39‑1 (MQ=255)
                                                              ttCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGtt                       <  1:586546/55‑1 (MQ=255)
                                                                ccATCCCGCCGTGCAGATCGAACTCCTTCGCCATCTGCTTGCCGAGCGTAGTCTGGTCGTCATGGAAcgcc     >  1:4016186/1‑70 (MQ=255)
                                                                ccATCCCGCCGTGCAGATCGAACTCCTTCGCCATCTGCTTGCCGAGCGTAGTCTGGTCGTCATGGAAcgcc     >  1:2925773/1‑70 (MQ=255)
                                                                ccATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTg                     <  1:1202950/55‑1 (MQ=255)
                                                                ccATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGaaa         >  1:4303145/1‑67 (MQ=255)
                                                                ccATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGaaa         >  1:1570089/1‑67 (MQ=255)
                                                                ccATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAAcgcg     >  1:309115/1‑71 (MQ=255)
                                                                ccATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAAcgcg     >  1:1689541/1‑71 (MQ=255)
                                                                ccATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAAcg       <  1:2869432/69‑1 (MQ=255)
                                                                ccATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAAcg       <  1:2862641/69‑1 (MQ=255)
                                                                ccATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAAcg       <  1:4278395/69‑1 (MQ=255)
                                                                  aTACCGCCATGTAGGCCAAATTCTTCGGCCATTTTCTTGcc                                 >  1:5351087/1‑41 (MQ=255)
                                                                  aTACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCAt             >  1:4436874/1‑61 (MQ=255)
                                                                  aTACCGCCATGTAGGCCAAATTCTTCCGCAATTTTCTTGcc                                 >  1:2218345/1‑41 (MQ=255)
                                                                   tACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGcgtc                          >  1:4178180/1‑47 (MQ=255)
                                                                   tACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGcgtc                          >  1:3174864/1‑47 (MQ=255)
                                                                   tACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGcgtc                          >  1:1820049/1‑47 (MQ=255)
                                                                   tACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGa           <  1:140151/62‑1 (MQ=255)
                                                                   tACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAACGCggg   >  1:266001/1‑70 (MQ=255)
                                                                   tACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAACGCGggc  <  1:4653536/71‑1 (MQ=255)
                                                                    aCCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAACGCggg   >  1:4202450/1‑69 (MQ=255)
                                                                    |                                                                     
ATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTAGGCCAAATTCTTCCGCCATTTTCTTGCCAAGCGTCGTTTGGTCGTCATGAAACGCGGGC  >  minE/2918734‑2918871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: