Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,350,232 G→A 52.0% R3C (CGC→TGC)  yjdF ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,350,2320GA52.0% ‑4.7 / 33.6 25R3C (CGC→TGC) yjdFconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base A (8/5);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAA  >  W3110S.gb/4350176‑4350296
                                                        |                                                                
cAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGtt                                                    >  1:1623571/1‑71 (MQ=255)
cAGGATTAAAGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGt                                                     <  1:134300/70‑1 (MQ=255)
    aTTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCa                                                <  1:584588/71‑1 (MQ=255)
     ttAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCa                                                >  1:409078/1‑70 (MQ=255)
       aaCGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAAt                                             >  1:751856/1‑71 (MQ=255)
       aaCGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAAt                                             >  1:204823/1‑71 (MQ=255)
        aCGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAATc                                            <  1:992912/71‑1 (MQ=255)
         cGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCa                                                <  1:741176/66‑1 (MQ=255)
             aGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCagg                                                                >  1:802036/1‑46 (MQ=255)
              gCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAATCAGGGa                                       >  1:1072377/1‑70 (MQ=255)
                gTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGaaaa                                    <  1:777623/71‑1 (MQ=255)
                 tCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAATCAGGGaaaa                                    <  1:979503/70‑1 (MQ=255)
                     tGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAATc                                            <  1:228607/58‑1 (MQ=255)
                     tGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAATc                                            <  1:508227/58‑1 (MQ=255)
                        gcTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGAc                             >  1:698249/1‑70 (MQ=255)
                            gtgtTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTAtt                         <  1:1148467/70‑1 (MQ=255)
                               ttAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTg                                                        <  1:755022/36‑1 (MQ=255)
                               ttAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGAc                             <  1:1364310/63‑1 (MQ=255)
                               ttAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGAc                             <  1:766196/63‑1 (MQ=255)
                                 aaGAATTAACGGCTTGAGTGTGCGGGTCATGGTGa                                                       >  1:1546183/1‑35 (MQ=255)
                                 aaGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAg                   >  1:1272841/1‑71 (MQ=255)
                                    aaTTAACGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAATCAGGGaaa                                     <  1:560224/50‑1 (MQ=255)
                                        aaCGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCg             >  1:247819/1‑70 (MQ=255)
                                        aaCGGCTTGAGTGTGCAGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGt            >  1:1440394/1‑71 (MQ=255)
                                            gCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATa                                  >  1:1443272/1‑45 (MQ=255)
                                                     tgCAGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACaa  >  1:295329/1‑68 (MQ=255)
                                                        |                                                                
CAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAA  >  W3110S.gb/4350176‑4350296

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: