Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 474,360 Δ1 bp coding (946/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,3590T.97.6% 190.5 / ‑1.5 42coding (947/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (17/24);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

GCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAAC  >  REL606/474286‑474430
                                                                         |                                                                       
gCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACtac                                                                       <  2:1312782/76‑3 (MQ=255)
         gATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGaa                                                             <  2:1366555/76‑1 (MQ=255)
          aTGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGaat                                                            >  2:785604/1‑76 (MQ=255)
          aTGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGaat                                                            <  1:351702/76‑1 (MQ=255)
             gTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGaataat                                                         >  1:888498/1‑76 (MQ=255)
             gTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGaataat                                                         <  1:302552/76‑1 (MQ=255)
             gTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGaataat                                                         >  1:701126/1‑76 (MQ=255)
             gTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGaataat                                                         >  1:6712/1‑76 (MQ=255)
               tAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCg                                                       >  1:1280555/1‑76 (MQ=255)
               tAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCg                                                       <  1:1314161/76‑1 (MQ=255)
               tAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCg                                                       >  2:596517/1‑76 (MQ=255)
                  tGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCa                                                    <  2:487285/76‑1 (MQ=255)
                  tGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCa                                                    <  2:378713/76‑1 (MQ=255)
                     cGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCg                                                 >  1:914909/1‑76 (MQ=255)
                        aCGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCg                                              >  1:734082/1‑76 (MQ=255)
                                aGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACt                                      <  2:892013/76‑1 (MQ=255)
                                aGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACt                                      <  2:751573/76‑1 (MQ=255)
                                   ggccAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGcc                                   <  2:36316/74‑1 (MQ=255)
                                   gACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGcc                                   <  2:1516351/76‑1 (MQ=255)
                                   gACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGcc                                   >  2:1478226/1‑76 (MQ=255)
                                      cAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGc                                <  2:431901/76‑1 (MQ=255)
                                        aaCAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCgg                              >  2:1522944/1‑76 (MQ=255)
                                         aCAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGt                             >  1:1354613/1‑76 (MQ=255)
                                          cAGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGtt                            >  1:1150206/1‑76 (MQ=255)
                                           aGTCGCACCAGTAAAAATGCGTCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTg                           <  2:591816/76‑1 (MQ=255)
                                           aGTCGCACCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTg                           <  2:1299199/76‑1 (MQ=255)
                                                cACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGa                       >  1:1314501/1‑76 (MQ=255)
                                                 aCCAGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAAt                     <  1:1579420/76‑1 (MQ=255)
                                                    agaaaaaaTGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAg                  <  1:657160/73‑1 (MQ=255)
                                                    aGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAg                  <  2:1365936/76‑1 (MQ=255)
                                                    aGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAg                  >  2:128006/1‑76 (MQ=255)
                                                    aGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAg                  <  1:1421114/76‑1 (MQ=255)
                                                    aGGAAAAATGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAg                  >  1:478845/1‑76 (MQ=255)
                                                       aaaaaTGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAAt               <  2:844706/76‑1 (MQ=255)
                                                       aaaaaTGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAAt               <  1:1507537/76‑1 (MQ=255)
                                                            tGCGGCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACg          >  2:439607/1‑76 (MQ=255)
                                                              cgtcTAACGGCTACCAAACAGAATAAGCGCCAGCGTCGCCAGCGCTGCCAGCGGTTGCTGAATCAGAATTAACGGa        <  1:1368163/73‑1 (MQ=255)
                                                               ggCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGAt       <  1:1157292/76‑1 (MQ=255)
                                                               ggCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGAt       >  2:226032/1‑76 (MQ=255)
                                                               ggCTAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGAt       <  1:641601/76‑1 (MQ=255)
                                                                 cTAACGGCTACCAAACAGAAAAATCGCCAGCGGCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCa     <  2:1232042/76‑1 (MQ=255)
                                                                  tAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCaa    <  1:921268/76‑1 (MQ=255)
                                                                  tAACGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCaa    >  1:1459622/1‑76 (MQ=255)
                                                                  tAACGACTACCAAACAGAAAAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCaa    <  2:770382/76‑1 (MQ=255)
                                                                    aCGACTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAaac  <  1:724010/76‑1 (MQ=255)
                                                                         |                                                                       
GCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAAC  >  REL606/474286‑474430

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