Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 505,179 C→T 68.3% L348L (CTG→TTG)  ushA → bifunctional UDP‑sugar hydrolase and 5'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb505,1790CT68.3% 48.5 / 34.2 41L348L (CTG→TTG) ushAbifunctional UDP‑sugar hydrolase and 5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (2/11);  new base T (12/16);  total (14/27)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAA  >  W3110S.gb/505125‑505225
                                                      |                                              
agcgagcgCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTc                                  >  1:154474/1‑69 (MQ=255)
agcgagcgCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTc                                  >  1:1162300/1‑69 (MQ=255)
agcgagcgCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTc                                  >  1:1795462/1‑69 (MQ=255)
agcgagcgCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTc                                  >  1:1390726/1‑69 (MQ=255)
agcgagcgCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTc                                  >  1:1243705/1‑69 (MQ=255)
agcgagcgCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTc                                  >  1:190013/1‑69 (MQ=255)
     gcgcGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACa                           <  1:991161/71‑1 (MQ=255)
     gcgcGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACa                           <  1:627576/71‑1 (MQ=255)
      cgcgTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACa                           <  1:1391169/70‑1 (MQ=255)
                           aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt     >  1:1302934/1‑71 (MQ=255)
                           aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt     >  1:1349099/1‑71 (MQ=255)
                           aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt     >  1:1871699/1‑71 (MQ=255)
                           aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt     >  1:1815333/1‑71 (MQ=255)
                           aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt     >  1:651091/1‑71 (MQ=255)
                           aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt     >  1:1663909/1‑71 (MQ=255)
                           aTCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt     >  1:1732515/1‑71 (MQ=255)
                           aTCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt     >  1:1672894/1‑71 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCATAGCGCAGCTGGAAGTGaa  <  1:550376/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:220491/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:932623/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:861552/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:123362/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:563650/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:407435/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:366099/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1554187/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1396281/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1445903/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1468656/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1028561/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1491074/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1848987/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:528330/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1483405/71‑1 (MQ=255)
                              gCTGAAAACCAGCAAACGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1652348/71‑1 (MQ=255)
                               cTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1329325/70‑1 (MQ=255)
                               cTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1173189/70‑1 (MQ=255)
                                                aTCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1127570/53‑1 (MQ=255)
                                                aTCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:937215/53‑1 (MQ=255)
                                                aTCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:347179/53‑1 (MQ=255)
                                                aTCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa  <  1:1954223/53‑1 (MQ=255)
                                                      |                                              
AGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAA  >  W3110S.gb/505125‑505225

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