Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,799,450 T→C 27.6% Y92Y (TAT→TAC proW → glycine betaine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,799,4500TC27.6% 66.3 / 25.5 65Y92Y (TAT→TACproWglycine betaine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (13/34);  new base C (13/5);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

GAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTGGG  >  minE/1799382‑1799482
                                                                    |                                
gAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATtat                                  >  1:2008206/1‑69 (MQ=255)
   gggATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATtata                                 <  1:1360746/67‑1 (MQ=255)
   gggATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATtata                                 <  1:2570272/67‑1 (MQ=255)
   gggATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATtata                                 <  1:2160746/67‑1 (MQ=255)
             gggTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGtt                     <  1:1753158/69‑1 (MQ=255)
             gggTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGtt                     <  1:2145836/69‑1 (MQ=255)
             gggTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGtt                     <  1:850141/69‑1 (MQ=255)
             gggTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCATCAACGGtt                     <  1:366456/69‑1 (MQ=255)
                 cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:2070328/69‑1 (MQ=255)
                 cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:2840038/69‑1 (MQ=255)
                 cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:1231882/69‑1 (MQ=255)
                 cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:3426235/69‑1 (MQ=255)
                 cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:1361933/69‑1 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCTTCAACGGTTTCCAGc               >  1:256411/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAg                >  1:3154430/1‑68 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:381918/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:3267939/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:2966784/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:3130476/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1072470/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:197805/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1784708/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1773373/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1726204/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1267579/1‑69 (MQ=255)
                   ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1477595/1‑69 (MQ=255)
                       ccATTTCCGTCCCGTCTTCCAGGGGGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:3155063/68‑1 (MQ=255)
                       ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:1294745/68‑1 (MQ=255)
                       ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:1295739/68‑1 (MQ=255)
                       ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:1657711/68‑1 (MQ=255)
                       ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:2151737/68‑1 (MQ=255)
                       ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:431994/68‑1 (MQ=255)
                       ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGGTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:3260701/68‑1 (MQ=255)
                            tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:2711696/1‑69 (MQ=255)
                            tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:2643039/1‑69 (MQ=255)
                            tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:1229880/1‑69 (MQ=255)
                            tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:3258793/1‑69 (MQ=255)
                            tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:1736385/1‑69 (MQ=255)
                            tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:178044/1‑69 (MQ=255)
                            tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctg       >  1:3022893/1‑68 (MQ=255)
                            tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctcc      >  1:2624381/1‑67 (MQ=255)
                                tcccgtcTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:1192398/69‑1 (MQ=255)
                                tcccgtcTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:743315/69‑1 (MQ=255)
                                tcccgtcTTCCAGGGCGTGCGCGGTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:3062834/69‑1 (MQ=255)
                                 cccgtcTTCCAGGGCTTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGc               <  1:601407/55‑1 (MQ=255)
                                 cccgtcTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGc               <  1:1605890/55‑1 (MQ=255)
                                 cccgtcTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGc               <  1:453319/55‑1 (MQ=255)
                                        tCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGtt                     <  1:906136/42‑1 (MQ=255)
                                               cGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCa              >  1:1395975/1‑42 (MQ=255)
                                               cGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCa              <  1:2568492/42‑1 (MQ=255)
                                               cGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCa              <  1:2515271/42‑1 (MQ=255)
                                               cGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCa              >  1:2209946/1‑42 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:368088/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:1128236/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:464369/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:2989008/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:2747516/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:2690359/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:2688848/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:2180151/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:2046439/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:1850369/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:1644230/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:1463773/50‑1 (MQ=255)
                                                   cgcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:1213964/50‑1 (MQ=255)
                                                    gcgTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:2870807/49‑1 (MQ=255)
                                                               gATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  >  1:2543280/1‑38 (MQ=255)
                                                               gATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  >  1:2275724/1‑38 (MQ=255)
                                                               gATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  >  1:2227828/1‑38 (MQ=255)
                                                                    |                                
GAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTGGG  >  minE/1799382‑1799482

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: