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

ACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTC  >  minE/904448‑904579
                                                                    |                                                                
acacTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGc                                                                <  1:2537710/71‑1 (MQ=255)
 cacTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCt                                                              <  1:201704/71‑1 (MQ=255)
   cTGAGATGGCGCTGATCATCACCATGGTGGCGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTAc                                                            <  1:2607593/71‑1 (MQ=255)
   cTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGc                                                                <  1:1387044/68‑1 (MQ=255)
   cTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTAc                                                            <  1:283462/71‑1 (MQ=255)
      agaTGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTc                                                         <  1:2637896/71‑1 (MQ=255)
      agaTGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTc                                                         <  1:949747/71‑1 (MQ=255)
            cgcTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACAcc                                                   <  1:2340250/71‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCTCCGTGTTTATGCTGAAAAGCTACTTCGa                                                       <  1:4138999/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTAc                                                            >  1:500985/1‑59 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTAc                                                            >  1:2544730/1‑59 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTAc                                                            >  1:1323503/1‑59 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTAc                                                            >  1:3289521/1‑59 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGa                                                       <  1:2220845/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGa                                                       <  1:3081884/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGa                                                       <  1:902251/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGa                                                       <  1:490716/64‑1 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATc                                                >  1:1308196/1‑71 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATc                                                >  1:1793190/1‑71 (MQ=255)
               tGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATc                                                >  1:3990977/1‑71 (MQ=255)
                 atcatcACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATccc                                              >  1:1524023/1‑71 (MQ=255)
                 atcatcACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATccc                                              >  1:637676/1‑71 (MQ=255)
                   catcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATcccc                                             >  1:132348/1‑70 (MQ=255)
                     tcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATcccc                                             <  1:156936/68‑1 (MQ=255)
                     tcaCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATcccc                                             <  1:2994407/68‑1 (MQ=255)
                       aCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGAt                                          >  1:2563429/1‑69 (MQ=255)
                       aCCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGAt                                          >  1:4156211/1‑69 (MQ=255)
                         cATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGaaaa                                                                  >  1:1652544/1‑44 (MQ=255)
                         cATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGaaaa                                                                  >  1:3333122/1‑44 (MQ=255)
                          aTGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGAt                                          >  1:2743860/1‑66 (MQ=255)
                             gTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCg                                   >  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)
                                   cAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAgc                            >  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)
                                    aaaCTCTGCCCACCGCCGTGTTTATGCTGAAAAGGTACTTCGa                                                       >  1:3130475/1‑43 (MQ=255)
                                    aaaCTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGa                                                       >  1:1710956/1‑43 (MQ=255)
                                    aaaCTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGa                                                       >  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)
                                      aCTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAgcggc                         <  1:3912401/71‑1 (MQ=255)
                                      aCTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAgcgg                          >  1:3016625/1‑70 (MQ=255)
                                       ctctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAgcggc                         >  1:3896271/1‑70 (MQ=255)
                                       ctctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAgcgg                          >  1:1009866/1‑69 (MQ=255)
                                        tctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACCCCATCCCCGATGAGATCGAAGAAGCCGCga                       >  1:2938343/1‑71 (MQ=255)
                                        tctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCga                       >  1:332859/1‑71 (MQ=255)
                                        tctGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCga                       >  1:1251521/1‑71 (MQ=255)
                                            cccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCgaagaa                              <  1:526848/60‑1 (MQ=255)
                                            cccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCgaagaa                              <  1:478692/60‑1 (MQ=255)
                                            cccACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCgaagaa                              <  1:1118032/60‑1 (MQ=255)
                                                      tgtTTATGCTGAAAAGCTACTTCGACACCATCCCCGATga                                        <  1:474167/40‑1 (MQ=255)
                                                      tgtTTATGCTGAAAAGCTACTTCGACACCATCCCCGATga                                        <  1:817656/40‑1 (MQ=255)
                                                        tttATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCaa       >  1:3921925/1‑71 (MQ=255)
                                                        tttATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCaa       >  1:3086150/1‑71 (MQ=255)
                                                        tttATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCa        <  1:2584041/70‑1 (MQ=255)
                                                          tATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCga                       >  1:1155916/1‑53 (MQ=255)
                                                            tGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATgatgg                 <  1:2979859/57‑1 (MQ=255)
                                                            tGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGt   <  1:4030717/71‑1 (MQ=255)
                                                             gCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTc  >  1:3226865/1‑71 (MQ=255)
                                                                  aaaGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCgatg                     >  1:285710/1‑47 (MQ=255)
                                                                  aaaGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCgatg                     >  1:1321726/1‑47 (MQ=255)
                                                                    |                                                                
ACACTGAGATGGCGCTGATCATCACCATGGTGACGCAAACTCTGCCCACCGCCGTGTTTATGCTGAAAAGCTACTTCGACACCATCCCCGATGAGATCGAAGAAGCGGCGATGATGGATGGCCTCAACCGTC  >  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: