Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 530,077 T→C 55.0% Y125C (TAT→TGT)  ybbB ← tRNA 2‑selenouridine synthase, selenophosphate‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb530,0770TC55.0% 3.4 / 22.1 20Y125C (TAT→TGT) ybbBtRNA 2‑selenouridine synthase, selenophosphate‑dependent
Reads supporting (aligned to +/- strand):  ref base T (4/5);  new base C (6/5);  total (10/10)
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

AATCAGCACTATCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATATGT  >  W3110S.gb/530010‑530139
                                                                   |                                                              
aaTCAGCACTATCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTACAAc                                                             >  1:1356031/1‑71 (MQ=255)
          aTCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTccac                                                   >  1:774984/1‑71 (MQ=255)
          aTCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTccac                                                   >  1:664428/1‑71 (MQ=255)
           tCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTccacc                                                  <  1:983925/71‑1 (MQ=255)
                 tttGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGa                                            >  1:527342/1‑71 (MQ=255)
                  ttGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGa                                            <  1:57896/70‑1 (MQ=255)
                              aaTAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATAATCAATcccc                                <  1:1587016/70‑1 (MQ=255)
                                tAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTtat                                                                <  1:81039/36‑1 (MQ=255)
                                tAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATAATCAATCCCCGCt                             >  1:1101007/1‑71 (MQ=255)
                                   tCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCg                                                           <  1:942119/38‑1 (MQ=255)
                                   tCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaat                                        <  1:1274922/57‑1 (MQ=255)
                                     gCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCtgcatg                        >  1:1324409/1‑71 (MQ=255)
                                     gCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGc                                               >  1:1028308/1‑48 (MQ=255)
                                     gCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGa                                            >  1:1285778/1‑51 (MQ=255)
                                        tGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCaa                     >  1:1283951/1‑71 (MQ=255)
                                        tGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATaatca                                      <  1:512688/54‑1 (MQ=255)
                                               cgGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCtt              <  1:651068/71‑1 (MQ=255)
                                                ggTCTGGCGCAGTGCATTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCttt             <  1:1270419/71‑1 (MQ=255)
                                                 gTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCttt             <  1:745494/70‑1 (MQ=255)
                                                        gcAGTGCCTTACAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAata     <  1:452506/71‑1 (MQ=255)
                                                             gCCTTACAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATAtgt  >  1:1611924/1‑69 (MQ=255)
                                                                 tataACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAata     <  1:1712723/62‑1 (MQ=255)
                                                                   |                                                              
AATCAGCACTATCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATATGT  >  W3110S.gb/530010‑530139

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: