Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 199,135 T→A 33.8% S207C (AGC→TGC)  metQ ← DL‑methionine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE199,1350TA33.8% 73.7 / 66.0 68S207C (AGC→TGC) metQDL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (17/28);  new base A (13/10);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

GTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCAGTTGCG  >  minE/199067‑199200
                                                                    |                                                                 
gTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctgg                                                                 <  1:3744343/71‑1 (MQ=255)
gTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctgg                                                                 <  1:3801964/71‑1 (MQ=255)
gTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctgg                                                                 <  1:2955412/71‑1 (MQ=255)
  ttACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCa                                                               >  1:3566682/1‑71 (MQ=255)
   tacgtacgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCa                                                               >  1:3267759/1‑70 (MQ=255)
   tacgtacgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTAAGGCCAATCTGGCTGGCa                                                               >  1:2275250/1‑70 (MQ=255)
     cgtacgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATa                                                             <  1:530381/70‑1 (MQ=255)
     cgtacgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATa                                                             <  1:4011808/70‑1 (MQ=255)
         cgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAggtgg                                                        <  1:3775403/71‑1 (MQ=255)
         cgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCAt                                                              >  1:3031980/1‑65 (MQ=255)
          ggggACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAggtggt                                                       <  1:2294344/71‑1 (MQ=255)
          ggggACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAggtggt                                                       <  1:2972531/71‑1 (MQ=255)
            ggACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGTAGTCAGGCCAATCTGGCAGGCATa                                                             <  1:2928344/63‑1 (MQ=255)
            ggACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATa                                                             <  1:2063653/63‑1 (MQ=255)
             gACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCAt                                                              <  1:2447980/61‑1 (MQ=255)
             gACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCAt                                                              <  1:569438/61‑1 (MQ=255)
             gACCCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCa                                                               >  1:3058600/1‑60 (MQ=255)
              aCTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAgg                                                           >  1:3630283/1‑63 (MQ=255)
              aCTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAgg                                                           >  1:1922939/1‑63 (MQ=255)
                 ctTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCAt                                                              <  1:187568/57‑1 (MQ=255)
                 ctTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCAt                                                              <  1:2069553/57‑1 (MQ=255)
                 ctTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCAt                                                              <  1:261707/57‑1 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctg                                                                  >  1:3943519/1‑52 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctg                                                                  >  1:1419900/1‑52 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctg                                                                  >  1:2601410/1‑52 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggct                                                                   <  1:2434338/51‑1 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGAt                                                 <  1:2959256/69‑1 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGAt                                                 <  1:1416274/69‑1 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGAt                                                 <  1:348024/69‑1 (MQ=255)
                    tATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACt                                             >  1:1099104/1‑71 (MQ=255)
                    tATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACt                                             >  1:3733224/1‑71 (MQ=255)
                    tATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACt                                             >  1:1640102/1‑71 (MQ=255)
                     aTCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATa                                                <  1:1576679/67‑1 (MQ=255)
                     aTCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTg                                            >  1:477322/1‑71 (MQ=255)
                     aTCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTg                                            >  1:3484109/1‑71 (MQ=255)
                     aTCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTg                                            >  1:3486341/1‑71 (MQ=255)
                      tCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTggcaggca                                                               >  1:1310819/1‑51 (MQ=255)
                        ttCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACAGc                                           >  1:3300623/1‑69 (MQ=255)
                                 taTACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCAt                                                              <  1:1506954/40‑1 (MQ=255)
                                 gATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCAt                                                              <  1:1367774/41‑1 (MQ=255)
                                 gATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCAt                                                              <  1:2210922/41‑1 (MQ=255)
                                 gATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCAt                                                              <  1:2765527/41‑1 (MQ=255)
                                 gATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCAt                                                              <  1:2339535/41‑1 (MQ=255)
                                 gATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCAt                                                              <  1:2917865/41‑1 (MQ=255)
                                   tACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATa                                                             >  1:3662187/1‑40 (MQ=255)
                                   tACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAg                                                            >  1:2459336/1‑41 (MQ=255)
                                   tACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAg                                                            >  1:2902255/1‑41 (MQ=255)
                                   tACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAg                                                            >  1:3510939/1‑41 (MQ=255)
                                   tACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAg                                                            >  1:663136/1‑41 (MQ=255)
                                   tACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAg                                                            >  1:104683/1‑41 (MQ=255)
                                    aCCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCag                                        <  1:1219190/60‑1 (MQ=255)
                                    aCCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCag                                        <  1:2149634/60‑1 (MQ=255)
                                    aCCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCag                                        <  1:337033/60‑1 (MQ=255)
                                     ccGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTgcgc                            <  1:1762485/71‑1 (MQ=255)
                                      cGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCagag                                      >  1:2899734/1‑60 (MQ=255)
                                      cGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCagag                                      >  1:2037236/1‑60 (MQ=255)
                                      cGTCTTTCGCCGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCagag                                      >  1:2162621/1‑59 (MQ=255)
                                          tttCGCCGGAGTCAGGCCAATCTGGCTGGCATAggtggt                                                       <  1:1311553/39‑1 (MQ=255)
                                          tttCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCa                                         >  1:2244827/1‑53 (MQ=255)
                                          tttCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGcgtcgt                       >  1:3050286/1‑71 (MQ=255)
                                          tttCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGcgtcgt                       >  1:278506/1‑71 (MQ=255)
                                          tttCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGcgtcgt                       >  1:814759/1‑71 (MQ=255)
                                          tttCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGcgtcgt                       >  1:1758979/1‑71 (MQ=255)
                                              gCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCa                    <  1:1445629/70‑1 (MQ=255)
                                              gCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCa                    <  1:1250369/70‑1 (MQ=255)
                                                cGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCagag                 <  1:3048637/71‑1 (MQ=255)
                                                    gTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAAcg             >  1:4081806/1‑71 (MQ=255)
                                                     tCAGGCCTATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCg                                    <  1:594350/47‑1 (MQ=255)
                                                     tCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCg                                    <  1:665986/47‑1 (MQ=255)
                                                        ggCCAATCTGGCTGGCATAGGTTGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGc         <  1:1149339/71‑1 (MQ=255)
                                                        ggCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGc         <  1:999962/71‑1 (MQ=255)
                                                        ggCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGc         <  1:1437369/71‑1 (MQ=255)
                                                        ggCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGc         <  1:2592113/71‑1 (MQ=255)
                                                            aaTCTGGCAGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGAGTCGTCCAGAGAACGCGGCAGtt     >  1:3315777/1‑71 (MQ=255)
                                                             aTCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGc         <  1:1441010/66‑1 (MQ=255)
                                                             aTCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCAGTTg    <  1:3750556/71‑1 (MQ=255)
                                                               ctggctggCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCAGTTGCg  >  1:556940/1‑71 (MQ=255)
                                                                    |                                                                 
GTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCAGTTGCG  >  minE/199067‑199200

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: