Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 530,077 T→C 15.4% 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,0770TC15.4% 126.2 / 16.7 52Y125C (TAT→TGT) ybbBtRNA 2‑selenouridine synthase, selenophosphate‑dependent
Reads supporting (aligned to +/- strand):  ref base T (25/19);  new base C (5/3);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

AATCAGCACTATCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATATGTGAGCGCTG  >  W3110S.gb/530010‑530147
                                                                   |                                                                      
aaTCAGCACTATCGGTTTTTGTGCCAGTTCAATAGTCGCATGAATCGCGGTCTGGCGCAGTGCCTtat                                                                        >  1:4360056/1‑68 (MQ=255)
     gCACTATCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCt                                                                <  1:4877436/71‑1 (MQ=255)
     gCACTATCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCt                                                                <  1:168656/71‑1 (MQ=255)
      cACTATCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCtt                                                               >  1:1529282/1‑71 (MQ=255)
           tCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTccac                                                           >  1:4785151/1‑70 (MQ=255)
             ggTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCTCAGTGCCTTATAACCGCCTTCCACCAg                                                        >  1:821643/1‑71 (MQ=255)
             ggTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAg                                                        <  1:1485498/71‑1 (MQ=255)
             ggTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAg                                                        >  1:1245429/1‑71 (MQ=255)
                       ccAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatca                                              <  1:2769546/71‑1 (MQ=255)
                       ccAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatca                                              <  1:1708128/71‑1 (MQ=255)
                       ccAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTACACCAGCGGATaatc                                               <  1:2578044/70‑1 (MQ=255)
                         aGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATaatcaat                                            <  1:4340643/71‑1 (MQ=255)
                         aGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATaatcaat                                            <  1:4334781/71‑1 (MQ=255)
                            taaaTAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCACCTTCCACCAGCGGATAATCAATccc                                         >  1:369493/3‑71 (MQ=255)
                            tCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATccc                                         >  1:615316/1‑71 (MQ=255)
                             cAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatcaa                                             >  1:2125759/1‑66 (MQ=255)
                             cAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatcaa                                             >  1:4273866/1‑66 (MQ=255)
                             cAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatcaa                                             >  1:165282/1‑66 (MQ=255)
                             cAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTAATAACCGCCTTCCACCAGCGGATAATCAATcccc                                        >  1:4953116/1‑71 (MQ=255)
                             cAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTAATAACCGCCTTCCACCAGCGGATAATCAATcccc                                        >  1:2181231/1‑71 (MQ=255)
                             cAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTAATAACCGCCTTCCACCAGCGGATAATCAATccc                                         >  1:3086884/1‑70 (MQ=255)
                              aaTAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCg                                       <  1:3630935/71‑1 (MQ=255)
                                tAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCg                                       >  1:3200569/1‑69 (MQ=255)
                                tAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCg                                       >  1:2168653/1‑69 (MQ=255)
                                 aGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaa                                                 >  1:285343/1‑58 (MQ=255)
                                 aGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaa                                                 >  1:3773060/1‑58 (MQ=255)
                                     gCCTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATAATCAATcccc                                        >  1:1442116/1‑63 (MQ=255)
                                       cTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatca                                              <  1:4278371/55‑1 (MQ=255)
                                       cTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCtgcatgca                              <  1:2130751/71‑1 (MQ=255)
                                       cTGAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATAATCAATCCCCGCtgcatgca                              <  1:4604425/71‑1 (MQ=255)
                                         gAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGc                                                       >  1:2026557/1‑44 (MQ=255)
                                         gAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCg                                                      >  1:585776/1‑45 (MQ=255)
                                         gAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATaa                                                 >  1:5035042/1‑50 (MQ=255)
                                         gAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATaa                                                 >  1:248353/1‑50 (MQ=255)
                                         gAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATaa                                                 >  1:2709256/1‑50 (MQ=255)
                                         gAATCGCGGTCTGGCGCAGTGCCTTACAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAAc                            >  1:207819/1‑71 (MQ=255)
                                           aTCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatca                                              <  1:1881911/51‑1 (MQ=255)
                                           aTCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatca                                              <  1:360941/51‑1 (MQ=255)
                                             cgcgGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatcaat                                            >  1:244613/1‑51 (MQ=255)
                                             cgcgGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatca                                              >  1:2248325/1‑49 (MQ=255)
                                             cgcgGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATaatca                                              >  1:331414/1‑49 (MQ=255)
                                                      gcgcAGTGCCTTATAACCGCCTTCCACCAGCGGATaa                                                 >  1:347533/1‑37 (MQ=255)
                                                      gcgcAGTGCCTTATAACCGCCTTCCACCAGCGGATaa                                                 >  1:5068301/1‑37 (MQ=255)
                                                      gcgcAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGca                  <  1:2592699/68‑1 (MQ=255)
                                                      gcgcAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGca                  <  1:4449566/68‑1 (MQ=255)
                                                       cgcAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAat              <  1:391659/71‑1 (MQ=255)
                                                            tGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGcac                 <  1:2421970/63‑1 (MQ=255)
                                                            tGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGcac                 <  1:4078586/63‑1 (MQ=255)
                                                            tGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATAtg           <  1:712633/69‑1 (MQ=255)
                                                            tGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATAtg           <  1:284474/69‑1 (MQ=255)
                                                             gCCTTATAACCGCCTTCCACCAGCGGATAATCAATCACCGCTGCATGCAACCAGCttt                     >  1:383452/1‑58 (MQ=255)
                                                              ccTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCtgca                                  >  1:383236/1‑44 (MQ=255)
                                                                ttATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAg                         <  1:4502296/51‑1 (MQ=255)
                                                                  ataACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATAtg           >  1:1078189/1‑63 (MQ=255)
                                                                  ataACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATAtg           >  1:375011/1‑63 (MQ=255)
                                                                  ataACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATATGTGAGCGCt   >  1:4652273/1‑71 (MQ=255)
                                                                   taACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATATGTGAGCGCt   >  1:3938337/1‑70 (MQ=255)
                                                                   taACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATATGTGAGCGCTg  <  1:1024450/71‑1 (MQ=255)
                                                                   taACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATATGTGAGCGCTg  <  1:935175/71‑1 (MQ=255)
                                                                   |                                                                      
AATCAGCACTATCGGTTTTTGTGCCAGTTCAATAGTCGCCTGAATCGCGGTCTGGCGCAGTGCCTTATAACCGCCTTCCACCAGCGGATAATCAATCCCCGCTGCATGCAACCAGCTTTGCACAATATGTGAGCGCTG  >  W3110S.gb/530010‑530147

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