Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 475,302 Δ1 bp coding (4/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
*REL606475,3020G.98.6% 344.4 / ‑2.3 70coding (4/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (46/23);  total (46/24)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAG  >  REL606/475232‑475371
                                                                      |                                                                     
gCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTcc                                                                 <  2:1255425/76‑1 (MQ=255)
gCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTcc                                                                 <  2:63243/76‑1 (MQ=255)
  cAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGt                                                               >  2:129579/1‑76 (MQ=255)
  cAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGt                                                               >  2:1909597/1‑76 (MQ=255)
   aGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGtc                                                              >  1:1169233/1‑76 (MQ=255)
    gCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGtct                                                             <  2:328777/76‑1 (MQ=255)
    gCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATCCCGtct                                                             >  2:697047/1‑76 (MQ=255)
      aTGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTcc                                                           >  2:1692306/1‑76 (MQ=255)
      aTGCAGAGAATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTcc                                                           <  1:190219/76‑1 (MQ=255)
          cGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCttt                                                        <  1:680223/76‑1 (MQ=255)
           gAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTc                                                      >  2:144309/1‑76 (MQ=255)
            aGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTcc                                                     <  1:592136/76‑1 (MQ=255)
            aGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTcc                                                     >  2:1534679/1‑76 (MQ=255)
            aGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTcc                                                     >  2:2190377/1‑76 (MQ=255)
               aTAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCtgg                                                  >  1:960149/1‑76 (MQ=255)
               aTAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCtgg                                                  >  1:701903/1‑76 (MQ=255)
               aTAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCtgg                                                  >  2:997368/1‑76 (MQ=255)
                tAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCtggt                                                 >  1:1349426/1‑76 (MQ=255)
                tAAAGGCGAGCACAAGGCCGACAACAATGGTGGTGATAAGCGGGGGGGCGTGATCATTCCGTCTCCTTTTCCtggt                                                 >  2:1419797/1‑76 (MQ=255)
                 aaaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCtggtg                                                >  2:1781450/1‑76 (MQ=255)
                 aaaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCtggtg                                                >  2:1303364/1‑76 (MQ=255)
                 aaaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCtggtg                                                >  2:782414/1‑76 (MQ=255)
                  aaGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCtggtgg                                               >  2:1415398/1‑76 (MQ=255)
                   aGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTggtggt                                              >  2:1749430/1‑76 (MQ=255)
                     gCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTa                                            <  1:1014445/76‑1 (MQ=255)
                       gAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTAtt                                          <  1:376074/76‑1 (MQ=255)
                       gAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTAtt                                          >  2:2236009/1‑76 (MQ=255)
                          cacaAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTc                                       <  1:1437704/76‑1 (MQ=255)
                          cacaAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTc                                       >  2:472274/1‑76 (MQ=255)
                           acaAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTcc                                      >  1:2131921/1‑76 (MQ=255)
                             aaGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAt                                    >  2:669047/1‑76 (MQ=255)
                             aaGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAt                                    >  2:1669038/1‑76 (MQ=255)
                              aGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAtt                                   <  2:1533030/76‑1 (MQ=255)
                              aGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCCTTTTCTGGTGGTTTTTTTCCAtt                                   >  1:1169193/1‑76 (MQ=255)
                                 ccgccAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAttttt                                >  1:261268/1‑76 (MQ=255)
                                 ccgccAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAttttt                                >  1:87771/1‑76 (MQ=255)
                                 ccgccAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAttttt                                <  2:1077917/76‑1 (MQ=255)
                                 ccgccAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAttttt                                >  1:1868754/1‑76 (MQ=255)
                                 ccgccAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCAttttt                                <  1:2160276/76‑1 (MQ=255)
                                  cgccAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTg                               >  2:1926211/1‑76 (MQ=255)
                                    ccAACAATGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGc                             <  1:78158/76‑1 (MQ=255)
                                     caacaaTGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGcc                            >  2:844288/1‑76 (MQ=255)
                                     caacaaTGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGcc                            >  1:1724298/1‑76 (MQ=255)
                                     caacaaTGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGcc                            >  1:1001397/1‑76 (MQ=255)
                                     caacaaTGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGcc                            >  1:514708/1‑76 (MQ=255)
                                     caacaaTGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGcc                            >  1:2188773/1‑76 (MQ=255)
                                       acaaTGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCgg                          >  1:618093/1‑76 (MQ=255)
                                        caaTGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCTATTTTTGGCCggg                         <  1:597092/76‑1 (MQ=255)
                                        caaTGGTGGTGATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCggg                         >  2:1293525/1‑76 (MQ=255)
                                              tggtgATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAAcc                   <  1:1909597/76‑1 (MQ=255)
                                               ggtgATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCa                  <  1:25172/76‑1 (MQ=255)
                                               ggtgATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCa                  >  1:1865064/1‑76 (MQ=255)
                                                gtgATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaa                 <  1:1308944/76‑1 (MQ=255)
                                                  gATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaaaa               >  1:1229807/1‑76 (MQ=255)
                                                  gATAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaaaa               >  1:1229828/1‑76 (MQ=255)
                                                   aTAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTTGCCGGGGAAACCAAAAt              >  1:1244327/1‑76 (MQ=255)
                                                   aTAAGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAAt              >  2:1486184/1‑76 (MQ=255)
                                                     aaGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTa            <  2:678289/76‑1 (MQ=255)
                                                     aaGCGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTa            >  2:1889463/1‑76 (MQ=255)
                                                        cGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAg         >  1:2191834/1‑76 (MQ=255)
                                                        cGGGGTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAg         >  1:1961605/1‑76 (MQ=255)
                                                         ggggTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAgg        <  2:297921/76‑1 (MQ=255)
                                                         ggggTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAgg        >  2:244344/1‑76 (MQ=255)
                                                         ggggTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAgg        >  2:1654016/1‑76 (MQ=255)
                                                          gggTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGt       >  1:6384/1‑76 (MQ=255)
                                                           ggTGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTa      <  1:1522062/76‑1 (MQ=255)
                                                             tGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAAt    <  1:941310/76‑1 (MQ=255)
                                                             tGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAAt    >  1:288548/1‑76 (MQ=255)
                                                             tGGCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGGAAt    >  1:676354/1‑76 (MQ=255)
                                                              ggCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATa   <  2:375690/76‑1 (MQ=255)
                                                              ggCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATa   <  1:887418/76‑1 (MQ=255)
                                                               gCGTGATCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAg  <  1:1271626/76‑1 (MQ=255)
                                                                      |                                                                     
GCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAG  >  REL606/475232‑475371

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: