Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 602,686 C→T 52.2% S39S (TCC→TCT ycaD → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE602,6860CT52.2% ‑0.3 / 60.5 69S39S (TCC→TCTycaDpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (3/30);  new base T (17/19);  total (20/49)
Fisher's exact test for biased strand distribution p-value = 5.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAC  >  minE/602620‑602751
                                                                  |                                                                 
gttgACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:819065/68‑1 (MQ=255)
gttgACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:1700311/68‑1 (MQ=255)
gttgACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:551492/68‑1 (MQ=255)
gttgACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:1597583/68‑1 (MQ=255)
gttgACTCTGGCGATTGCGGTGCTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:60939/68‑1 (MQ=255)
          gCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAgg                                                         <  1:190209/67‑1 (MQ=255)
          gCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAgg                                                         <  1:3025909/67‑1 (MQ=255)
           cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:1897125/1‑69 (MQ=255)
           cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:1780813/1‑69 (MQ=255)
           cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:137314/1‑69 (MQ=255)
           cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:2576949/1‑69 (MQ=255)
           cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:2394228/1‑69 (MQ=255)
           cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAgg                                                         <  1:3155191/66‑1 (MQ=255)
              ttGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCATGTAGGCg                                                   <  1:1168070/69‑1 (MQ=255)
              ttGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCg                                                   <  1:1454816/69‑1 (MQ=255)
              ttGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCg                                                   <  1:2609432/69‑1 (MQ=255)
              ttGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCg                                                   <  1:2722527/69‑1 (MQ=255)
               tGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCg                                                   <  1:760859/68‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:2552897/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:3013623/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:2948499/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:3014087/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:2571926/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:2760935/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:310347/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1807842/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1759511/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:124087/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:635476/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1603494/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1448261/69‑1 (MQ=255)
                   gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1477355/69‑1 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACATATGTCTACATgg                                                             >  1:2502728/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:1810881/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:287165/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:290306/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:1404711/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:1596747/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:2534738/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:554598/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:362381/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:2055150/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:2053069/1‑47 (MQ=255)
                          aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATg                                                              >  1:118598/1‑46 (MQ=255)
                           tACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATc                                      <  1:862047/69‑1 (MQ=255)
                           tACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATc                                      <  1:1821421/69‑1 (MQ=255)
                            acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:3353072/69‑1 (MQ=255)
                            acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:678920/69‑1 (MQ=255)
                            acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:326535/69‑1 (MQ=255)
                            acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:3255362/69‑1 (MQ=255)
                            acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:3116120/69‑1 (MQ=255)
                            acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:255836/69‑1 (MQ=255)
                            acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:2653621/69‑1 (MQ=255)
                               cTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:424965/68‑1 (MQ=255)
                               cTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:1434836/68‑1 (MQ=255)
                                 cGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTAttt                                <  1:848157/69‑1 (MQ=255)
                                           tGGCTCGCCCAGGATCACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATc                                      <  1:559428/53‑1 (MQ=255)
                                                   ccATGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1756026/69‑1 (MQ=255)
                                                   ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:2166486/69‑1 (MQ=255)
                                                   ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1200046/69‑1 (MQ=255)
                                                   ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1200260/69‑1 (MQ=255)
                                                   ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:2937887/69‑1 (MQ=255)
                                                   ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:2629312/69‑1 (MQ=255)
                                                   ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1653057/69‑1 (MQ=255)
                                                   ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1339968/69‑1 (MQ=255)
                                                    cAGGAACACATGTCCACATGGCTGGTAGGCGTTGTCAGCTCATcc                                     <  1:674842/45‑1 (MQ=255)
                                                              tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:2781375/37‑1 (MQ=255)
                                                              tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:1643644/37‑1 (MQ=255)
                                                              tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:297423/37‑1 (MQ=255)
                                                              tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:2013775/37‑1 (MQ=255)
                                                              tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:2145882/37‑1 (MQ=255)
                                                               gTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAc  >  1:374292/1‑69 (MQ=255)
                                                               gTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAc  >  1:2775003/1‑69 (MQ=255)
                                                               gTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAc  >  1:1775879/1‑69 (MQ=255)
                                                                  |                                                                 
GTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAC  >  minE/602620‑602751

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