Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 429,850 T→C 47.2% 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,8500TC47.2% 1.8 / 75.6 53Y8H (TAT→CAT) yajDconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/14);  new base C (15/10);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGC  >  W3110S.gb/429793‑429915
                                                         |                                                                 
tGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAg                                                       >  1:1870039/1‑70 (MQ=255)
  cgcgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGa                                                      <  1:4316799/69‑1 (MQ=255)
  cgcgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGaa                                                     <  1:2828678/70‑1 (MQ=255)
   gcgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAg                                                   <  1:4234895/71‑1 (MQ=255)
             aCATTTCATCCAGGTAACTTTTCATCGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATcgc                                         >  1:1679996/1‑71 (MQ=255)
              cATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATcgc                                         >  1:807436/1‑70 (MQ=255)
               aTTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGaa                                                     <  1:951937/57‑1 (MQ=255)
                 ttCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaa                                      >  1:3070266/1‑70 (MQ=255)
                  tCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaaa                                     >  1:2512757/1‑70 (MQ=255)
                    aTCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAg                                  <  1:4974266/71‑1 (MQ=255)
                       cAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCAATCGCGAAAAAGCa                                >  1:1344234/1‑70 (MQ=255)
                         ggTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGt                                                          >  1:4718356/1‑42 (MQ=255)
                          gTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCa                                <  1:1874492/67‑1 (MQ=255)
                          gTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCa                                <  1:4790585/67‑1 (MQ=255)
                          gTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAg                                                       <  1:711315/44‑1 (MQ=255)
                          gTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTa                             >  1:1174988/1‑70 (MQ=255)
                          gTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTa                             >  1:3558196/1‑70 (MQ=255)
                           tAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaa                           <  1:3068294/71‑1 (MQ=255)
                            aaCTTTTCATGGCTATGATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaa                           >  1:2224367/1‑70 (MQ=255)
                            aaCTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGaa                                                     >  1:4541843/1‑44 (MQ=255)
                            aaCTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGaa                                                     >  1:881925/1‑44 (MQ=255)
                            aaCTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGaa                                                     >  1:115530/1‑44 (MQ=255)
                            aaCTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGaa                                                     >  1:1967657/1‑44 (MQ=255)
                            aaCTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGaa                                                     >  1:2118556/1‑44 (MQ=255)
                            aaCTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaa                           >  1:4476173/1‑70 (MQ=255)
                             aCTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                          <  1:515112/70‑1 (MQ=255)
                             aCTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                          <  1:1494255/70‑1 (MQ=255)
                             aCTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                          <  1:3755299/70‑1 (MQ=255)
                              cTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAt                        >  1:1373591/1‑71 (MQ=255)
                               ttttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATc                       >  1:2492415/1‑71 (MQ=255)
                               ttttCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAt                        >  1:2373161/1‑70 (MQ=255)
                                tttCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGaa                                                     >  1:4960112/1‑40 (MQ=255)
                                 ttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaaa                                     >  1:3331134/1‑55 (MQ=255)
                                 ttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaa                                      >  1:1005188/1‑54 (MQ=255)
                                     tGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAg                                  <  1:3444037/54‑1 (MQ=255)
                                     tGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCg                 <  1:388764/71‑1 (MQ=255)
                                      ggCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGt                >  1:2052414/1‑71 (MQ=255)
                                        cTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTgg              <  1:2541812/71‑1 (MQ=255)
                                         tATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCt                                              <  1:2436813/38‑1 (MQ=255)
                                          atcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTggg             >  1:697258/1‑70 (MQ=255)
                                            catcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCt          >  1:4070478/1‑71 (MQ=255)
                                             atcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATc                   <  1:3733130/61‑1 (MQ=255)
                                             atcCCAAAAAACCAT‑‑GCGGTTAGAAAGCGGCTATCGCGAAAAAg                                  >  1:2200935/1‑44 (MQ=37)
                                              tcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaa                           <  1:832712/52‑1 (MQ=255)
                                              tcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaa                           <  1:4343022/52‑1 (MQ=255)
                                               cccAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCg                 >  1:4645064/1‑61 (MQ=255)
                                                 cAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                          <  1:1863127/50‑1 (MQ=255)
                                                 cAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa                          <  1:1343516/50‑1 (MQ=255)
                                                   aaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTc    <  1:2567007/70‑1 (MQ=255)
                                                   aaaaaCCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCAtt                              >  1:2953596/1‑44 (MQ=255)
                                                   aaaaaCCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTgg              >  1:4579098/1‑60 (MQ=255)
                                                   aaaaaCCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCt          <  1:3901221/64‑1 (MQ=255)
                                                    aaaaCCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGc  <  1:4160786/71‑1 (MQ=255)
                                                         |                                                                 
TGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGC  >  W3110S.gb/429793‑429915

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: