Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 505,179 C→T 34.6% 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,1790CT34.6% 59.7 / 47.4 52L348L (CTG→TTG) ushAbifunctional UDP‑sugar hydrolase and 5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (22/12);  new base T (12/6);  total (34/18)
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

GGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATGGT  >  W3110S.gb/505111‑505244
                                                                    |                                                                 
gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                                                                   <  1:1543444/69‑1 (MQ=255)
gggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                                                                   >  1:3040989/1‑69 (MQ=255)
 ggAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCt                                                                  >  1:1301540/1‑69 (MQ=255)
    aGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                                                               <  1:1925370/69‑1 (MQ=255)
    aGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                                                               <  1:1822575/69‑1 (MQ=255)
     gACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGc                                                                   >  1:496256/1‑64 (MQ=255)
          gAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCAcc                                                         >  1:2829540/1‑69 (MQ=255)
                 gagcgCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                                                               >  1:784765/1‑56 (MQ=255)
                 gagcgCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                                                               <  1:162327/56‑1 (MQ=255)
                            ttACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATGAc                                                          <  1:305918/50‑1 (MQ=255)
                             tACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCAc                                                          <  1:2814126/49‑1 (MQ=255)
                               cacTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAg                                     >  1:200506/1‑68 (MQ=255)
                               cacTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAg                                     >  1:3408348/1‑68 (MQ=255)
                               cacTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAg                                     >  1:2121672/1‑68 (MQ=255)
                                 cTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAgcgc                                  >  1:1584443/1‑69 (MQ=255)
                                 cTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAgcgc                                  >  1:1202953/1‑69 (MQ=255)
                                 cTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAgcgc                                  >  1:2528656/1‑69 (MQ=255)
                                  tCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCaa                                       <  1:1832932/63‑1 (MQ=255)
                                  tCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCa                                 >  1:61088/1‑69 (MQ=255)
                                  tCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCa                                 >  1:125444/1‑69 (MQ=255)
                                   ccTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGtt                                                               <  1:3257120/38‑1 (MQ=255)
                                   ccTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACTGTTCCAGAACAAAGGCAAAGCGCAg                                >  1:1421642/1‑69 (MQ=255)
                                   ccTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAg                                >  1:832666/1‑69 (MQ=255)
                                    cTGAAATCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGc                               >  1:3087641/1‑69 (MQ=255)
                                       aaaTCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAg                                <  1:1459016/65‑1 (MQ=255)
                                       aaaTCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCt                              >  1:2840140/1‑67 (MQ=255)
                                         aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACa                                              <  1:2349735/49‑1 (MQ=255)
                                         aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACa                                              >  1:570732/1‑49 (MQ=255)
                                         aTCGCTGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACa                                              >  1:1982671/1‑49 (MQ=255)
                                         aTCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCa                                        >  1:1561716/1‑55 (MQ=255)
                                          tCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTgg                            >  1:2324579/1‑66 (MQ=255)
                                           cGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGt                        <  1:1217615/69‑1 (MQ=255)
                                            gCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGaa                                                >  1:3334260/1‑44 (MQ=255)
                                            gCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACaa                                             >  1:2413751/1‑47 (MQ=255)
                                              tGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                     <  1:1684807/69‑1 (MQ=255)
                                              tGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                     <  1:3368135/69‑1 (MQ=255)
                                              tGAAAACCAGCAAATGATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                     <  1:2060133/69‑1 (MQ=255)
                                              tGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                     <  1:3665319/69‑1 (MQ=255)
                                              tGAAAACCAGCAAATGATCTCGCTGTTATCACCATTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGaa                     <  1:2631602/69‑1 (MQ=255)
                                                    ccAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAg                >  1:1519447/1‑68 (MQ=255)
                                                            tGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCa       <  1:1646748/69‑1 (MQ=255)
                                                            tGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCa       >  1:32182/1‑69 (MQ=255)
                                                            tGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCa       >  1:2537520/1‑69 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTgg                            <  1:1574790/47‑1 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTgg                            <  1:1280128/47‑1 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATa                 >  1:2098780/1‑58 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATa                 >  1:525285/1‑58 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATa                 >  1:3478155/1‑58 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATa                 >  1:2719850/1‑58 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATa                 >  1:2057623/1‑58 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATa                 >  1:249619/1‑58 (MQ=255)
                                                             gATCTCGTTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATa                 >  1:20607/1‑58 (MQ=255)
                                                              aTCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCaa      >  1:370280/1‑68 (MQ=255)
                                                                 tcGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATGgt  >  1:551748/1‑69 (MQ=255)
                                                                 tcGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATGgt  >  1:1192822/1‑69 (MQ=255)
                                                                    |                                                                 
GGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATGGT  >  W3110S.gb/505111‑505244

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