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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,3870AG20.4% 277.9 / 65.0 142G103G (GGT→GGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (68/45);  new base G (20/9);  total (88/54)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 144 total aligned reads displayed.

CCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTG  >  minE/2919318‑2919453
                                                                     |                                                                  
ccATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATAc                                                                   >  1:2680694/1‑71 (MQ=255)
ccATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATAc                                                                   >  1:402945/1‑71 (MQ=255)
ccATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATAc                                                                   >  1:5014182/1‑71 (MQ=255)
 cATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATa                                                                    <  1:4792160/69‑1 (MQ=255)
  aTACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCg                                                                 <  1:2367806/71‑1 (MQ=255)
  aTACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCg                                                                 <  1:1697541/71‑1 (MQ=255)
    aCCGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGt                                                                >  1:1745219/1‑70 (MQ=255)
    aCCGGCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTc                                                               <  1:4298337/71‑1 (MQ=255)
     ccGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTc                                                               <  1:3940701/70‑1 (MQ=255)
     ccGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                              >  1:4518782/1‑71 (MQ=255)
       ggCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATa                                                            >  1:3176107/1‑71 (MQ=255)
       ggCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATa                                                            <  1:2234140/71‑1 (MQ=255)
       ggCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATa                                                            >  1:4488857/1‑71 (MQ=255)
        gCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATAc                                                           <  1:3546077/71‑1 (MQ=255)
        gCACGCCCGCATACTGCGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATa                                                            >  1:5183585/1‑70 (MQ=255)
         cACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACa                                                          >  1:10161/1‑71 (MQ=255)
             cTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGc                                                       >  1:400834/1‑70 (MQ=255)
             cTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCg                                                      <  1:4496304/71‑1 (MQ=255)
             cTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCg                                                      <  1:5055038/71‑1 (MQ=255)
                gCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                     <  1:5155006/69‑1 (MQ=255)
                  gTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGAcc                                                   >  1:4983560/1‑69 (MQ=255)
                   tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAg                                                >  1:2112056/1‑71 (MQ=255)
                   tATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATAATGAATACCGTCATACATGCGACCAAg                                                >  1:5192670/1‑71 (MQ=255)
                    aTTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAg                                                >  1:1579204/1‑70 (MQ=255)
                    aTTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAg                                                >  1:776690/1‑70 (MQ=255)
                     ttCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                              >  1:1888607/1‑55 (MQ=255)
                      tCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa                                              >  1:5674224/1‑70 (MQ=255)
                       ccGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                     >  1:3163553/1‑62 (MQ=255)
                       ccGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcg                                            >  1:832260/1‑71 (MQ=255)
                       ccGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAc                                             <  1:74845/70‑1 (MQ=255)
                       ccGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAc                                             <  1:949554/70‑1 (MQ=255)
                        cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACa                                                          >  1:4034677/1‑56 (MQ=255)
                        cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACa                                                          >  1:975375/1‑56 (MQ=255)
                        cGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcgc                                           <  1:4914136/71‑1 (MQ=255)
                        cGCCAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAc                                             <  1:1889663/69‑1 (MQ=255)
                         gCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGt                                                                <  1:305113/49‑1 (MQ=255)
                         gCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGt                                                                <  1:2428923/49‑1 (MQ=255)
                           cAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                              >  1:1321866/1‑49 (MQ=255)
                           cAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCa                                                              >  1:550620/1‑49 (MQ=255)
                           cAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCggg                                        <  1:4073582/71‑1 (MQ=255)
                           cAGCGTTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAAcccgcg                                        <  1:4149179/71‑3 (MQ=255)
                            agcgtTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACATCCGCCCGAGAACCCgcgc                                       <  1:2994542/68‑1 (MQ=255)
                             gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCggg                                        >  1:1101333/1‑69 (MQ=255)
                             gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCggg                                        >  1:846323/1‑69 (MQ=255)
                             gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCggg                                        >  1:3615736/1‑69 (MQ=255)
                             gtgtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCggg                                        >  1:1392164/1‑69 (MQ=255)
                               gtTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACAtc                                                        >  1:2186047/1‑50 (MQ=255)
                               gtTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACAtc                                                        >  1:4166330/1‑50 (MQ=255)
                               gtTTCGACCACTTCCTGGCCGTGACCGCGATACTGAATGCCGTCATACAtc                                                        >  1:1706309/1‑50 (MQ=255)
                               gtTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgt                                    >  1:115740/1‑71 (MQ=255)
                                tttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGa                                                     <  1:1687128/53‑1 (MQ=255)
                                tttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtg                                   >  1:261763/1‑71 (MQ=255)
                                tttCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgt                                    <  1:2736729/70‑1 (MQ=255)
                                tttCGACAATCTCCTGACAATAGCCGCGATACTGAATACCGTCATAc                                                           <  1:1045940/47‑1 (MQ=255)
                                  tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa                                              >  1:137224/1‑58 (MQ=255)
                                  tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTc                                 >  1:161096/1‑71 (MQ=255)
                                  tCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTc                                 >  1:4688723/1‑71 (MQ=255)
                                   cGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCa                                              >  1:3830768/1‑57 (MQ=255)
                                    gACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCa                                      >  1:4847672/1‑64 (MQ=255)
                                    gACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                               <  1:2301822/71‑1 (MQ=255)
                                    gACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCtt                               <  1:5474177/71‑1 (MQ=255)
                                       aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCAcgc                                           >  1:2663031/1‑56 (MQ=255)
                                       aaTCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTaa                            <  1:1606965/71‑1 (MQ=255)
                                             cTGACCATAGCCGCTATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                      >  1:4619687/1‑71 (MQ=255)
                                             cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                          <  1:1424075/67‑1 (MQ=255)
                                             cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                          >  1:1537469/1‑67 (MQ=255)
                                             cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                        >  1:2821532/1‑69 (MQ=255)
                                             cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                      >  1:2343272/1‑71 (MQ=255)
                                             cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGAGTCTTTAATCGAct                      >  1:772278/1‑71 (MQ=255)
                                              tGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActc                     >  1:5559379/1‑71 (MQ=255)
                                                aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtgt                                  >  1:2735544/1‑56 (MQ=255)
                                                 ccATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                        >  1:2899789/1‑65 (MQ=255)
                                                   aTAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                >  1:612447/1‑71 (MQ=255)
                                                    tAGCGGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                <  1:1430247/70‑1 (MQ=255)
                                                    tAGCCGCGATACTGAATACCGTCATACATGCGACCAAg                                                <  1:3016552/38‑1 (MQ=255)
                                                    tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                >  1:3135233/1‑70 (MQ=255)
                                                      gCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                        >  1:2471840/1‑60 (MQ=255)
                                                       ccGCGATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATgcc             >  1:3837458/1‑68 (MQ=255)
                                                       ccGCGATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATgcc             >  1:118750/1‑68 (MQ=255)
                                                          cgATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGAct                      >  1:2446098/1‑58 (MQ=255)
                                                          cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                   <  1:3285013/61‑1 (MQ=255)
                                                          cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                   <  1:1306294/61‑1 (MQ=255)
                                                          cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt         <  1:705228/71‑1 (MQ=255)
                                                            aTACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                          >  1:4883388/1‑52 (MQ=255)
                                                             tACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTaa                            >  1:134308/1‑49 (MQ=255)
                                                             tACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtg                                   <  1:1824800/42‑1 (MQ=255)
                                                             tACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                          <  1:75454/51‑1 (MQ=255)
                                                             tACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt       <  1:5450936/70‑1 (MQ=255)
                                                             tACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt       <  1:548513/70‑1 (MQ=255)
                                                              aCTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg     >  1:5162771/1‑71 (MQ=255)
                                                               cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtgt                                  <  1:4325460/41‑1 (MQ=255)
                                                               cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAgtgt                                  <  1:736227/41‑1 (MQ=255)
                                                               cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                        >  1:2246479/1‑51 (MQ=255)
                                                               cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                        >  1:3373813/1‑51 (MQ=255)
                                                               cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg     <  1:63655/70‑1 (MQ=255)
                                                               cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg     <  1:3782015/70‑1 (MQ=255)
                                                               cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGc    >  1:5073876/1‑71 (MQ=255)
                                                                   aTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa              <  1:4920252/57‑1 (MQ=255)
                                                                    tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                   <  1:1953064/51‑1 (MQ=255)
                                                                    tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt       >  1:3357756/1‑63 (MQ=255)
                                                                     |                                                                  
CCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTG  >  minE/2919318‑2919453

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: