Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 475,187 C→T G40D (GGT→GAT)  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,1870CT97.4% 140.9 / ‑3.4 39G40D (GGT→GAT) ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/17);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GCCCAGTTCAGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCG  >  REL606/475114‑475255
                                                                         |                                                                    
gCCCAGTTCAGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACAtcc                                                                    <  1:324308/76‑1 (MQ=255)
     gTTCAGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTaa                                                               <  1:1373949/76‑1 (MQ=255)
     gTTCAGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTaa                                                               <  1:318695/76‑1 (MQ=255)
      ttCAGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAAc                                                              >  1:1465822/1‑76 (MQ=255)
       tCAGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACa                                                             >  2:1285581/1‑76 (MQ=255)
        cAGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAg                                                            <  1:492233/76‑1 (MQ=255)
         aGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGa                                                           >  2:1268058/1‑76 (MQ=255)
           ccAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGata                                                         >  2:1351499/1‑76 (MQ=255)
              gTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATcc                                                      <  2:1510902/76‑1 (MQ=255)
               ttCCGGGGCAAGCTTGGTATCGGCAATAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATccc                                                     <  1:1413808/76‑1 (MQ=255)
               ttCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATccc                                                     >  2:1218189/1‑76 (MQ=255)
                  cGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCccacc                                                  <  2:1326095/76‑1 (MQ=255)
                  cGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCccacc                                                  <  1:529831/76‑1 (MQ=255)
                       cAAGCTGGGTAGCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAgagg                                             <  2:922337/76‑1 (MQ=255)
                            ttGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGaaa                                        >  1:953049/1‑76 (MQ=255)
                             tGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAAt                                       <  1:958598/76‑1 (MQ=255)
                               gTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATAc                                     <  1:1016291/76‑1 (MQ=255)
                                 aTCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGCGGAGAAATACGt                                   >  2:1158219/1‑76 (MQ=255)
                                 aTCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGt                                   <  1:740812/76‑1 (MQ=255)
                                   cGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAg                                 >  1:1223738/1‑76 (MQ=255)
                                   cGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAg                                 >  1:720466/1‑76 (MQ=255)
                                       aacaaAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTa                             <  2:1156887/76‑1 (MQ=255)
                                        acaaAGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTAt                            >  2:194426/1‑76 (MQ=255)
                                           aaGCCCGGAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTgg                         >  2:378822/1‑76 (MQ=255)
                                                 ggAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCa                   >  2:106952/1‑76 (MQ=255)
                                                 ggAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCa                   <  1:758113/76‑1 (MQ=255)
                                                 ggAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCa                   >  2:696953/1‑76 (MQ=255)
                                                  gAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCAt                  <  1:925651/76‑1 (MQ=255)
                                                  gAGTGAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCAt                  <  1:462932/76‑1 (MQ=255)
                                                      gAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCg              >  2:1003279/1‑76 (MQ=255)
                                                      gAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCg              >  2:991475/1‑76 (MQ=255)
                                                      gAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCg              >  2:926919/1‑76 (MQ=255)
                                                      gAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCg              >  2:893512/1‑76 (MQ=255)
                                                      gAATGGTCCTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCg              >  2:119935/1‑76 (MQ=255)
                                                             ccTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATaa       >  1:627793/1‑76 (MQ=255)
                                                              cTGCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATaaa      <  1:1351499/76‑1 (MQ=255)
                                                                gCCAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAgg    >  2:108868/1‑76 (MQ=255)
                                                                 ccAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGc   >  1:416708/1‑76 (MQ=255)
                                                                  cAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCg  <  1:52641/76‑1 (MQ=255)
                                                                  cAGCACATCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCg  <  1:893837/76‑1 (MQ=255)
                                                                         |                                                                    
GCCCAGTTCAGCCAGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCG  >  REL606/475114‑475255

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