Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 474,210 Δ1 bp coding (1096/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,2100G.96.4% 121.5 / ‑0.0 28coding (1096/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (13/14);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

CCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGC  >  REL606/474144‑474282
                                                                  |                                                                        
ccAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGc                                                                >  1:209876/1‑76 (MQ=255)
  agtagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGGTTTTCCGGCCCGCgg                                                              >  2:1402843/1‑76 (MQ=255)
  agtaggGCGAACAGTACGGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCgg                                                              <  1:468713/76‑1 (MQ=255)
     agtGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAg                                                           <  1:1121518/76‑1 (MQ=255)
     agtGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAg                                                           >  2:1715758/1‑76 (MQ=255)
      gtGCGAACAGTACCGGGTGGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGt                                                          <  1:487832/76‑1 (MQ=255)
        gCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTa                                                         <  1:1588874/76‑1 (MQ=255)
         cGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTaaa                                                       <  1:1737909/76‑1 (MQ=255)
          gAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAAt                                                      >  2:221564/1‑76 (MQ=255)
            aCAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTc                                                    >  1:332062/1‑76 (MQ=255)
             cAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCa                                                   <  1:187889/76‑1 (MQ=255)
              aGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCaa                                                  >  1:1788708/1‑76 (MQ=255)
                 aCCGGGTTGAGCATAATCGACAGGATCGCCCCCGCCAGTACCAGGCTTTTCCGGCCTCCCGCAGGAAACACAAtac                                               >  2:1479318/1‑74 (MQ=255)
                  cggggTTGAGCATACTCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGcc                                              <  1:717981/74‑1 (MQ=255)
                  ccGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGcc                                              <  1:1499670/76‑1 (MQ=255)
                   cGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCa                                             >  2:493956/1‑76 (MQ=255)
                   cGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCa                                             >  2:847788/1‑76 (MQ=255)
                         gAGCATAGGCGTCAGGATCCCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCa                                       <  1:113133/76‑1 (MQ=255)
                            cATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTc                                    >  1:1295215/1‑76 (MQ=255)
                              tAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTccc                                  <  1:15345/76‑1 (MQ=255)
                                 ccGACAGGAGCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGcc                               <  1:739271/75‑1 (MQ=255)
                                    aCAGGATCCCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAgg                            <  1:941778/76‑1 (MQ=255)
                                      aGGATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGAt                          >  2:23280/1‑76 (MQ=255)
                                       ggATCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATa                         <  1:1754352/76‑1 (MQ=255)
                                         aTCGCCCCTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATaaa                       <  1:987067/76‑1 (MQ=255)
                                             ccccTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAAcgcg                   <  1:1193102/76‑1 (MQ=255)
                                                cTGCCAGTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAAc                >  1:1393954/1‑76 (MQ=255)
                                                      gTACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCa          >  2:52670/1‑76 (MQ=255)
                                                       tACCAGGTTTTTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCaa         >  1:19701/1‑76 (MQ=255)
                                                               ttttTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTgcgc  >  2:599217/1‑76 (MQ=255)
                                                                  |                                                                        
CCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGC  >  REL606/474144‑474282

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