Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 517,941 A→G 23.1% S145P (TCA→CCA)  ybbO ← predicted oxidoreductase with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb517,9410AG23.1% 45.4 / 11.7 26S145P (TCA→CCA) ybbOpredicted oxidoreductase with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (12/8);  new base G (4/2);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGT  >  W3110S.gb/517871‑518006
                                                                      |                                                                 
gCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCgatga                                                                   <  1:3426357/71‑1 (MQ=255)
         gcATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGGTGtcatca                                                           >  1:2575702/1‑70 (MQ=255)
         gcATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATcac                                                          <  1:2356448/71‑1 (MQ=255)
          cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGtc                                                               <  1:2976653/65‑1 (MQ=255)
            tatTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGt                                                                <  1:2937683/62‑1 (MQ=255)
                     gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAAt                                                       >  1:248719/1‑62 (MQ=255)
                     gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAAt                                                       >  1:3663530/1‑62 (MQ=255)
                         cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGt                                                                <  1:417806/49‑1 (MQ=255)
                             aGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGatgatg                                                                 >  1:243907/1‑44 (MQ=255)
                             aGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGGTGTCATCACAATACGCCCTTCACCGTGCg                                      >  1:1018358/1‑71 (MQ=255)
                             aGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATa                                                      >  1:122681/1‑55 (MQ=255)
                               cgcCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATcaca                                                         <  1:335514/50‑1 (MQ=255)
                                gcCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGt                                    >  1:3413030/1‑70 (MQ=255)
                                   aCGACCCGGCGTGGAGATTAATCCCATCACCgatgat                                                                  >  1:1750440/1‑37 (MQ=255)
                                   aCGACCCGGCGTGGAGATTAATCCCATCACCgatgat                                                                  >  1:949544/1‑37 (MQ=255)
                                   aCGACCCGGCGTGGAGATTAATCCCATCACCgatgat                                                                  >  1:2946804/1‑37 (MQ=255)
                                           gCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGt                        <  1:1642686/71‑1 (MQ=255)
                                                gagaTTAATCCCATCACCGATGGTGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAg                   >  1:2896437/1‑71 (MQ=255)
                                                     tAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCa              >  1:2356688/1‑71 (MQ=255)
                                                      aaTCCCATCACCGATGGTGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCa              <  1:1668565/70‑1 (MQ=255)
                                                              cACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCggg                         >  1:2837362/1‑51 (MQ=255)
                                                              cACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCggg                         >  1:1314809/1‑51 (MQ=255)
                                                               aCCGATGGTGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg    >  1:710081/1‑71 (MQ=255)
                                                               aCCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg    >  1:2184134/1‑71 (MQ=255)
                                                               aCCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg    >  1:2069286/1‑71 (MQ=255)
                                                               aCCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg    >  1:3615067/1‑71 (MQ=255)
                                                               aCCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg    >  1:1781127/1‑71 (MQ=255)
                                                               aCCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg    >  1:977756/1‑71 (MQ=255)
                                                                ccGATGGTGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTgg   <  1:1965999/71‑1 (MQ=255)
                                                                 cGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGgt  <  1:3851160/71‑1 (MQ=255)
                                                                 cGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGgt  <  1:612293/71‑1 (MQ=255)
                                                                      |                                                                 
GCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGT  >  W3110S.gb/517871‑518006

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