Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 165,160 A→G 15.5% F772F (TTT→TTC glnD ← uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE165,1600AG15.5% 207.2 / 29.4 84F772F (TTT→TTCglnDuridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (45/26);  new base G (8/5);  total (53/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCC  >  minE/165091‑165229
                                                                     |                                                                     
tAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCaa                                                                      <  1:191152/71‑1 (MQ=255)
tAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCaa                                                                      <  1:1479125/71‑1 (MQ=255)
 aaTTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCgaa                                                                     >  1:1999818/1‑71 (MQ=255)
 aaTTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCgaa                                                                     >  1:2528351/1‑71 (MQ=255)
 aaTTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCaaa                                                                     >  1:823044/1‑71 (MQ=255)
  aTTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAAc                                                                    >  1:1043238/1‑71 (MQ=255)
  aTTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAAc                                                                    >  1:1080932/1‑71 (MQ=255)
     tGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCg                                                                  <  1:3705251/70‑1 (MQ=255)
         gggTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATaa                                                             <  1:1321850/71‑1 (MQ=255)
            ttGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCGAAcc                                                                   >  1:1866039/1‑62 (MQ=255)
                 gACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGa                                                     <  1:3512547/71‑1 (MQ=255)
                   cGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACg                                                   <  1:3401362/71‑1 (MQ=255)
                    ggggCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCGAAc                                                                    >  1:2307930/1‑53 (MQ=255)
                    ggggCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCGAACCGAATAACCTCATGACGa                                                  >  1:2473234/1‑71 (MQ=255)
                    ggggCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTc                                                         <  1:163446/64‑1 (MQ=255)
                    ggggCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACgc                                                  >  1:1158142/1‑70 (MQ=255)
                     gggCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGaa                                                                <  1:3188378/56‑1 (MQ=255)
                           cGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGc                                             <  1:2178426/69‑1 (MQ=255)
                                gCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGa                                                                 <  1:9142/44‑1 (MQ=255)
                                gCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGc                                      >  1:788891/1‑71 (MQ=255)
                                 cTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGAt                                                 >  1:956412/1‑59 (MQ=255)
                                 cTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCg                                     <  1:710533/71‑1 (MQ=255)
                                   gCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGa                                          >  1:1404954/1‑64 (MQ=255)
                                    ccAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAAc                                                            <  1:3572076/45‑1 (MQ=255)
                                    ccAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCgg                                    <  1:1195173/69‑1 (MQ=255)
                                    ccAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCgg                                    <  1:2389576/69‑1 (MQ=255)
                                     cAGCTACTTTGCGTCAGTACTTGCTCCAGACCGAACCGAATAACCTCATGACGATCTGCGGACAGCg                                     >  1:197932/1‑67 (MQ=255)
                                     cAGCTACCTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGtt                                 <  1:4037525/71‑1 (MQ=255)
                                       gCTACTTTGCGTCAGTACTTGCTCCAGACCGAACCGAATAACCTCa                                                        >  1:4019877/1‑46 (MQ=255)
                                       gCTACTTTGCGTCAGTACTTGCTCCAGACCGAACCGAATAACCTCa                                                        >  1:3124682/1‑46 (MQ=255)
                                        cTACTTTGCGTCAGTACTTGCTCCAGACCGAACCg                                                                  <  1:1978446/35‑1 (MQ=255)
                                        cTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGc                               >  1:1189604/1‑70 (MQ=255)
                                         tACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCa                             >  1:1086948/1‑71 (MQ=255)
                                           cTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAAcc                                                           <  1:3262625/39‑1 (MQ=255)
                                            tttGCGTCAGTACTTGCTCCAGACCGAACCGAATAACCTc                                                         <  1:4061944/40‑1 (MQ=255)
                                            tttGCGTCAGTACTTGCTCCAGACCAAACCGAATaa                                                             >  1:3875807/1‑36 (MQ=255)
                                            tttGCGTCAGTACTTGCTCCAGACCAAACCGAATaa                                                             >  1:1418816/1‑36 (MQ=255)
                                             ttGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCAt                            >  1:520477/1‑68 (MQ=255)
                                             ttGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCAt                            >  1:2405855/1‑68 (MQ=255)
                                              tGCGTCAGTACTTGCTCCAGACCGAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCggg                        >  1:829070/1‑71 (MQ=255)
                                              tGCGTCAGTACTTGCTCCAGACCGAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCg                          <  1:2024253/69‑1 (MQ=255)
                                              tGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCAt                                                       >  1:647517/1‑40 (MQ=255)
                                              tGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTg                                              >  1:2296781/1‑49 (MQ=255)
                                              tGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCggg                        >  1:3716622/1‑71 (MQ=255)
                                               gCGTCAGTACTTGCTCCAGACCGAACCGAATAACCTCATg                                                      <  1:1501609/40‑1 (MQ=255)
                                               gCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCggg                        <  1:718976/70‑1 (MQ=255)
                                                 gTCAGTACTTGCTCCAGACCGAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTc                     <  1:879939/71‑1 (MQ=255)
                                                    aGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCgg                                    >  1:2246532/1‑53 (MQ=255)
                                                    aGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTc                     >  1:1314075/1‑68 (MQ=255)
                                                    aGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTc                     >  1:1832863/1‑68 (MQ=255)
                                                    aGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAg                  >  1:2871551/1‑71 (MQ=255)
                                                     gTACTTGCTCCAGACCAAACCGAATAACCTCATGACg                                                   <  1:3519115/37‑1 (MQ=255)
                                                      tACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCg                                            <  1:1020182/43‑1 (MQ=255)
                                                      tACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAg                                       <  1:1954848/48‑1 (MQ=255)
                                                       aCTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGtt                                 >  1:104517/1‑53 (MQ=255)
                                                       aCTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCAc               >  1:2501692/1‑71 (MQ=255)
                                                       aCTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCAc               >  1:3526504/1‑71 (MQ=255)
                                                        cTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCg                                     >  1:912531/1‑48 (MQ=255)
                                                        cTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGcc                              >  1:732501/1‑55 (MQ=255)
                                                        cTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGcc                              >  1:3254840/1‑55 (MQ=255)
                                                        cTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGcc                              >  1:2109353/1‑55 (MQ=255)
                                                        cTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCAc               >  1:1821098/1‑70 (MQ=255)
                                                         ttGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGAc                                         <  1:737705/43‑1 (MQ=255)
                                                          tGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGa             >  1:2356810/1‑70 (MQ=255)
                                                           gCTCCAGACCGAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGAt            >  1:110264/1‑70 (MQ=255)
                                                           gCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGAc                                         >  1:2897214/1‑41 (MQ=255)
                                                           gCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCggg                                   >  1:3981773/1‑47 (MQ=255)
                                                           gCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGc                 >  1:856705/1‑65 (MQ=255)
                                                            cTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATa           <  1:4011876/70‑1 (MQ=255)
                                                             tCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATa           >  1:3507984/1‑69 (MQ=255)
                                                             tCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATa           >  1:524350/1‑69 (MQ=255)
                                                               cAGACCAACCCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCa                   >  1:3454120/1‑59 (MQ=255)
                                                               cAGACCAAACCGAATAACCTCATGACGATCTGCGGACAg                                       >  1:3946739/1‑39 (MQ=255)
                                                               cAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATa           <  1:380857/67‑1 (MQ=255)
                                                               cAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAgg       >  1:2995590/1‑71 (MQ=255)
                                                                aGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAg        <  1:2769988/69‑1 (MQ=255)
                                                                aGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAg        <  1:1525584/69‑1 (MQ=255)
                                                                 gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCggg                                   <  1:3880923/41‑1 (MQ=255)
                                                                 gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTc                     >  1:904916/1‑55 (MQ=255)
                                                                 gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTc                     >  1:2084564/1‑55 (MQ=255)
                                                                 gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCa                   >  1:2653523/1‑57 (MQ=255)
                                                                 gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCa                   >  1:1290153/1‑57 (MQ=255)
                                                                 gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTa     >  1:3114407/1‑71 (MQ=255)
                                                                 gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTa     >  1:1952879/1‑71 (MQ=255)
                                                                 gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTa     >  1:821253/1‑71 (MQ=255)
                                                                   ccAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGt      >  1:1626608/1‑68 (MQ=255)
                                                                   ccAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATc   <  1:2697189/71‑1 (MQ=255)
                                                                    cAAACCGAATAACCTCATGACGGTCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATcc  <  1:2733132/71‑1 (MQ=255)
                                                                    cAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATc   >  1:1503416/1‑70 (MQ=255)
                                                                     aaaCCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGt      <  1:2204244/66‑1 (MQ=255)
                                                                     |                                                                     
TAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCC  >  minE/165091‑165229

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: