Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,944,962 T→C pseudogene (507/729 nt) yhjQ → ECK3519:JW5941:b3534; cell division protein; chromosome partitioning ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,944,9620TC99.8% 1452.2 / NA 623pseudogene (507/729 nt)yhjQECK3519:JW5941:b3534; cell division protein; chromosome partitioning ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (0/621);  minor base G (0/1);  total (0/623)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 623 total aligned reads displayed.

AGCGCTGCCGGATGGCGCACATATTTTGATTAATGACTTCCGTATTGGCAGTCAGGTTCAGGA  >  W3110S.gb/3944911‑3944973
                                                   |           
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTTGCAGCCAGGTTCAGGa  <  1:886603/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCGGCCAGGTTCAGGa  <  1:688750/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:61266/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:729687/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:720083/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:703405/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:695083/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:673612/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:666842/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:659514/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:646139/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:630218/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:625909/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:742703/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:605114/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:585137/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:577866/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:568040/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:548854/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:537646/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:533314/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:528252/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:519124/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:1000038/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:499747/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:891684/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:993471/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:987316/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:975847/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:969219/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:9601/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:952423/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:940136/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:924585/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:914136/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:90815/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:899827/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:735879/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:874739/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:841017/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:830710/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:820770/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:806183/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:792375/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:773155/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:763520/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:757368/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:751995/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:5105/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:164299/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:285057/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:277530/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:26352/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:245827/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:233697/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:226617/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:220246/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:210761/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:195571/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:183987/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:172919/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:479184/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:15189/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:148944/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:14239/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:135962/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:130587/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:117848/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:107672/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:1033091/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:1022310/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:1010554/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:373621/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:483735/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:469352/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:466337/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:460380/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:45129/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:443728/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:437409/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:420406/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:412136/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:397491/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:390073/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:383086/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:298294/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:3634/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:355892/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:352128/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:339603/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:325041/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:315371/62‑1 (MQ=255)
 gcgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:309079/62‑1 (MQ=255)
 gcgcCGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:780633/62‑1 (MQ=255)
  cgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:862023/61‑1 (MQ=255)
  cgcTGCCGGATGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:851795/61‑1 (MQ=255)
         ggATGGTGCACATATTTTGATTAATTACTTCCTTATTGGCAGCCAGGTTCAGGa  <  1:254021/54‑1 (MQ=37)
           aTGGTGCACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:272448/52‑1 (MQ=39)
                gcACATATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:32970/47‑1 (MQ=255)
                   caTATTTTGATTAATGACTTCCGTATTGGCAGCCAGGTTCAGGa  <  1:15682/44‑1 (MQ=255)
                                                   |           
AGCGCTGCCGGATGGCGCACATATTTTGATTAATGACTTCCGTATTGGCAGTCAGGTTCAGGA  >  W3110S.gb/3944911‑3944973

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: