Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 474,270 C→T E346K (GAG→AAG)  ybaL ← predicted transporter with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606474,2700CT92.9% 88.6 / ‑0.5 29E346K (GAG→AAG) ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (15/11);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

AGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGA  >  REL606/474197‑474341
                                                                         |                                                                       
aGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTAc                                                                       >  2:54871/1‑76 (MQ=255)
    ccAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAAt                                                                   <  1:823600/76‑1 (MQ=255)
        gTTTTGTCCGTCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCCGGAAAAACGCGCACTAACCAATCTgc                                                               >  1:234387/1‑76 (MQ=255)
        gTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTgc                                                               >  1:825754/1‑76 (MQ=255)
        gTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTgc                                                               <  1:4315/76‑1 (MQ=255)
            tGTCCGGCCGGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCa                                                           <  2:713137/76‑1 (MQ=255)
             gTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAg                                                          >  1:883024/1‑76 (MQ=255)
                  gCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTgg                                                     <  2:909496/76‑1 (MQ=255)
                  gCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTgg                                                     >  2:191915/1‑76 (MQ=255)
                      ccGGCATGAAATTCAATTCCATTCCCAGCCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTggcggc                                                 <  2:671408/75‑1 (MQ=255)
                         tcAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGAt                                              >  1:663446/2‑76 (MQ=255)
                               aaTTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTaa                                        <  2:1027037/76‑1 (MQ=255)
                                aTTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGTTGGTTAAt                                       >  2:447177/1‑76 (MQ=255)
                                     aTGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGGCga                                  >  2:457638/1‑75 (MQ=255)
                                     aTGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGt                                  >  2:56598/1‑76 (MQ=255)
                                        ccATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACg                               <  1:1178921/76‑1 (MQ=255)
                                             cccAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTggg                          <  1:702166/76‑1 (MQ=255)
                                             cccAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTggg                          >  2:314188/1‑76 (MQ=255)
                                             cccAGTCCCGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGGTAATGCCGTACGTTggg                          >  1:372048/1‑76 (MQ=255)
                                                    ccGCCAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGAcc                   >  2:163277/1‑76 (MQ=255)
                                                        cAGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAAc               <  1:1043438/76‑1 (MQ=255)
                                                         aGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACa              <  1:1196251/76‑1 (MQ=255)
                                                         aGGATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACa              <  1:251190/76‑1 (MQ=255)
                                                          ggATAAACGCGAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAg             <  1:1186749/76‑1 (MQ=255)
                                                                 cgcgAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCAcc      >  2:729270/1‑76 (MQ=255)
                                                                    naaCTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGGGACCAAACAGTCGCACCAgg   >  1:462509/2‑76 (MQ=255)
                                                                    gAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAgg   >  1:918011/1‑76 (MQ=255)
                                                                    gAACTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAgg   <  1:528189/76‑1 (MQ=255)
                                                                     aaCTTACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGa  >  1:725760/1‑76 (MQ=255)
                                                                         |                                                                       
AGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGA  >  REL606/474197‑474341

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: