Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 474,186 T→A I374F (ATC→TTC)  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,1860TA97.6% 131.3 / NA 42I374F (ATC→TTC) ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (17/23);  minor base G (0/1);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGG  >  REL606/474113‑474256
                                                                         |                                                                      
tCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACg                                                                      <  1:1389946/76‑1 (MQ=255)
     cGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCt                                                                 <  2:766637/76‑1 (MQ=255)
      gTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTg                                                                <  2:634895/76‑1 (MQ=255)
        ttCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGcc                                                              <  1:274989/76‑1 (MQ=255)
           ggTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGt                                                           <  2:1141142/76‑1 (MQ=255)
           ggTCTCCGCCAGATATTTCCCCAGTAGTGCAAACAGTACGGGGTTGAGCATAATCGACAGGACCGCCCCTGCCAGt                                                           <  2:1395704/76‑1 (MQ=255)
             tCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTAc                                                         <  1:1715727/76‑1 (MQ=255)
               ttCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCa                                                       >  1:1138496/1‑76 (MQ=255)
                    cAGATATTCCGCAAGTAGTACGACAGGTACAGGGTTGAGCATAACGAAAAGGAACGCCCCTGCCAGTACCAGGttt                                                  <  2:557697/76‑1 (MQ=255)
                      gATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTg                                                <  2:244930/76‑1 (MQ=255)
                          tttCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCg                                            <  1:56671/76‑1 (MQ=255)
                           ttCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCgg                                           <  2:1351759/76‑1 (MQ=255)
                            tctcCAGTGGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGc                                          <  1:1007671/76‑1 (MQ=255)
                             ctcCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGcc                                         <  1:1306100/76‑1 (MQ=255)
                               ccAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTg                                       >  1:1533838/1‑76 (MQ=255)
                               ccAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTg                                       <  2:1519068/76‑1 (MQ=255)
                                 agtagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCg                                     >  1:1191635/1‑76 (MQ=255)
                                 agtagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCg                                     >  1:424091/1‑76 (MQ=255)
                                 agtagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCg                                     >  1:572860/1‑76 (MQ=255)
                                  gtagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCgg                                    >  1:166908/1‑76 (MQ=255)
                                  gtagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCgg                                    <  2:777159/76‑1 (MQ=255)
                                   tggtGCGAACAGTCCCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGc                                   <  1:456364/74‑1 (MQ=255)
                                   tagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGc                                   >  2:462269/1‑76 (MQ=255)
                                   tagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGc                                   >  1:1531418/1‑76 (MQ=255)
                                   taggGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGc                                   <  2:166302/76‑1 (MQ=255)
                                      tGCGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAga                                >  1:229723/1‑75 (MQ=255)
                                       gCGAACAGTCCCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTa                               <  2:650963/76‑1 (MQ=255)
                                        cGAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCAAGTACCAGGTTTTGTCCGCCCTGCGCCAGTaa                              >  2:189319/1‑76 (MQ=255)
                                         gAACAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTaaa                             <  1:1487325/76‑1 (MQ=255)
                                           aCAGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAAtt                           <  2:296373/76‑1 (MQ=255)
                                             aGTACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGGAAATTCa                         >  2:1229161/1‑76 (MQ=255)
                                               tgccGGGTTGAGCATAATCGACAGGAGCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAAt                       <  1:556104/74‑1 (MQ=255)
                                               tACCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAAt                       <  1:420079/76‑1 (MQ=255)
                                                cccGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATg                      <  2:802250/75‑1 (MQ=255)
                                                aCCGGGTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTATCCGGCCTGCGGCAGTAAATTCAATg                      >  1:796998/1‑76 (MQ=255)
                                                 ccGGGTTGAGCATAATCGACAGAAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGc                     <  2:1659392/76‑1 (MQ=255)
                                                     gTTGAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCAtt                 >  2:1203902/1‑76 (MQ=255)
                                                        gAGCATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTccc              <  2:1283831/76‑1 (MQ=255)
                                                           cATAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGt           >  1:861686/1‑76 (MQ=255)
                                                            aTAATCGACGGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTc          >  1:1336885/1‑76 (MQ=255)
                                                            aTAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTc          <  2:374538/76‑1 (MQ=255)
                                                             tAATCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTcc         >  2:1483529/1‑76 (MQ=255)
                                                                tCGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGc      <  1:1127264/76‑1 (MQ=255)
                                                                 cGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGcc     <  2:229723/76‑1 (MQ=255)
                                                                 cGACAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGcc     <  1:364723/76‑1 (MQ=255)
                                                                    cAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTGCCGCCAgg  >  2:1505561/1‑76 (MQ=255)
                                                                    cAGGAACGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAgg  <  1:1299651/76‑1 (MQ=255)
                                                                         |                                                                      
TCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGG  >  REL606/474113‑474256

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