Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 904,516 A→T 28.8% S163C (AGC→TGC)  ycjP → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE904,5160AT28.8% 112.3 / 60.7 73S163C (AGC→TGC) ycjPpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base T (19/2);  total (47/26)
Fisher's exact test for biased strand distribution p-value = 2.98e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

ACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTC  >  minE/904448‑904579
                                                                    |                                                                
acacTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGc                                                                <  1:2537710/71‑1 (MQ=255)
 cacTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑t                                                              <  1:201704/71‑1 (MQ=255)
   cTGAGATGGCGCTGATCATCACCATGGTGGCGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TAc                                                            <  1:2607593/71‑1 (MQ=255)
   cTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGc                                                                <  1:1387044/68‑1 (MQ=255)
   cTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TAc                                                            <  1:283462/71‑1 (MQ=255)
      agaTGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTc                                                         <  1:2637896/71‑1 (MQ=255)
      agaTGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTc                                                         <  1:949747/71‑1 (MQ=255)
            cgcTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACAcc                                                   <  1:2340250/71‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCTCCGTGTTTATGCTGAAAAGC‑TACTTCGa                                                       <  1:4138999/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TAc                                                            >  1:500985/1‑59 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TAc                                                            >  1:2544730/1‑59 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TAc                                                            >  1:1323503/1‑59 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TAc                                                            >  1:3289521/1‑59 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGa                                                       <  1:2220845/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGa                                                       <  1:3081884/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGa                                                       <  1:902251/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGa                                                       <  1:490716/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATc                                                >  1:1308196/1‑71 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATc                                                >  1:1793190/1‑71 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATc                                                >  1:3990977/1‑71 (MQ=255)
                 atcatcACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATccc                                              >  1:1524023/1‑71 (MQ=255)
                 atcatcACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATccc                                              >  1:637676/1‑71 (MQ=255)
                   catcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATcccc                                             >  1:132348/1‑70 (MQ=255)
                     tcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATcccc                                             <  1:156936/68‑1 (MQ=255)
                     tcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATcccc                                             <  1:2994407/68‑1 (MQ=255)
                       aCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGAt                                          >  1:2563429/1‑69 (MQ=255)
                       aCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGAt                                          >  1:4156211/1‑69 (MQ=255)
                         cATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGaaaa                                                                  >  1:1652544/1‑44 (MQ=255)
                         cATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGaaaa                                                                  >  1:3333122/1‑44 (MQ=255)
                          aTGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGAt                                          >  1:2743860/1‑66 (MQ=255)
                             gTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCg                                   >  1:2101955/1‑70 (MQ=255)
                              acatgagaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:1703422‑M2/23‑41 (MQ=255)
                              acatgagaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:318038‑M2/23‑41 (MQ=255)
                              acatgagaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:3884928‑M2/23‑41 (MQ=255)
                              acatgagaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:1121767‑M2/23‑41 (MQ=255)
                               catgagaagaacggcatagtgCGTGTTTATGCTTAAATgcg                                                               <  1:2105507‑M2/20‑2 (MQ=255)
                               catgagaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:476092‑M2/22‑40 (MQ=255)
                               catgagaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:1186597‑M2/22‑40 (MQ=255)
                                atgagaagaacggcatagtgCGTGTTTATGCTTAAATGCGTacttat                                                        <  1:1398037‑M2/27‑3 (MQ=255)
                                  gagaagaacggcatagtgCGTGTTTATGCTTAAATgcg                                                               >  1:3816963‑M2/19‑37 (MQ=255)
                                  gagaagaacggcatagtgCGTGTTTATGCTTAAATgcg                                                               >  1:449586‑M2/19‑37 (MQ=255)
                                  gagaagaacggcatagtgCGTGTTTATGCTTAAATgcg                                                               >  1:237312‑M2/19‑37 (MQ=255)
                                  gagaagaacggcatagtgCGTGTTTATGCTTAAATgcg                                                               >  1:1361002‑M2/19‑37 (MQ=255)
                                  gagaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:1476330‑M2/19‑37 (MQ=255)
                                   cAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAgc                            >  1:3407769/1‑71 (MQ=255)
                                   agaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:99047‑M2/18‑36 (MQ=255)
                                   agaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:1636201‑M2/18‑36 (MQ=255)
                                   agaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:119418‑M2/18‑36 (MQ=255)
                                   agaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:3256942‑M2/18‑36 (MQ=255)
                                   agaagaacggcatagtgCGTGTTTATGCTTAAATGc                                                                >  1:3974867‑M2/18‑36 (MQ=255)
                                    aaaCTCTGCCCACCGCCGTGTTTATGCTGAAAAGG‑TACTTCGa                                                       >  1:3130475/1‑43 (MQ=255)
                                    aaaCTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGa                                                       >  1:1710956/1‑43 (MQ=255)
                                    aaaCTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGa                                                       >  1:271354/1‑43 (MQ=255)
                                     aagaacggcatagtgCGTGTTTATGCTTAAATgcg                                                               >  1:806065‑M2/16‑34 (MQ=255)
                                     aagaacggcatagtgCGTGTTTATGCTTAAATgcg                                                               >  1:1454644‑M2/16‑34 (MQ=255)
                                     aagaacggcatagtgCGTGTTTATGCTTAAATgcg                                                               >  1:617631‑M2/16‑34 (MQ=255)
                                      aCTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAgcggc                         <  1:3912401/71‑1 (MQ=255)
                                      aCTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAgcgg                          >  1:3016625/1‑70 (MQ=255)
                                       ctctGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAgcggc                         >  1:3896271/1‑70 (MQ=255)
                                       ctctGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAgcgg                          >  1:1009866/1‑69 (MQ=255)
                                        tctGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACCCCATCCCCGATGAGATCGAAGAAGCCGCga                       >  1:2938343/1‑71 (MQ=255)
                                        tctGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCga                       >  1:332859/1‑71 (MQ=255)
                                        tctGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCga                       >  1:1251521/1‑71 (MQ=255)
                                            cccACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCgaagaa                              <  1:526848/60‑1 (MQ=255)
                                            cccACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCgaagaa                              <  1:478692/60‑1 (MQ=255)
                                            cccACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCgaagaa                              <  1:1118032/60‑1 (MQ=255)
                                                      tgtTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATga                                        <  1:474167/40‑1 (MQ=255)
                                                      tgtTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATga                                        <  1:817656/40‑1 (MQ=255)
                                                        tttATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCaa       >  1:3921925/1‑71 (MQ=255)
                                                        tttATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCaa       >  1:3086150/1‑71 (MQ=255)
                                                        tttATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCa        <  1:2584041/70‑1 (MQ=255)
                                                          tATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCga                       >  1:1155916/1‑53 (MQ=255)
                                                            tGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATgatgg                 <  1:2979859/57‑1 (MQ=255)
                                                            tGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGt   <  1:4030717/71‑1 (MQ=255)
                                                             gCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTc  >  1:3226865/1‑71 (MQ=255)
                                                                  aaaGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCgatg                     >  1:285710/1‑47 (MQ=255)
                                                                  aaaGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCgatg                     >  1:1321726/1‑47 (MQ=255)
                                                                    |                                                                
ACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGC‑TACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTC  >  minE/904448‑904579

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: ‑