Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 199,601 T→C 7.4% K51K (AAA→AAG metQ ← DL‑methionine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE199,6010TC 7.4% 229.5 / 6.7 82K51K (AAA→AAGmetQDL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (45/30);  new base C (4/2);  total (50/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CGTTTGGCAGAACATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTT  >  minE/199534‑199671
                                                                   |                                                                      
cGTTTGGCAGAACATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACtt                                                                       <  1:2010151/69‑1 (MQ=255)
    tGGCAGAACATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTgcg                                                                 >  1:4789289/1‑71 (MQ=255)
    tGGCAGAACATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTgcg                                                                 >  1:2401569/1‑71 (MQ=255)
      gCAGAACATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTgcgc                                                                <  1:3550483/70‑1 (MQ=255)
        aGAACATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCTCAACTTTCTGCGCGAc                                                             <  1:2215594/71‑1 (MQ=255)
            cATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACt                                                                        >  1:5668060/1‑56 (MQ=255)
            cATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACt                                                                        >  1:4177215/1‑56 (MQ=255)
            cATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACt                                                                        >  1:5157729/1‑56 (MQ=255)
            cATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCt                                                         >  1:4673807/1‑71 (MQ=255)
             aTAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCt                                                         <  1:1337058/70‑1 (MQ=255)
               aGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACtt                                                                       >  1:358721/1‑54 (MQ=255)
               aGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACtt                                                                       >  1:2920160/1‑54 (MQ=255)
               aGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCt                                                                    >  1:2191452/1‑57 (MQ=255)
               aGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCt                                                                    >  1:4783052/1‑57 (MQ=255)
                gTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCa                                                      <  1:2439409/70‑1 (MQ=255)
                gTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCa                                                      <  1:4693274/70‑1 (MQ=255)
                gTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCa                                                      <  1:3098313/70‑1 (MQ=255)
                   gTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTc                                                                     >  1:1118941/1‑52 (MQ=255)
                   gTTGAAGGTTACCAGCTCAACGTCCAGGCCATATATGTCTTTCGCAACTTTc                                                                     >  1:135150/1‑52 (MQ=255)
                    ttGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCaa                                                     <  1:3199998/67‑1 (MQ=255)
                      gAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACttt                                                                      >  1:4082564/1‑48 (MQ=255)
                       aaGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACctgctg                                              >  1:498431/1‑71 (MQ=255)
                       aaGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACctgctg                                              >  1:3125641/1‑71 (MQ=255)
                          gTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTc                                           <  1:237018/71‑1 (MQ=255)
                          cctACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTc                                           >  1:44007/3‑71 (MQ=255)
                           ttACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTgc                                                                  <  1:1820407/47‑1 (MQ=255)
                           ttACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTgc                                                                  <  1:444054/47‑1 (MQ=255)
                             aCCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTgcgc                                                                >  1:3396982/1‑47 (MQ=255)
                             aCCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTgcgc                                                                >  1:2015192/1‑47 (MQ=255)
                             aCCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACctg                                                 >  1:1912610/1‑62 (MQ=255)
                              ccAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGcg                                       >  1:4296210/1‑70 (MQ=255)
                               cAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGcgc                                      <  1:4101298/71‑3 (MQ=255)
                               cAGCTCAACGTCCAGGCCATATATCACTTTCGCAACATTCTGCCCGACTTCTGCAACCTGCTGtt                                            >  1:4735878/1‑65 (MQ=255)
                                agttCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGcc                                     <  1:1546961/68‑1 (MQ=255)
                                    cAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACt                                                            >  1:4723275/1‑44 (MQ=255)
                                    cAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACcaac                                  >  1:3109691/1‑70 (MQ=255)
                                    cAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACcaac                                  >  1:4920648/1‑70 (MQ=255)
                                    cAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCcat                                   >  1:2335649/1‑68 (MQ=255)
                                    cAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCcat                                   >  1:661592/1‑68 (MQ=255)
                                     aaCGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACcaac                                  <  1:553301/69‑1 (MQ=255)
                                      aCGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTACTGCAACCTGCTGTTCGGCACCAACAAt                               <  1:2688514/71‑1 (MQ=255)
                                         tCCAGGCCATATTTGTCTTTCGCAACTTTCTTCGCGAc                                                             <  1:2473881/38‑1 (MQ=255)
                                         tCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGAc                                                             <  1:3847435/38‑1 (MQ=255)
                                         tCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGAc                                                             <  1:989113/38‑1 (MQ=255)
                                         tCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCAc                            <  1:3783933/71‑1 (MQ=255)
                                         tCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCAc                            >  1:3192328/1‑71 (MQ=255)
                                           cAGGCCTTATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGc                          >  1:732471/1‑71 (MQ=255)
                                           cAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTc                                           <  1:5571799/54‑1 (MQ=255)
                                           cAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTc                                           <  1:2996453/54‑1 (MQ=255)
                                            aGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGcc                         <  1:946684/71‑1 (MQ=255)
                                             ggCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAAc                                                    >  1:624659/1‑43 (MQ=255)
                                             ggCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAAc                                                    >  1:4501396/1‑43 (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)
                                                    tttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAAc                                                    >  1:5120446/1‑36 (MQ=255)
                                                    tttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAAc                                                    >  1:548530/1‑36 (MQ=255)
                                                    tttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAAc                                                    >  1:2795265/1‑36 (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)
                                                     ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAAc                                                    >  1:1296223/1‑35 (MQ=255)
                                                     ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACctgctg                                              <  1:2149686/41‑1 (MQ=255)
                                                     ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACctgctg                                              <  1:904095/41‑1 (MQ=255)
                                                     ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACctgctg                                              <  1:5334055/41‑1 (MQ=255)
                                                     ttGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACctgctg                                              <  1:4268615/41‑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:3724639/1‑71 (MQ=255)
                                                      tGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAAtg               >  1:3717946/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)
                                                                   |                                                                      
CGTTTGGCAGAACATAGTCGTTGAAGGTTACCAGCTCAACGTCCAGGCCATATTTGTCTTTCGCAACTTTCTGCGCGACTTCTGCAACCTGCTGTTCGGCACCAACAATCACGCCGACTTTAATGTGGTTTGGATCTT  >  minE/199534‑199671

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: