Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,721,565 T→C 25.0% intergenic (‑249/‑25) ydhH ← / → slyB conserved hypothetical protein/outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,721,5650TC25.0% 46.0 / 14.6 28intergenic (‑249/‑25)ydhH/slyBconserved hypothetical protein/outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (13/8);  new base C (4/3);  total (17/11)
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

GTAGAACATGTAACCATGGATTCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCAATGGTAGGTCTGTC  >  W3110S.gb/1721498‑1721630
                                                                   |                                                                 
gTAGAACATGTAACCATGGATTCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTc                                                                 <  1:1634499/70‑1 (MQ=255)
      cATGTAACCATGGATTCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTg                                                           >  1:1347073/1‑70 (MQ=255)
       aTGTAACCATGGATTCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGaa                                                         >  1:1062585/1‑71 (MQ=255)
          tAACCATGGATTCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTCTCATGGTGAACAg                                                      <  1:167263/71‑1 (MQ=255)
                  gATTCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATg                                                              >  1:428156/1‑55 (MQ=255)
                  gATTCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTCTCATGGTGAACAGGAGAtttt                                              >  1:1725341/1‑71 (MQ=255)
                     tCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCa                                                                <  1:213122/50‑1 (MQ=255)
                            tGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCaa                                          <  1:895178/65‑1 (MQ=255)
                            tGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGAtt                                     >  1:78837/1‑70 (MQ=255)
                              ccATATACTTTGACCATGAGAGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTaaa                                  <  1:634075/71‑1 (MQ=255)
                                atatACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACa                                                       >  1:913490/1‑48 (MQ=255)
                                atatACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAAcg                                >  1:1412211/1‑71 (MQ=255)
                                 tataCTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATgg                                                             >  1:1516739/1‑41 (MQ=255)
                                 tataCTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGAtt                                     <  1:1492573/65‑1 (MQ=255)
                                  ataCTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAAcg                                <  1:878076/69‑1 (MQ=255)
                                   taCTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAAcgcg                              >  1:1511810/1‑70 (MQ=255)
                                          aCCATGAGGGATGCTTGCGTGGCGTCTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGgtt                      <  1:299113/71‑1 (MQ=255)
                                              tGAGGGATGCTTGCGTGGCGTCTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGttt                  >  1:570867/1‑71 (MQ=255)
                                              tGAGGGATGCTTGCGTGGCGTCTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGttt                  >  1:1218772/1‑71 (MQ=255)
                                                 gggATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAATGCGTATTGGTTGTTTCaa               >  1:1297252/1‑71 (MQ=255)
                                                 gggATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCaa               >  1:221143/1‑71 (MQ=255)
                                                     tGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCAATGGt           <  1:340297/71‑1 (MQ=255)
                                                       cTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGgttgt                    <  1:267463/60‑1 (MQ=255)
                                                          gCGTGGCGTCTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCa                <  1:432470/61‑1 (MQ=255)
                                                              ggCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCAATGGTAGgtctgtc  >  1:242723/1‑71 (MQ=255)
                                                               gCGTCTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCAATGGTAGgtct     >  1:1618176/1‑67 (MQ=255)
                                                                cGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCAATGGTAGgt       >  1:1671012/1‑64 (MQ=255)
                                                                cGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCAATGGTAGgt       >  1:942106/1‑64 (MQ=255)
                                                                   |                                                                 
GTAGAACATGTAACCATGGATTCACATATGCCATATACTTTGACCATGAGGGATGCTTGCGTGGCGTTTCATGGTGAACAGGAGATTTTTCAATGATTAAACGCGTATTGGTTGTTTCAATGGTAGGTCTGTC  >  W3110S.gb/1721498‑1721630

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