Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 505,179 C→T 100% 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,1790CT77.8% 49.3 / 11.7 27L348L (CTG→TTG) ushAbifunctional UDP‑sugar hydrolase and 5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base T (8/13);  total (8/19)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATGG  >  W3110S.gb/505113‑505243
                                                                  |                                                                
gAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGt                                                               <  1:2204998/70‑1 (MQ=255)
     cGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCACCAAATGATCTCGTTGTTATCAc                                                         <  1:1767810/71‑1 (MQ=255)
                    cgTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCAc                                                         <  1:1538559/56‑1 (MQ=255)
                      tGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGc                                        <  1:2026867/71‑1 (MQ=255)
                        cTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCa                                       >  1:2603192/1‑70 (MQ=255)
                        cTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCa                                       >  1:2727539/1‑70 (MQ=255)
                           tACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAg                                    <  1:263454/70‑1 (MQ=255)
                               cTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGaa                                               <  1:1132095/55‑1 (MQ=255)
                                tCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGc                              >  1:552548/1‑71 (MQ=255)
                                tCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaa                                            <  1:508422/57‑1 (MQ=255)
                                  cTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTg                            >  1:181399/1‑71 (MQ=255)
                                          gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                    <  1:1574540/71‑1 (MQ=255)
                                          gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                    <  1:2275332/71‑1 (MQ=255)
                                          gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                    <  1:2303267/71‑1 (MQ=255)
                                          gCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                    <  1:2676280/71‑1 (MQ=255)
                                          gCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                    <  1:213314/71‑1 (MQ=255)
                                           cTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaaa                   <  1:230227/71‑1 (MQ=255)
                                              aaaaCCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATa                <  1:1598660/71‑1 (MQ=255)
                                              aaaaCCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaaaa                  <  1:1446504/69‑1 (MQ=255)
                                                 aCCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGc             >  1:929544/1‑71 (MQ=255)
                                                     gCAAATGATCTCGCTGTTATCACCGTTCCAGAACaaa                                           <  1:2382537/37‑1 (MQ=255)
                                                       aaaTGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAAcc       <  1:1553395/71‑1 (MQ=255)
                                                         aTGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                    <  1:2682698/56‑1 (MQ=255)
                                                         aTGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                    <  1:318951/56‑1 (MQ=255)
                                                          tGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaaaa                  >  1:2635358/1‑57 (MQ=255)
                                                            aTCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATgg  >  1:2435661/1‑71 (MQ=255)
                                                            aTCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATgg  >  1:1922110/1‑71 (MQ=255)
                                                                  |                                                                
GAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATGG  >  W3110S.gb/505113‑505243

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: