Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 200,183 A→T 25.6% G88G (GGT→GGA metI ← DL‑methionine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE200,1830AT25.6% 162.5 / 60.4 113G88G (GGT→GGAmetIDL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (45/39);  new base T (13/16);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

Only 100 of 118 total aligned reads displayed.

CCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGG  >  minE/200125‑200248
                                                          |                                                                 
ccATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           <  1:3214049/67‑1 (MQ=255)
ccATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           <  1:454383/67‑1 (MQ=255)
ccATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATg                                                         <  1:2320318/69‑1 (MQ=255)
ccATACGGGCAATAAACGGTGCTGCACCAACGGTTAAAGGAACAATCGCTGCCTGCAAACCAATCGATg                                                         <  1:3380867/69‑1 (MQ=255)
 cATACGGGCAATAAACGGTTCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           <  1:1844182/66‑1 (MQ=255)
     cGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTAcc                                                     >  1:3125384/1‑68 (MQ=255)
     cGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTAcc                                                     >  1:2662390/1‑68 (MQ=255)
     cGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTAcc                                                     >  1:2445890/1‑68 (MQ=255)
       ggCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTG‑‑‑GCAAACCAATCGATGTACCGACAAt                                               <  1:3017627/69‑1 (MQ=255)
          aaTAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAAt                                               <  1:1123501/69‑1 (MQ=255)
          aaTAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAAt                                               <  1:712778/69‑1 (MQ=255)
             aaaCGGTGCTGCACTAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           >  1:2815068/1‑54 (MQ=255)
             aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           >  1:1461093/1‑54 (MQ=255)
             aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           >  1:446263/1‑54 (MQ=255)
             aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           >  1:345130/1‑54 (MQ=255)
             aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           >  1:1145661/1‑54 (MQ=255)
             aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           >  1:1347336/1‑54 (MQ=255)
             aaaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGa                                                           >  1:733147/1‑54 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGt                                                        >  1:2341084/1‑56 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGt                                                        >  1:2096489/1‑56 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGt                                                        >  1:2609779/1‑56 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGt                                                        >  1:2757215/1‑56 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGt                                                        >  1:779453/1‑56 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGt                                                        >  1:848012/1‑56 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGt                                                        >  1:975278/1‑56 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAAcg                                           >  1:535060/1‑69 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAAcg                                           >  1:2147814/1‑69 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAAcg                                           >  1:939316/1‑69 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAAcg                                           >  1:2990593/1‑69 (MQ=255)
              aaCGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCCGCAAACCAATCGATGTACCGACAATAAcg                                           >  1:3272196/1‑69 (MQ=255)
                         aCCAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:2910722/69‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:3350985/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:2772438/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:3048347/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:283613/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:31016/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:2513036/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:3016413/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:1042138/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:3259580/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:2193388/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:1113317/68‑1 (MQ=255)
                          ccAACGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:678401/68‑1 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:152229/66‑1 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:1983819/66‑1 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGaa                              >  1:59427/1‑68 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:2742998/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:1742325/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:2409258/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:180123/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:1567067/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:1412970/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:1310443/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:294124/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:2970909/1‑69 (MQ=255)
                            aaCGGTTAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGGTAAACGGAat                             >  1:3213155/1‑69 (MQ=255)
                            aaCGGATAACGGAACAATCGCTGCCTGCAATCCAATCGATGTACCGACAATAACGCGGCTAAACGGaa                              >  1:2566647/1‑68 (MQ=255)
                             aCGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAAt                                               <  1:1061729/50‑1 (MQ=255)
                               ggTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAatca                           >  1:1932395/1‑68 (MQ=255)
                                         aCAATCGCTGCCTTCAAACCAATCGATGTACCGACAAt                                               >  1:892570/1‑38 (MQ=255)
                                         aCAATCGCTGCCTGCAAACCAATCGATGTGCCGACAATAACGCGGGTAAACGGAATCATCCATAcaa                  >  1:468270/1‑67 (MQ=255)
                                         aCAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATAcaa                  >  1:3203826/1‑67 (MQ=255)
                                         aCAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATAcaa                  >  1:3333688/1‑67 (MQ=255)
                                         aCAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATAcaa                  >  1:2234253/1‑67 (MQ=255)
                                         aCAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATAca                   >  1:1715366/1‑66 (MQ=255)
                                              cGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGAt           <  1:1319605/69‑1 (MQ=255)
                                              cGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGAt           <  1:2804898/69‑1 (MQ=255)
                                              cGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGAt           <  1:3028680/69‑1 (MQ=255)
                                              cGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGAt           <  1:811864/69‑1 (MQ=255)
                                              cGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGAt           <  1:3076880/69‑1 (MQ=255)
                                               gCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATa          >  1:2502479/1‑69 (MQ=255)
                                               gCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATa          >  1:2111651/1‑69 (MQ=255)
                                               gCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATa          >  1:190748/1‑69 (MQ=255)
                                               gCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATa          >  1:3257404/1‑69 (MQ=255)
                                               gCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATa          >  1:458907/1‑69 (MQ=255)
                                               gCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATa          >  1:1870155/1‑69 (MQ=255)
                                               gCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATa          >  1:1157682/1‑69 (MQ=255)
                                               gCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATa          >  1:155563/1‑69 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAAGGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:846677/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:3324819/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:424226/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:419201/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:1310115/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:3340911/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:2430266/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:3194099/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:2550180/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:1242451/69‑1 (MQ=255)
                                                  gcctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:2915814/69‑1 (MQ=255)
                                                   cctgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATg       <  1:1688350/68‑1 (MQ=255)
                                                     tgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:3406579/41‑1 (MQ=255)
                                                     tgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:1960929/41‑1 (MQ=255)
                                                     tgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:2993368/41‑1 (MQ=255)
                                                     tgcAAACCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:2141722/41‑1 (MQ=255)
                                                       cAAACCAATCGATGTACCGACAATAACGCGGGTAAACgg                                <  1:1089364/39‑1 (MQ=255)
                                                         aaCCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACgg  <  1:781056/67‑1 (MQ=255)
                                                         aaCCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACgg  <  1:1982049/67‑1 (MQ=255)
                                                         aaCCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACgg  <  1:2612001/67‑1 (MQ=255)
                                                         aaCCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACgg  <  1:2785172/67‑1 (MQ=255)
                                                         aaCCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACgg  <  1:3048247/67‑1 (MQ=255)
                                                          |                                                                 
CCATACGGGCAATAAACGGTGCTGCACCAACGGTTAACGGAACAATCGCTGCCTGCAAACCAATCGATGTACCGACAATAACGCGGGTAAACGGAATCATCCATACAAGCAAGATAATGAACGG  >  minE/200125‑200248

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: