Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 473,924 A→C V461G (GTC→GGC)  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
*REL606473,9240AC96.7% 224.4 / ‑1.2 60V461G (GTC→GGC) ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (26/32);  total (27/33)
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.31e-01

GCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAGGCGA  >  REL606/473861‑473996
                                                                |                                                                        
gCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCtctc                                                              <  2:81371/76‑1 (MQ=255)
 cacaTTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCTCGCTCTCg                                                             >  1:400967/1‑76 (MQ=255)
  acaTTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGc                                                            <  2:990588/76‑1 (MQ=255)
  acaTTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGc                                                            <  2:812877/76‑1 (MQ=255)
  acaTTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGc                                                            >  1:348389/1‑76 (MQ=255)
  acaTTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGc                                                            <  2:1129595/76‑1 (MQ=255)
  acaGTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGc                                                            <  2:557000/76‑1 (MQ=255)
   ccatCCAGGAGCCCCAGTTGCCATAATTTCTCCGTCCCCCGCATCCCCCAATCCGGCGCGCCCCCCGCGCTCTCGc                                                            <  2:288951/73‑1 (MQ=255)
   caTTCCAGATGCGCCAGTTGCATAATTTCTTCGTGCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCa                                                           <  2:743352/76‑1 (MQ=255)
    aTTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGcgg                                                          >  1:1335543/1‑74 (MQ=255)
      tCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCt                                                        <  2:517057/76‑1 (MQ=255)
      tCCAGATGCGCCAGTTGCATAATTACTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCt                                                        <  2:378425/76‑1 (MQ=255)
       ccAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCtc                                                       >  1:1249266/1‑76 (MQ=255)
       ccAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCtc                                                       <  2:811792/76‑1 (MQ=255)
         aGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCtcat                                                     <  1:1338504/76‑1 (MQ=255)
         aGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCtcat                                                     <  1:383839/76‑1 (MQ=255)
         aGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCtcat                                                     >  1:1054773/1‑76 (MQ=255)
          gATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCtcatc                                                    <  2:1502424/76‑1 (MQ=255)
            tGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCaa                                                  >  1:643845/1‑76 (MQ=255)
             gcgcCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAac                                                 <  1:956408/76‑1 (MQ=255)
             gcgcCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAac                                                 <  2:491642/76‑1 (MQ=255)
             gcgcCAGTTGCATAATGTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAac                                                 <  2:1631278/76‑1 (MQ=255)
              cgcCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAaca                                                <  1:1558795/76‑1 (MQ=255)
                  aGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACggg                                            <  1:121996/76‑1 (MQ=255)
                    ttGCAAAATTCCTGCGTTCGCCGCACTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGtt                                          <  2:124842/76‑1 (MQ=255)
                     tGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTc                                         <  2:1170213/76‑1 (MQ=255)
                     tGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTc                                         >  1:66214/1‑76 (MQ=255)
                     tGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGAGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTc                                         >  2:54236/1‑76 (MQ=255)
                        cATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGt                                       >  2:1204095/1‑76 (MQ=255)
                        cATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGt                                       >  2:794124/1‑76 (MQ=255)
                         aTAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTg                                      <  1:1486561/76‑1 (MQ=255)
                         aTAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTg                                      <  1:1816340/76‑1 (MQ=255)
                         aTAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTg                                      <  1:1330007/76‑1 (MQ=255)
                            aTTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACg                                   <  1:1873758/76‑1 (MQ=255)
                               tcttcGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGtct                                <  2:55230/76‑1 (MQ=255)
                               tcttcGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGtct                                <  2:723851/76‑1 (MQ=255)
                                cttcGTTCGCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGtctc                               <  2:518945/76‑1 (MQ=255)
                                  tcgttcgCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCaa                             >  2:65975/1‑76 (MQ=255)
                                  tcgttcgCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCaa                             >  1:37502/1‑76 (MQ=255)
                                   cgttcgCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATAAACACGGGTTCGTGACGTCTCAAt                            >  1:1002391/1‑76 (MQ=255)
                                    gttcgCCGCATTGCCCAATACTGCGCGGCCCCCGTGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATc                           <  2:773250/76‑1 (MQ=255)
                                      tcgCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCGATcac                         <  2:1351511/76‑1 (MQ=255)
                                      tcgCCGCATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATcac                         >  2:1646796/1‑76 (MQ=255)
                                       cgccgcATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATcacc                        >  2:600977/1‑76 (MQ=255)
                                       cgccgcATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATcacc                        >  2:1729609/1‑76 (MQ=255)
                                       cgccgcATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATcacc                        >  1:678202/1‑76 (MQ=255)
                                       cgccgcATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATcacc                        >  2:951192/1‑76 (MQ=255)
                                        gccgcATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATcacca                       >  1:799374/1‑76 (MQ=255)
                                        gccgcATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGCTCGTGACGTCTCAATcacca                       >  1:1023499/1‑76 (MQ=255)
                                           gcATTGCCCAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGTCGTCTCAATCaccaccg                    >  2:518285/1‑75 (MQ=255)
                                                  ccAATTCTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAat             <  2:255949/76‑1 (MQ=255)
                                                  ccAATACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAat             >  1:1060050/1‑76 (MQ=255)
                                                    aaTACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCCTGACGTCTCAATCACCACCAGCGGAatat           <  2:946176/76‑1 (MQ=255)
                                                     aTACTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATc          >  2:1377828/1‑76 (MQ=255)
                                                       aCTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCag        >  1:401644/1‑76 (MQ=255)
                                                       aCTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCag        <  1:1148991/76‑1 (MQ=255)
                                                       aCTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCag        <  2:939216/76‑1 (MQ=255)
                                                       aCTGCGCGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCag        <  2:1010301/76‑1 (MQ=255)
                                                          gcgcGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAgg     <  1:833607/76‑1 (MQ=255)
                                                          gcgcGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAgg     >  2:1683683/1‑76 (MQ=255)
                                                          gcgcGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAgg     <  1:1096197/76‑1 (MQ=255)
                                                          gcgcGGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAgg     >  2:971439/1‑76 (MQ=255)
                                                             cgGCCCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAGGCGa  >  2:112516/1‑76 (MQ=255)
                                                                |                                                                        
GCACATTCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGCCGCATTGCCCAATACTGCGCGGACCCCGCGCTCTCGCAGCTCATCAACACGGGTTCGTGACGTCTCAATCACCACCAGCGGAATATCAGAGGCGA  >  REL606/473861‑473996

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: