Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,396,732 A→G 58.6% Y152Y (TAT→TAC corA ← magnesium/nickel/cobalt transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,396,7320AG58.6% 33.7 / 56.7 82Y152Y (TAT→TACcorAmagnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base G (29/19);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTCGTCGTACTCATCGCCCTGATGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACTCG  >  minE/2396666‑2396798
                                                                  |                                                                  
cgcCTCGTCGTACTCATCGCCCTGATGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCtata                                                                  <  1:3156722/69‑1 (MQ=255)
cgcCTCGTCGTACTCATCGCCCTGATGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCtata                                                                  <  1:1838559/69‑1 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:2931996/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:2929664/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:2582398/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:3224385/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:334415/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:390435/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:2094385/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:2514215/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:2505345/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:737141/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:754465/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:811367/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:2407193/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:3058610/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCGGCTCAACTGCTACAGGTCGCTGTAGata                                                              >  1:1081559/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACCCAGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:671670/1‑69 (MQ=255)
    tcgtcgTACTCATCGCCCTGATGCCCTTCCATAATTACACGGCTCAACTGCTCCAGGTCGCTGTAGata                                                              >  1:500293/1‑69 (MQ=255)
               aTCGCCCTGATGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTc                                                   >  1:2553830/1‑69 (MQ=255)
                 cGCCCTGATGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCa                                                  <  1:1722650/68‑1 (MQ=255)
                 cGCCCTGATGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCa                                                  <  1:1592267/68‑1 (MQ=255)
                 cGCCCTGATGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCa                                                  <  1:837447/68‑1 (MQ=255)
                        aTGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCaa                                          >  1:986752/1‑69 (MQ=255)
                        aTGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCa                                           >  1:369299/1‑68 (MQ=255)
                          gCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCaa                                          <  1:2295926/67‑1 (MQ=255)
                          gCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCaa                                          <  1:1559378/67‑1 (MQ=255)
                          gCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCaa                                          <  1:564487/67‑1 (MQ=255)
                             cTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTGt                                      >  1:1089965/1‑68 (MQ=255)
                             cTTACATAATCACCCGGCTAAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTg                                       >  1:445365/1‑67 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATc                                                >  1:2564914/1‑57 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:3124588/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:3177351/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:2699624/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:708574/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:722376/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:2515749/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:1019320/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:1688257/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:1344265/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:163063/1‑69 (MQ=255)
                              ttCCATAATTACCCGGCTCAACTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACTGTTc                                    >  1:1244628/1‑69 (MQ=255)
                                     aTCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTc                                                         <  1:2172299/41‑1 (MQ=255)
                                     aTCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTc                                                         <  1:1671076/41‑1 (MQ=255)
                                     aTCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTc                                                         <  1:2992693/41‑1 (MQ=255)
                                     aTCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTc                                                         <  1:1916684/41‑1 (MQ=255)
                                     aTCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCa                                                  <  1:2237127/48‑1 (MQ=255)
                                     aTCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTg                             >  1:768883/1‑69 (MQ=255)
                                     aTCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTg                             >  1:1365775/1‑69 (MQ=255)
                                     aTCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTg                             >  1:2866103/1‑69 (MQ=255)
                                             gCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCa                                                  <  1:3410839/40‑1 (MQ=255)
                                                 taCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:157784/45‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:3320862/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:3225860/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:3368898/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:3119334/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:1453782/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:2910468/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:2909890/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:13449/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:281052/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:2706658/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:2085263/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:2122296/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:1195589/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:1012127/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:762847/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCCAACt                                        <  1:2281144/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATCTGCAAACt                                        <  1:3020960/46‑1 (MQ=255)
                                                 aaCTGCTCCAGGTCGCTGTAGATATTTTCAATTTCATATGCCAACt                                        <  1:1409183/46‑1 (MQ=39)
                                                    tGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATcc              <  1:3298439/69‑1 (MQ=255)
                                                    tGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATcc              <  1:1520204/69‑1 (MQ=255)
                                                           ggTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCagcagc        <  1:1013580/68‑1 (MQ=255)
                                                           ggTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCagcagc        <  1:215692/68‑1 (MQ=255)
                                                           ggTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCagcagc        <  1:2499003/68‑1 (MQ=255)
                                                           ggTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCCAACAGATCCagcagc        <  1:1640898/68‑1 (MQ=255)
                                                            gTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCagcagc        <  1:1138325/67‑1 (MQ=255)
                                                              cGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACt    >  1:2709450/1‑69 (MQ=255)
                                                              cGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACt    >  1:1236200/1‑69 (MQ=255)
                                                              cGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACt    >  1:122176/1‑69 (MQ=255)
                                                              cGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACt    >  1:2483874/1‑69 (MQ=255)
                                                              cGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACt    >  1:82330/1‑69 (MQ=255)
                                                              cGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACt    >  1:972373/1‑69 (MQ=255)
                                                              cGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACt    >  1:2359661/1‑69 (MQ=255)
                                                                 tataGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACTCg  <  1:152247/68‑1 (MQ=255)
                                                                 tataGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACTCg  <  1:2907124/68‑1 (MQ=255)
                                                                 tataGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACTCg  <  1:2091040/68‑1 (MQ=255)
                                                                 tataGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACTCg  <  1:1179002/68‑1 (MQ=255)
                                                                 tataGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACTCg  <  1:2378532/68‑1 (MQ=255)
                                                                  |                                                                  
CGCCTCGTCGTACTCATCGCCCTGATGCCCTTCCATAATCACCCGGCTCAACTGCTCCAGGTCGCTATAGATATTTTCAATTTCATCTGCCAACTGTTCGATTTTGGTTTCGAACAGATCCAGCAGCAACTCG  >  minE/2396666‑2396798

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: