Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 474,656 Δ1 bp coding (650/1677 nt) 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,6560C.94.6% 149.3 / 5.0 38coding (650/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (10/25);  minor base G (0/2);  total (11/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CAACCCCTAACGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGCGCTGCGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGAT  >  REL606/474591‑474729
                                                                 |                                                                         
cAACCCCTAACGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATcc                                                               <  2:1177415/76‑1 (MQ=255)
 aaCCCCTAACGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCa                                                              <  1:1145046/76‑1 (MQ=255)
 aaCCCCTAACGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCa                                                              <  1:1451911/76‑1 (MQ=255)
      cTAACGCCAGCGCCGGCGCCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCa                                                         <  1:1011849/76‑1 (MQ=255)
       tAACGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCAc                                                        >  2:108722/1‑76 (MQ=255)
       tAACGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCAGTCTGGTGCCATAATCCACGGCAc                                                        >  2:1271387/1‑76 (MQ=255)
        aaCGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCAcc                                                       <  2:441708/76‑1 (MQ=255)
        aaCGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCAcc                                                       <  2:1277721/76‑1 (MQ=255)
         aCGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCAccc                                                      >  2:94988/1‑75 (MQ=255)
           gccagcgccagcACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGa                                                    <  1:2980/76‑1 (MQ=255)
            ccagcgccagcACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGAc                                                   <  2:1724323/76‑1 (MQ=255)
            ccagcgccagcACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGAc                                                   >  2:1670531/1‑76 (MQ=255)
             cagcgccagcACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACg                                                  <  1:1090611/76‑1 (MQ=255)
                cgccagcACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCg                                               <  2:1100226/76‑1 (MQ=255)
                cgccagcACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAAGCCACGGCACCAGACGGCg                                               <  2:1485549/76‑1 (MQ=255)
                     ccACCGCCAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTa                                          <  1:1063464/75‑1 (MQ=255)
                       aCCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTAcc                                        <  2:1400290/76‑1 (MQ=255)
                        ccGACGGGGTAAACAGCTCGCGAAGACCGGTGGCCGTCGGGTGCCATAACCCACGGCACCGGACGGCGACCTACCa                                       <  2:573834/76‑1 (MQ=255)
                        ccGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCa                                       <  2:1200440/76‑1 (MQ=255)
                        ccGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCa                                       <  2:652588/76‑1 (MQ=255)
                           aCAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCATGGCACCAGACGGCGACCTACCAGca                                    >  2:1402716/1‑76 (MQ=255)
                           aCAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGca                                    <  2:91279/76‑1 (MQ=255)
                           aCAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGca                                    <  2:652847/76‑1 (MQ=255)
                            cAGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGGGCCATAATCCACGGCACCAGACGGCGACCTACCAGcat                                   >  2:665252/1‑76 (MQ=255)
                             aGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGcatc                                  <  2:1147626/76‑1 (MQ=255)
                             aGGGTAAACAGCTCGCGAGAACCGGTTGCCGTCTGGGGCCATAATCCAAGGCACCAGACCGCGGACTACCAAcatc                                  >  2:702457/1‑76 (MQ=255)
                                 tAAACAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATaa                              >  2:1708694/1‑76 (MQ=255)
                                   aaCAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATg                            <  1:1307190/76‑1 (MQ=255)
                                     cggCTCGCGAGACCCGGTTGCCGTCGGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGc                          <  1:73623/74‑1 (MQ=255)
                                     cAGCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGc                          >  2:1510096/1‑76 (MQ=255)
                                       gCTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGa                        <  1:81577/76‑1 (MQ=255)
                                        cTCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGAt                       <  1:1660510/76‑1 (MQ=255)
                                         tCGCGAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATa                      <  1:1520726/76‑1 (MQ=255)
                                         gcgcgAAAACCGGTTGCCGTCTGGGGCCATAATCCACGGCCCCAGACGGCGACCTACCAGCATCATAATGGCGATa                      <  2:882217/75‑1 (MQ=255)
                                          cgcgAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATaa                     >  2:6488/1‑76 (MQ=255)
                                          cgcgAGAACCGGTTGCCGTCTTGTGCCCTAAACCACCGCCCCCGACGGCGACCCACCCGCATCATAAAAGCGATaa                     >  2:1553654/1‑76 (MQ=255)
                                           gcgAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATaaa                    <  1:524767/76‑1 (MQ=255)
                                           gcgAGAACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATaaa                    <  1:1221071/76‑1 (MQ=255)
                                             gagaACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCCGCATCACAAAGGCGATaacaa                  >  2:1455654/1‑73 (MQ=255)
                                              agaACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGc                 >  2:21734/1‑76 (MQ=255)
                                              agaACCGGTTGCCGTCTGGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGc                 >  2:1296216/1‑76 (MQ=255)
                                               gaACCGGTTGCCGTCTGGGGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAAACGCGATAAATGCg                >  2:1358454/1‑76 (MQ=255)
                                                               tGCGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGAt  >  1:26324/1‑76 (MQ=255)
                                                                 |                                                                         
CAACCCCTAACGCCAGCGCCAGCACCGACAGGGTAAACAGCTCGCGAGAACCGGTTGCCGCGCTGCGTGCCATAATCCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGAT  >  REL606/474591‑474729

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