Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 199,631 G→T 7.6% A41A (GCC→GCA metQ ← DL‑methionine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE199,6310GT 7.6% 247.6 / 4.8 79A41A (GCC→GCAmetQDL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (48/24);  new base T (4/2);  total (52/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCA  >  minE/199564‑199695
                                                                   |                                                                
ccAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGcg                                                               >  1:4296210/1‑70 (MQ=255)
 cAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGcgc                                                              <  1:4101298/71‑3 (MQ=255)
  agttCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGcc                                                             <  1:1546961/68‑1 (MQ=255)
      cAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACcaac                                                          >  1:3109691/1‑70 (MQ=255)
      cAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACcaac                                                          >  1:4920648/1‑70 (MQ=255)
      cAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCcat                                                           >  1:661592/1‑68 (MQ=255)
      cAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCcat                                                           >  1:2335649/1‑68 (MQ=255)
       aaCGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACcaac                                                          <  1:553301/69‑1 (MQ=255)
        aCGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTACTGCAACCTGCTGTTCGGCACCAACAAt                                                       <  1:2688514/71‑1 (MQ=255)
           tCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCAc                                                    <  1:3783933/71‑1 (MQ=255)
           tCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCAc                                                    >  1:3192328/1‑71 (MQ=255)
             cAGGCCTTATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGc                                                  >  1:732471/1‑71 (MQ=255)
              aGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGcc                                                 <  1:946684/71‑1 (MQ=255)
               ggCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCa                                                     >  1:2362217/1‑66 (MQ=255)
                 ccATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACg                                                   >  1:4402586/1‑66 (MQ=255)
                  cATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCgg                                                                 >  1:2581423/1‑51 (MQ=255)
                  cATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCgg                                                                 >  1:2293087/1‑51 (MQ=255)
                  cATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACt                                             >  1:2912895/1‑71 (MQ=255)
                  cATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACt                                             >  1:4797009/1‑71 (MQ=255)
                   atatTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGAc                                              >  1:1435920/1‑69 (MQ=255)
                   atatTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCAGCACCAACAATCACGCCGACtt                                            <  1:5226942/71‑1 (MQ=255)
                    tatTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACttt                                           <  1:4599779/71‑1 (MQ=255)
                     atTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACttt                                           <  1:4870568/70‑1 (MQ=255)
                      tttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGc                                                                >  1:3245063/1‑48 (MQ=255)
                      tttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCa                                                               >  1:3503704/1‑49 (MQ=255)
                      tttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCACCAATc                                                      >  1:1993699/1‑58 (MQ=255)
                      tttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTaa                                         >  1:4701123/1‑71 (MQ=255)
                      tttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTa                                          <  1:4414412/70‑1 (MQ=255)
                       ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAAt                                        >  1:5197455/1‑71 (MQ=255)
                       ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAAt                                        >  1:3959560/1‑71 (MQ=255)
                       ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAAt                                        >  1:4436240/1‑71 (MQ=255)
                       ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAAt                                        >  1:2851138/1‑71 (MQ=255)
                        tGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAAtg                                       >  1:3717946/1‑71 (MQ=255)
                        tGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAAtg                                       >  1:3724639/1‑71 (MQ=255)
                          tCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGa                                               >  1:3354070/1‑61 (MQ=255)
                          tCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAAtgtg                                     <  1:2101692/71‑1 (MQ=255)
                           cTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTgg                                    >  1:2120852/1‑71 (MQ=255)
                            tttCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATc                                                      >  1:291920/1‑52 (MQ=255)
                              tCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTgg                                    >  1:5155188/1‑68 (MQ=255)
                              tCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGtt                                  <  1:345309/70‑1 (MQ=255)
                               cGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTg                                >  1:121920/1‑71 (MQ=255)
                                 cAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACcaac                                                          <  1:1999374/43‑1 (MQ=255)
                                     tttCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCtt                          >  1:1769013/1‑71 (MQ=255)
                                       tCTGCGCGCCTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCtttt                        >  1:5476276/1‑71 (MQ=255)
                                       tCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACg                                                   >  1:1042255/1‑44 (MQ=255)
                                       tCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCtttt                        >  1:5141071/1‑71 (MQ=255)
                                        cTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttttt                       >  1:5647715/1‑71 (MQ=255)
                                        cTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCtttt                        <  1:3818655/70‑1 (MQ=255)
                                          gcgcGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTg                                <  1:1765505/60‑1 (MQ=255)
                                           cgcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttttt                       <  1:4991990/68‑1 (MQ=255)
                                           cgcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCa                     <  1:5305873/70‑1 (MQ=255)
                                           cgcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCa                     <  1:4383731/70‑1 (MQ=255)
                                           cgcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCa                     <  1:4940091/70‑1 (MQ=255)
                                            gcgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCg                                                >  1:3226176/1‑42 (MQ=255)
                                             cgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATcc                  >  1:2961040/1‑71 (MQ=255)
                                             cgACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATcc                  >  1:476340/1‑71 (MQ=255)
                                              gACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCt                 <  1:2163639/71‑1 (MQ=255)
                                              gACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCt                 >  1:2461757/1‑71 (MQ=255)
                                              gACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATcgct                         >  1:2968153/1‑60 (MQ=255)
                                                  tCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGAcc             >  1:2699028/1‑71 (MQ=255)
                                                  tCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGAc              >  1:2299905/1‑70 (MQ=255)
                                                   cTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCg            >  1:2364959/1‑71 (MQ=255)
                                                   cTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCg            >  1:5150116/1‑71 (MQ=255)
                                                    tGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCg            >  1:3483311/1‑70 (MQ=255)
                                                    tGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCg            >  1:450319/1‑70 (MQ=255)
                                                    tGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCg            >  1:697501/1‑70 (MQ=255)
                                                     gCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttt                         <  1:1587460/56‑1 (MQ=255)
                                                     gCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttt                         <  1:62238/56‑1 (MQ=255)
                                                     gCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttt                         <  1:2766342/56‑1 (MQ=255)
                                                       agccTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTc                      <  1:3873994/55‑1 (MQ=255)
                                                        aCCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTg                >  1:1013324/1‑62 (MQ=255)
                                                         ccTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCAt                    >  1:4334295/1‑57 (MQ=255)
                                                             ctgTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTAcca  <  1:4673989/71‑1 (MQ=255)
                                                                ttCGGCACCAACAATCACGCCGACTTTAATGTGGt                                   <  1:1829858/35‑1 (MQ=255)
                                                                ttCGGCACCAACAATCACGCCGACTTTAATGTGGt                                   <  1:1744081/35‑1 (MQ=255)
                                                                  cGGCACCAACAATCACGCCGACTTTAATGTGGTTTgg                               >  1:1157424/1‑37 (MQ=255)
                                                                  cGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGa                              >  1:1412544/1‑38 (MQ=255)
                                                                  cGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGa                              >  1:3885790/1‑38 (MQ=255)
                                                                  cGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttt                         >  1:218466/1‑43 (MQ=255)
                                                                  cGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttt                         >  1:5657249/1‑43 (MQ=255)
                                                                  cGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCttt                         >  1:2723486/1‑43 (MQ=255)
                                                                   |                                                                
CCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTTTTTCATCCTGACCGCAGCCTACCA  >  minE/199564‑199695

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: