Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 475,286 A→C L7R (CTT→CGT)  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
*REL606475,2860AC98.4% 254.5 / ‑3.4 62L7R (CTT→CGT) ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (45/16);  total (45/17)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAA  >  REL606/475213‑475358
                                                                         |                                                                        
agaAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGc                                                                        >  1:480477/1‑76 (MQ=255)
agaAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGGGGCGAGAAGc                                                                        >  2:262062/1‑76 (MQ=255)
 gaAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAcgcg                                                                       <  2:111172/76‑1 (MQ=255)
  aaaTACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCgg                                                                      >  1:1345464/1‑76 (MQ=255)
  aaaTACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCgg                                                                      >  2:859715/1‑76 (MQ=255)
  aaaTACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCgg                                                                      >  2:1090903/1‑76 (MQ=255)
  aaaTACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCgg                                                                      >  1:1962042/1‑76 (MQ=255)
   aaTACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCggg                                                                     >  1:1842372/1‑76 (MQ=255)
    aTACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCgggg                                                                    <  2:964705/76‑1 (MQ=255)
    aTACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCgggg                                                                    >  2:802486/1‑76 (MQ=255)
      aCGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTg                                                                  <  2:402680/76‑1 (MQ=255)
         tAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCg                                                               <  1:1994757/76‑1 (MQ=255)
          aGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCgg                                                              >  2:1319178/1‑75 (MQ=255)
           gTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTg                                                             >  1:1988047/1‑76 (MQ=255)
            tttATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGa                                                            <  1:364933/76‑1 (MQ=255)
            tttATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGa                                                            >  1:1953207/1‑76 (MQ=255)
               aTTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGc                                                         >  2:354091/1‑76 (MQ=255)
               aTTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGc                                                         >  1:949568/1‑76 (MQ=255)
                ttGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCa                                                        >  1:1938843/1‑76 (MQ=255)
                  ggCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCAtt                                                      <  1:458750/76‑1 (MQ=255)
                  ggCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCAtt                                                      >  2:1506394/1‑76 (MQ=255)
                  ggCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCAtt                                                      >  2:877138/1‑76 (MQ=255)
                   gCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTc                                                     >  2:581348/1‑76 (MQ=255)
                    ccAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTcc                                                    >  1:1810745/1‑76 (MQ=255)
                    ccAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTcc                                                    >  2:840746/1‑76 (MQ=255)
                    ccAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTcc                                                    >  1:1391274/1‑76 (MQ=255)
                             cGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCttt                                           <  2:1113612/76‑1 (MQ=255)
                              gAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCtttt                                          >  2:1894683/1‑76 (MQ=255)
                              gAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCtttt                                          <  1:2015494/76‑1 (MQ=255)
                              gAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGAGGCGTGATGCATTCTGGCTCCtttt                                          >  1:481201/1‑76 (MQ=255)
                                 gATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCt                                       <  1:405836/76‑1 (MQ=255)
                                  aTAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCtg                                      >  2:1888750/1‑76 (MQ=255)
                                  aTAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCtg                                      >  1:1754625/1‑76 (MQ=255)
                                  aTAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCtg                                      <  2:1480007/76‑1 (MQ=255)
                                    aanGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCtggt                                    >  2:1425096/1‑76 (MQ=255)
                                     acngCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCtggtg                                   >  2:1423407/4‑76 (MQ=255)
                                     aaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCtggtg                                   >  2:135965/1‑76 (MQ=255)
                                      aGGCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCtggtgg                                  >  2:243427/1‑76 (MQ=255)
                                       ggCGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTggtggt                                 >  2:795578/1‑76 (MQ=255)
                                         cGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTa                               <  2:293067/76‑1 (MQ=255)
                                         cGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTa                               >  2:1724704/1‑76 (MQ=255)
                                         cGAGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTa                               >  2:1012835/1‑76 (MQ=255)
                                           aGCACAAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGttnct                             <  2:1426681/76‑4 (MQ=255)
                                             cacaAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGt                           >  2:492497/1‑76 (MQ=255)
                                             cacaAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGt                           >  2:1897892/1‑76 (MQ=255)
                                               caAGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTcc                         >  2:1477501/1‑76 (MQ=255)
                                                aaGGCCGCCAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCa                        >  2:288268/1‑76 (MQ=255)
                                                   gccgccAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAttt                     >  1:1549912/1‑76 (MQ=255)
                                                   gccgccAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAttt                     >  1:241666/1‑76 (MQ=255)
                                                       ccAACAATGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTgg                 >  1:718776/1‑76 (MQ=255)
                                                        caacaaTGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGc                >  1:1517140/1‑76 (MQ=255)
                                                         aacaaTGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGcc               <  2:457318/76‑1 (MQ=255)
                                                          acaaTGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCg              >  2:1765159/1‑76 (MQ=255)
                                                            aaTGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCggg            >  1:889365/1‑76 (MQ=255)
                                                            aaTGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCggg            >  1:87700/1‑76 (MQ=255)
                                                            aaTGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCggg            >  1:852073/1‑76 (MQ=255)
                                                             aTGGTGGTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGa           <  2:651774/76‑1 (MQ=255)
                                                              tggtggTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaa          >  2:167558/1‑76 (MQ=255)
                                                              tggtggTGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaa          <  1:435957/76‑1 (MQ=255)
                                                               ggtggtGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaa         <  1:1125506/76‑1 (MQ=255)
                                                               ggtggtGATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaa         <  2:1217199/76‑1 (MQ=255)
                                                                gtggtgATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaaa        >  1:82246/1‑76 (MQ=255)
                                                                  ggtgATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAAcc      <  2:610911/76‑1 (MQ=255)
                                                                    tgATACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaa    >  2:471101/1‑76 (MQ=255)
                                                                      aTACGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaaaa  >  2:543422/1‑76 (MQ=255)
                                                                         |                                                                        
AGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAA  >  REL606/475213‑475358

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