Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 429,850 T→C 12.8% Y8H (TAT→CAT)  yajD → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb429,8500TC12.8% 116.6 / 10.6 47Y8H (TAT→CAT) yajDconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base C (3/3);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

ATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGT  >  W3110S.gb/429781‑429918
                                                                     |                                                                    
aTCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTa                                                                     <  1:825594/71‑1 (MQ=255)
   ttGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAGAAAACCATgc                                                                  >  1:3606780/1‑71 (MQ=255)
   ttGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATgc                                                                  >  1:3591787/1‑71 (MQ=255)
          aCTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTa                                                           <  1:2226495/71‑1 (MQ=255)
            tGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACt                                                                      >  1:3154111/1‑58 (MQ=255)
             gcgcGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGa                                                         >  1:971546/1‑70 (MQ=255)
              cgcgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGaaa                                                       <  1:3673827/71‑1 (MQ=255)
                cgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGc                                                     <  1:1222542/71‑1 (MQ=255)
                        tACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTAt                                                                    >  1:2012115/1‑48 (MQ=255)
                            tttCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTCTgc                                                                  <  1:3755596/46‑1 (MQ=255)
                            tttCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaa                                         <  1:3464396/71‑1 (MQ=255)
                            tttCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGa                                          <  1:3406462/70‑1 (MQ=255)
                             ttCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaa                                         >  1:2183843/1‑70 (MQ=255)
                              tCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATgcga                                                                >  1:4232701/1‑45 (MQ=255)
                              tCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaaa                                        >  1:1816327/1‑70 (MQ=255)
                               cATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGaa                                                        <  1:672156/53‑1 (MQ=255)
                                   cAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGaa                                                        >  1:1969374/1‑49 (MQ=255)
                                   cAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCAt                                  >  1:2698364/1‑71 (MQ=255)
                                     ggTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTa                                <  1:77216/71‑1 (MQ=255)
                                         aCTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                             <  1:733302/70‑1 (MQ=255)
                                          cTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAg                                                          <  1:1363402/40‑1 (MQ=255)
                                           ttttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAg                                                          <  1:174498/39‑1 (MQ=255)
                                            tttCATGGCTATCATCCCAAAAAACCATGCGCGGTTa                                                           >  1:2197346/1‑37 (MQ=255)
                                            tttCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATc                          >  1:2459062/1‑70 (MQ=255)
                                             ttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTa                        >  1:12301/1‑71 (MQ=255)
                                                  ggCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                             >  1:2190100/1‑61 (MQ=255)
                                                   gCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAt                           >  1:354154/1‑62 (MQ=255)
                                                   gCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGt                   <  1:697617/70‑1 (MQ=255)
                                                    cTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTgg                 <  1:2700958/71‑1 (MQ=255)
                                                     tATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                             <  1:3815336/58‑1 (MQ=255)
                                                       tcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGt                   <  1:3835020/66‑1 (MQ=255)
                                                       tcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGt               >  1:54412/1‑70 (MQ=255)
                                                       tcatcCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAg                                     <  1:2290532/48‑1 (MQ=255)
                                                        catcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaaaa                                       >  1:3637423/1‑45 (MQ=255)
                                                         atcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaa                               >  1:2815250/1‑52 (MQ=255)
                                                         atcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAt                           >  1:334933/1‑56 (MQ=255)
                                                            ccAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCgg         >  1:1270309/1‑71 (MQ=255)
                                                              aaaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATc                          <  1:2230464/52‑1 (MQ=255)
                                                               aaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAg                                     <  1:3421931/40‑1 (MQ=255)
                                                               aaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCg      <  1:3969514/71‑1 (MQ=255)
                                                                aaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGc     >  1:3414034/1‑71 (MQ=255)
                                                                 aaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                             >  1:1527694/1‑46 (MQ=255)
                                                                  aaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTg   <  1:1010578/71‑1 (MQ=255)
                                                                   aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGt  <  1:3588822/71‑1 (MQ=255)
                                                                   aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGt  <  1:1276157/71‑1 (MQ=255)
                                                                   aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGt  <  1:216365/71‑1 (MQ=255)
                                                                   aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGt  >  1:1787017/1‑71 (MQ=255)
                                                                    cTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCt             >  1:2937988/1‑59 (MQ=255)
                                                                     |                                                                    
ATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGT  >  W3110S.gb/429781‑429918

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: