Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,368,720 T→C 10.4% N232S (AAC→AGC)  dusC ← tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,368,7200TC10.4% 126.7 / 6.2 48N232S (AAC→AGC) dusCtRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base C (2/3);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

CAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCG  >  Exported/1368659‑1368789
                                                             |                                                                     
caaagcaaCCACCTCCGTCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAatat                                                              <  1:4645443/71‑1 (MQ=255)
 aaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATAtt                                                             <  1:3566500/71‑1 (MQ=255)
 aaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATAtt                                                             <  1:3560962/71‑1 (MQ=255)
 aaagcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGCTGGGAATAtt                                                             <  1:3550892/71‑1 (MQ=255)
    gcaaCCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAg                                                          >  1:615858/1‑71 (MQ=255)
       accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAgcgc                                                       >  1:1686451/1‑71 (MQ=255)
       accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAgcgc                                                       >  1:609355/1‑71 (MQ=255)
       accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAgcgc                                                       >  1:3085608/1‑71 (MQ=255)
       accaccTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAgcgc                                                       >  1:4170527/1‑71 (MQ=255)
              ccGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGa                                                >  1:3466078/1‑71 (MQ=255)
               cGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGAc                                               <  1:3153817/71‑1 (MQ=255)
               cGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGAc                                               <  1:2893738/71‑1 (MQ=255)
                    aCGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGcc                                                      <  1:1224757/59‑1 (MQ=255)
                      ggCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAatcat                                        >  1:4613431/1‑71 (MQ=255)
                        cATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGaa                                                                 <  1:3192622/44‑1 (MQ=255)
                          ttCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGaa                                                                 <  1:369348/42‑1 (MQ=255)
                            cgcgGTTCGTTATATTTTACCACCCGGCTCAGGTTggg                                                                   >  1:20790/1‑38 (MQ=255)
                            cgcgGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCa                                             >  1:4567984/1‑60 (MQ=255)
                            cgcgGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCa                                             >  1:1912201/1‑60 (MQ=255)
                            cgcgGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCACGCGACCAATCATCACTg                                   >  1:2841584/1‑70 (MQ=255)
                                  tCGTTATATTTTACCACCCGGCTCAGGCTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCg                              >  1:3916813/1‑69 (MQ=255)
                                  tCGTTATATTTTACCACCCGGCTCAGGCTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCg                              >  1:652390/1‑69 (MQ=255)
                                    gTTATATTTTACCACCCGGCTCAGGCTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGc                          <  1:3255524/71‑1 (MQ=255)
                                     ttATATTTTACCACCCGGCTCAGGCTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGc                          <  1:1239040/70‑1 (MQ=255)
                                       atatTTTACCACCCGGCTCAGTTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGc                       >  1:955076/1‑71 (MQ=255)
                                       atatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGc                       >  1:3707906/1‑71 (MQ=255)
                                       atatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGc                       >  1:4670018/1‑71 (MQ=255)
                                       atatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGc                       >  1:4704090/1‑71 (MQ=255)
                                       atatTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGc                       >  1:1111625/1‑71 (MQ=255)
                                          ttttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTg                     <  1:3589003/70‑1 (MQ=255)
                                          ttttACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGa                    >  1:3273622/1‑71 (MQ=255)
                                              accaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCg                 <  1:839636/70‑1 (MQ=255)
                                               ccaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGcc               <  1:2013356/71‑1 (MQ=255)
                                                   ccGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATg            <  1:2984417/70‑1 (MQ=255)
                                                     ggCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCAt         <  1:1129491/71‑1 (MQ=255)
                                                     ggCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCAt         <  1:1708137/71‑1 (MQ=255)
                                                     ggCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCAt         <  1:173425/71‑1 (MQ=255)
                                                     ggCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCAt         <  1:2459821/71‑1 (MQ=255)
                                                      gCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCAt         <  1:2609386/70‑1 (MQ=255)
                                                       cTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCAttg       <  1:2572924/71‑1 (MQ=255)
                                                       cTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCAttg       <  1:635157/71‑1 (MQ=255)
                                                          aGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGc                                  <  1:2789551/41‑1 (MQ=255)
                                                          aGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGc                                  <  1:3047227/41‑1 (MQ=255)
                                                            gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGt                                >  1:2786194/1‑41 (MQ=255)
                                                            gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGt                                >  1:4273668/1‑41 (MQ=255)
                                                            gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGt                                >  1:4111821/1‑41 (MQ=255)
                                                            gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGt                                >  1:945770/1‑41 (MQ=255)
                                                            gTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTgcg  >  1:1294738/1‑71 (MQ=255)
                                                             |                                                                     
CAAAGCAACCACCTCCGGCCACGGCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCG  >  Exported/1368659‑1368789

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: