Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,799,423 C→T 45.9% F83F (TTC→TTT proW → glycine betaine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,799,4230CT45.9% 11.9 / 59.4 74F83F (TTC→TTTproWglycine betaine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base T (13/21);  total (30/44)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

CCCGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTGGG  >  minE/1799360‑1799482
                                                               |                                                           
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:1328872/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:401891/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:2686485/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:48207/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:506369/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:1500843/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:2010604/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:2410659/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:778833/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:3210719/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:844553/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:2447540/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:84751/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:1053212/67‑1 (MQ=255)
cccGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:18247/67‑1 (MQ=255)
    cTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGt                                                     >  1:2362269/1‑68 (MQ=255)
    cTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGt                                                     >  1:2036069/1‑68 (MQ=255)
    cTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGt                                                     >  1:1175239/1‑68 (MQ=255)
    cTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTg                                                    >  1:2587917/1‑69 (MQ=255)
    cTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTg                                                    >  1:2341496/1‑69 (MQ=255)
    cTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTg                                                    >  1:81080/1‑69 (MQ=255)
    cTCGACAGTTGCGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTg                                                    >  1:1399301/1‑69 (MQ=255)
                gTCACTGAAGGGATCGACTGGGTCGTTACCCATTTTCGTCCCGTCTTTCAg                                                          <  1:71253/51‑1 (MQ=255)
                     tGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATta                                   <  1:92125/69‑1 (MQ=255)
                      gAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATtat                                  >  1:2008206/1‑69 (MQ=255)
                      gAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATta                                   <  1:445343/68‑1 (MQ=255)
                      gAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATta                                   >  1:2355491/1‑68 (MQ=255)
                         gggATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATtata                                 <  1:2160746/67‑1 (MQ=255)
                         gggATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATtata                                 <  1:1360746/67‑1 (MQ=255)
                         gggATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATtata                                 <  1:2570272/67‑1 (MQ=255)
                                   gggTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGtt                     <  1:2145836/69‑1 (MQ=255)
                                   gggTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGtt                     <  1:1753158/69‑1 (MQ=255)
                                   gggTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGtt                     <  1:850141/69‑1 (MQ=255)
                                   gggTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCATCAACGGtt                     <  1:366456/69‑1 (MQ=255)
                                       cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:2840038/69‑1 (MQ=255)
                                       cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:3426235/69‑1 (MQ=255)
                                       cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:2070328/69‑1 (MQ=255)
                                       cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:1361933/69‑1 (MQ=255)
                                       cGTTACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCa                 <  1:1231882/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:1072470/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1267579/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1477595/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1726204/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1773373/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:3267939/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:1784708/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:3130476/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:2966784/1‑69 (MQ=255)
                                         ttACCCATTTTCGTCCCGTCTTTCAGGGCGTGCGCGTTCCGGTTGATTACATCCTCAACGGTTTCCAGc               >  1:197805/1‑69 (MQ=255)
                                             ccATTTCCGTCCCGTCTTCCAGGGGGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:3155063/68‑1 (MQ=255)
                                             ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:2151737/68‑1 (MQ=255)
                                             ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:1294745/68‑1 (MQ=255)
                                             ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:1295739/68‑1 (MQ=255)
                                             ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:431994/68‑1 (MQ=255)
                                             ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:1657711/68‑1 (MQ=255)
                                             ccATTTCCGTCCCGTCTTCCAGGGCGTGCGCGGTCCGGTTGATTATATCCTCAACGGTTTCCAGCAAt            <  1:3260701/68‑1 (MQ=255)
                                                  tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:2643039/1‑69 (MQ=255)
                                                  tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:1229880/1‑69 (MQ=255)
                                                  tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:2711696/1‑69 (MQ=255)
                                                  tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:1736385/1‑69 (MQ=255)
                                                  tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:178044/1‑69 (MQ=255)
                                                  tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctgc      >  1:3258793/1‑69 (MQ=255)
                                                  tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctg       >  1:3022893/1‑68 (MQ=255)
                                                  tCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATtgctcc      >  1:2624381/1‑67 (MQ=255)
                                                      tcccgtcTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:743315/69‑1 (MQ=255)
                                                      tcccgtcTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:1192398/69‑1 (MQ=255)
                                                      tcccgtcTTCCAGGGCGTGCGCGGTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTggg  <  1:3062834/69‑1 (MQ=255)
                                                       cccgtcTTCCAGGGCTTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGc               <  1:601407/55‑1 (MQ=255)
                                                       cccgtcTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGc               <  1:453319/55‑1 (MQ=255)
                                                       cccgtcTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGc               <  1:1605890/55‑1 (MQ=255)
                                                              tCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGtt                     <  1:906136/42‑1 (MQ=255)
                                                               |                                                           
CCCGCTCGACAGTTGGGTCACTGAAGGGATCGACTGGGTCGTTACCCATTTCCGTCCCGTCTTCCAGGGCGTGCGCGTTCCGGTTGATTATATCCTCAACGGTTTCCAGCAATTGCTGCTGGG  >  minE/1799360‑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: