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

CCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGT  >  REL606/474242‑474384
                                                                         |                                                                     
cccAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTgga                                                                     <  1:647683/76‑2 (MQ=255)
 ccAGTCCCGCCAGGATAGACGCGGACTCACCAATCTGCGCCAGGCTTGCGGCGATGGTTAATGCCGTAGGTTggag                                                                    >  1:71859/1‑74 (MQ=255)
 ccAGTCCCGCCAGGATAAACGCGAACTCACGAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTggag                                                                    >  1:1303712/1‑74 (MQ=255)
 ccAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTggag                                                                    >  1:1150785/1‑74 (MQ=255)
       ccGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCa                                                             <  1:385728/76‑1 (MQ=255)
              gATAAACGCGACCCCACCAACCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTc                                                      <  1:522704/76‑1 (MQ=255)
              gATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTc                                                      <  1:337533/76‑1 (MQ=255)
                tAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTTGAGTGACCAAACAGTCGc                                                    <  1:977348/76‑1 (MQ=255)
                tAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGc                                                    <  1:1081235/76‑1 (MQ=255)
                 aaaCGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCa                                                   >  2:468262/1‑76 (MQ=255)
                  aaCGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCAc                                                  <  2:1034499/76‑1 (MQ=255)
                      cgAACTCACCAATNTGCGCCAGGCTGGCGGCGATGGTTAATGCCGAACGTTGGAGTGACCAAACAGTCGCACCAgg                                              <  2:1298864/76‑1 (MQ=255)
                       gAACTCACCGATCTGCGCCAGGCTGGCGGCGATCGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGa                                             >  2:1228655/1‑76 (MQ=255)
                         aCTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGaaa                                           <  1:479366/76‑1 (MQ=255)
                          cTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGaaaa                                          <  1:774788/76‑1 (MQ=255)
                            cACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAAt                                        >  2:1156138/1‑76 (MQ=255)
                              ccAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGc                                      >  2:69950/1‑76 (MQ=255)
                                aaTCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCgg                                    >  1:737854/1‑76 (MQ=255)
                                   cTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTa                                 >  2:687284/1‑76 (MQ=255)
                                   cTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTa                                 <  2:396979/76‑1 (MQ=255)
                                     gcgcCAGGCTGGCGGCGATGGTTAATGCCGCACGTTGGAGAGACCCAACAGTCGCACCAGGANAAATGCGGCTAAc                               >  2:1308281/1‑76 (MQ=255)
                                      cgcCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACg                              <  2:975468/76‑1 (MQ=255)
                                              tGGCGGCGATGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCa                      <  1:1263207/76‑1 (MQ=255)
                                               ggcggcGATGGTTAATGTAGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAANGCGGCTAACGACTTACCaa                     <  1:1328930/76‑1 (MQ=255)
                                                   gcgATGGTTAATGCCGTACGTTGGAGTAACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAg                 <  1:572578/76‑1 (MQ=255)
                                                      aTGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGaa               >  2:443404/1‑76 (MQ=255)
                                                      aTGGTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGaat              >  2:693441/1‑76 (MQ=255)
                                                        ggTTAATGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGaataa            >  1:235763/1‑76 (MQ=255)
                                                             aTGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGcc       >  1:1200838/1‑76 (MQ=255)
                                                             aTGCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGcc       <  1:951054/76‑1 (MQ=255)
                                                               gCCGTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAg     <  2:85105/76‑1 (MQ=255)
                                                                 cGTACGTTGGAGTGACCAAACAGTCGCAACAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCCCTAGCg   >  1:1408842/1‑76 (MQ=255)
                                                                  gTACGTTGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGAGt  >  2:409477/1‑76 (MQ=255)
                                                                         |                                                                     
CCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGT  >  REL606/474242‑474384

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