Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 474,837 C→T E157K (GAA→AAA)  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,8370CT95.9% 174.5 / ‑2.5 49E157K (GAA→AAA) ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (28/19);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAG  >  REL606/474764‑474909
                                                                         |                                                                        
cccACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTca                                                                        >  1:498411/1‑76 (MQ=255)
  cacaTCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTcaca                                                                      >  2:709147/1‑76 (MQ=255)
  cacaTCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTcaca                                                                      >  1:1692201/1‑76 (MQ=255)
   acaTCTCCCTTTTCCATCTTCCCTGCCACTCCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACaa                                                                     <  2:605202/76‑1 (MQ=255)
   acaTCGCCCTTTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACaa                                                                     <  2:1273839/76‑1 (MQ=255)
    caTCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAAt                                                                    >  2:551396/1‑76 (MQ=255)
         cccTGTCCCACCATCCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATcaacc                                                               <  1:661762/76‑1 (MQ=255)
          ccTGTTCCATCATTCTTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATcaacca                                                              <  1:432141/76‑1 (MQ=255)
          ccTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATcaacca                                                              >  2:1503972/1‑76 (MQ=255)
           cTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATcaaccaa                                                             >  2:1698272/1‑76 (MQ=255)
            tGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATcaaccaac                                                            >  1:1800032/1‑76 (MQ=255)
              ttCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCg                                                          >  1:1537480/1‑76 (MQ=255)
                ccATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGAt                                                        >  2:1039837/1‑76 (MQ=255)
                 catcatTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATg                                                       >  2:372507/1‑76 (MQ=255)
                   tgatTCCTGCCACTGGGGGCAGCAACACCAGCGTCAGAGCCATTGCCAGGTCTTTCACAATCAGCCAACCGATGGc                                                     <  1:23043/74‑1 (MQ=255)
                   tcatTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGc                                                     >  1:1610303/1‑76 (MQ=255)
                    catTCCTGCCACTGCGGGCAGCAACACCAGCGTCGGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCg                                                    >  1:1773223/1‑76 (MQ=255)
                    catTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTGCCAGGTCTTTCACAATCAACCAACCGATGGCg                                                    <  1:228115/76‑1 (MQ=255)
                    catTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCg                                                    >  2:1210859/1‑76 (MQ=255)
                     atTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGa                                                   >  1:1877038/1‑76 (MQ=255)
                     atTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGa                                                   >  2:524974/1‑76 (MQ=255)
                       cccTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGAtt                                                 <  1:763079/75‑1 (MQ=255)
                       cccTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGAtt                                                 <  1:937878/75‑1 (MQ=255)
                        ccTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGAttt                                                <  2:1492597/76‑1 (MQ=255)
                         ctcccACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTg                                               <  2:877638/73‑1 (MQ=255)
                          tGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGc                                              >  2:232641/1‑76 (MQ=255)
                           gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGcc                                             >  2:1597885/1‑76 (MQ=255)
                           gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGcc                                             >  2:376645/1‑76 (MQ=255)
                           gCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGcc                                             >  2:668628/1‑76 (MQ=255)
                                   gggCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGa                                     >  1:367746/1‑76 (MQ=255)
                                      cagcaACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTg                                  >  1:197320/1‑76 (MQ=255)
                                         caACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCa                               >  2:183233/1‑76 (MQ=255)
                                         caACACCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCa                               >  1:843823/1‑76 (MQ=255)
                                           acacCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaat                             <  1:182064/76‑1 (MQ=255)
                                             cccAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaatta                           <  1:1148009/75‑1 (MQ=255)
                                             acCAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaatta                           <  2:1087875/76‑1 (MQ=255)
                                              ccAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaattaa                          <  2:942495/76‑1 (MQ=255)
                                               cAGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaattaat                         <  1:1643037/76‑1 (MQ=255)
                                                aGCGTCAGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCaattaatt                        <  2:860980/76‑1 (MQ=255)
                                                      aGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGtt                  >  2:1517054/1‑76 (MQ=255)
                                                      aGAACCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGtt                  >  1:1599146/1‑76 (MQ=255)
                                                         aCCATTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGttctt               <  2:1507046/76‑1 (MQ=255)
                                                            aTTACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCaa            >  1:1462790/1‑76 (MQ=255)
                                                            aTTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCaa            >  1:279742/1‑76 (MQ=255)
                                                             ttACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAg           >  2:280565/1‑76 (MQ=255)
                                                              tACCAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGt          >  1:1850437/1‑76 (MQ=255)
                                                                ccAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTgc        >  1:987213/1‑76 (MQ=255)
                                                                ccAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTgc        <  1:709147/76‑1 (MQ=255)
                                                                 cAGGTCTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTgcg       <  2:340892/76‑1 (MQ=255)
                                                                      cTTTCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAg  <  1:1125347/76‑1 (MQ=255)
                                                                         |                                                                        
CCCACATCGCCCTGTTCCATCATTCCTGCCACTGCGGGCAGCAACACCAGCGTCAGAACCATTACCAGGTCTTCCACAATCAACCAACCGATGGCGATTTGCCCACGCTGACTGTCAATTAATTGCCGTTCTTCAAGTGCGCGCAG  >  REL606/474764‑474909

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