Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,335,995 C→T 21.9% pseudogene (3078/4605 nt)
pseudogene (3078/4104 nt)
yfaS ←
yfaS ←
ECK2220:JW2222+JW2221:b4500; hypothetical protein
ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,335,9950CT21.9% 64.5 / 13.8 32pseudogene (3078/4605 nt)
pseudogene (3078/4104 nt)
yfaS
yfaS
ECK2220:JW2222+JW2221:b4500; hypothetical protein
ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment
Reads supporting (aligned to +/- strand):  ref base C (20/5);  new base T (3/4);  total (23/9)
Fisher's exact test for biased strand distribution p-value = 7.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGT  >  W3110S.gb/2335927‑2336062
                                                                    |                                                                   
tCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCAt                                                                   >  1:590533/1‑71 (MQ=255)
  actactTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGa                                                                 >  1:1140610/1‑71 (MQ=255)
                    atTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACATTTTTCTTCTGTTc                                                 <  1:362603/69‑1 (MQ=255)
                    atTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCaa                                               >  1:1734834/1‑71 (MQ=255)
                     ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACATTTTTCTGCTGTTc                                                 <  1:2174786/68‑1 (MQ=255)
                     ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACATTTTTCTGCTGTTc                                                 <  1:60908/68‑1 (MQ=255)
                     ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAAc                                              >  1:1237497/1‑71 (MQ=255)
                      tGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCaa                                               >  1:1209894/1‑69 (MQ=255)
                       gCTCGGCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACATTTTTctgc                                                      >  1:1736067/1‑61 (MQ=255)
                       gCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACgg                                            >  1:2752998/1‑71 (MQ=255)
                       gCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACgg                                            >  1:144327/1‑71 (MQ=255)
                        cTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGTATGACATTTTTCTGCTGTTCAACgg                                            <  1:380001/70‑1 (MQ=255)
                              tGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACgg                                            <  1:872644/64‑1 (MQ=255)
                                    ggCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATc                               >  1:1925949/1‑71 (MQ=255)
                                         cATCAGCGCGTTATCGCCACCACCGAGCATGATATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCcaca                           >  1:2623125/1‑70 (MQ=255)
                                          aTCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCaa                                               >  1:2627200/1‑49 (MQ=255)
                                          aTCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACaa                          >  1:676540/1‑70 (MQ=255)
                                             aGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCaa                                               >  1:67286/1‑46 (MQ=255)
                                             aGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCaa                                               >  1:796301/1‑46 (MQ=255)
                                             aGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGa                      >  1:1031309/1‑71 (MQ=255)
                                              gcgcGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGAc                     <  1:1585039/71‑1 (MQ=255)
                                                     aTCGCCACCACCGAGTATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTtg              >  1:2409432/1‑71 (MQ=255)
                                                     aTCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAAc                                              >  1:1178255/1‑39 (MQ=255)
                                                     aTCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATc                               <  1:1096261/54‑1 (MQ=255)
                                                     aTCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTtg              >  1:2009063/1‑71 (MQ=255)
                                                     aTCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTtg              >  1:1697181/1‑71 (MQ=255)
                                                        gCCACCACCGAGCATGACATTTTTCTGCTGTTCaa                                               >  1:1173232/1‑35 (MQ=255)
                                                        gCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGc                                          >  1:974684/1‑40 (MQ=255)
                                                         ccaccaCCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCt          <  1:1551733/71‑1 (MQ=255)
                                                            ccaccGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAAt       <  1:321181/71‑1 (MQ=255)
                                                                 gAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGt  >  1:671739/1‑71 (MQ=255)
                                                                   gCATGACATTTTTCTGCTGTTCAACGGGCCAGGTGTTAt                                >  1:157805/1‑39 (MQ=255)
                                                                    |                                                                   
TCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGT  >  W3110S.gb/2335927‑2336062

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