Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,727,498 A→G 30.8% T244T (ACT→ACC purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,727,4980AG30.8% 83.4 / 45.5 94T244T (ACT→ACCpurLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (40/25);  new base G (11/18);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 4.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCGATAACCCAGTCGGC  >  minE/1727436‑1727560
                                                              |                                                              
ccTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTAttt                                                          >  1:293882/1‑69 (MQ=255)
   tCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTg                                                       <  1:648307/69‑1 (MQ=255)
   tCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTg                                                       <  1:3202163/69‑1 (MQ=255)
   tCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTg                                                       <  1:2426482/69‑1 (MQ=255)
       ttACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:1861824/69‑1 (MQ=255)
       ttACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:1816331/69‑1 (MQ=255)
       ttACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:1988204/69‑1 (MQ=255)
       ttACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:2089111/69‑1 (MQ=255)
       ttACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:2094103/69‑1 (MQ=255)
       ttACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:791225/69‑1 (MQ=255)
       ttACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:2639849/69‑1 (MQ=255)
       ttACGGCGGCGTTATCTTTAAAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:324657/69‑1 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                                                 >  1:847802/1‑69 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                                                 >  1:1994026/1‑69 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                                                 >  1:1966230/1‑69 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                                                 >  1:593140/1‑69 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                                                 >  1:2103482/1‑69 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                                                 >  1:1730375/1‑69 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcatc                                                 >  1:866837/1‑69 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcat                                                  >  1:2877035/1‑68 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcat                                                  >  1:3104646/1‑68 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatcat                                                  >  1:799949/1‑68 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatc                                                    >  1:2464182/1‑66 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatc                                                    >  1:265102/1‑66 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatc                                                    >  1:2921309/1‑66 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatc                                                    >  1:2081112/1‑66 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGataat                                                  >  1:3041343/1‑68 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGat                                                     >  1:2184566/1‑65 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGat                                                     >  1:1357066/1‑65 (MQ=255)
         aCGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTAATTTTGata                                                    >  1:2607980/1‑65 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              <  1:2573931/69‑1 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              <  1:739945/69‑1 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              <  1:425142/69‑1 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              <  1:909033/69‑1 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              <  1:1315488/69‑1 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              <  1:665539/69‑1 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              <  1:3157482/69‑1 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              <  1:1995283/69‑1 (MQ=255)
            gcggcgTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTg                                              >  1:1848697/1‑69 (MQ=255)
                cgTTATCTTTATCAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatca                                                   <  1:2357112/60‑1 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGatc                                                    >  1:550915/1‑55 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCg                                       >  1:1459250/1‑68 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:2886570/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:125265/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:2887157/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:446968/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:3279848/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:3311711/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:2295700/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:2034295/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:1637001/1‑69 (MQ=255)
                    aTCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGa                                      >  1:848091/1‑69 (MQ=255)
                               cAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGgctgc                            <  1:2557372/68‑1 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGat                                                     >  1:90758/1‑35 (MQ=25)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:3142353/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:44021/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:1168218/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:2951969/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:2606877/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:774261/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:2174095/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:2070427/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:903088/1‑68 (MQ=255)
                                       aCGTGGTCTGGCGTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCAccg                   >  1:1103342/1‑68 (MQ=255)
                                            aTCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCGAt              <  1:556997/69‑1 (MQ=255)
                                            aTCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCGAt              <  1:2741503/69‑1 (MQ=255)
                                            aTCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCGAt              <  1:1630729/69‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:3338136/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:913536/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:330577/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:276468/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:1288209/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:1365691/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:1514339/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:1569854/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:1785623/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:2065498/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:2490660/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:2819495/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:2832330/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:2886599/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:2936407/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:3018224/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:3122211/67‑1 (MQ=255)
                                                   gTGGTTTCGAAGGTGTTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCGTCGATAAtcc         <  1:3271764/67‑1 (MQ=255)
                                                        ttCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGgctgc                            <  1:2561360/43‑1 (MQ=255)
                                                         tCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCg                >  1:567220/1‑54 (MQ=255)
                                                         tCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCg                >  1:3286868/1‑54 (MQ=255)
                                                         tCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCg                >  1:2504459/1‑54 (MQ=255)
                                                         tCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCg                >  1:2141451/1‑54 (MQ=255)
                                                         tCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCg                >  1:3341856/1‑54 (MQ=255)
                                                         tCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCg                >  1:1297068/1‑54 (MQ=255)
                                                         tCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCGATAACCCAGTCGGc  >  1:2903214/1‑68 (MQ=255)
                                                         tCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCGATAACCCAGTCGGc  >  1:3052325/1‑68 (MQ=255)
                                                              |                                                              
CCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACGTGATCTGGCGTGGTTTCGAAAGTATTTTTGATCATCTTGAACAGCGATTTCGGCTGCTGTTCACCATCGATAACCCAGTCGGC  >  minE/1727436‑1727560

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: