Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,225,948 G→C 69.2% V129V (GTG→GTC dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,225,9480GC69.2% 99.7 / 61.3 78V129V (GTG→GTCdusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base G (8/16);  new base C (25/29);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAC  >  Exported/2225887‑2226011
                                                             |                                                               
tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATg                                                             >  1:1665820/1‑66 (MQ=255)
tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATg                                                             >  1:221103/1‑66 (MQ=255)
       ttCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTc                                                 <  1:1925402/71‑1 (MQ=255)
          aGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCt                                               <  1:2927601/70‑1 (MQ=255)
              ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTAc                                          >  1:4624023/1‑71 (MQ=255)
              ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTAc                                          >  1:1336361/1‑71 (MQ=255)
              ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTAc                                          >  1:5237067/1‑71 (MQ=255)
              ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAc                                          >  1:3468253/1‑71 (MQ=255)
                ctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAccc                                        >  1:1295283/1‑71 (MQ=255)
                 ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGt                                             >  1:1795504/1‑65 (MQ=255)
                   gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATg                                                             >  1:40647/1‑47 (MQ=255)
                   gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATg                                                             >  1:5041518/1‑47 (MQ=255)
                   gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAg                                                          >  1:5185048/1‑50 (MQ=255)
                   gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAg                                                          >  1:629643/1‑50 (MQ=255)
                   gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAg                                                          >  1:5028377/1‑50 (MQ=255)
                    cAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGt                                             >  1:3917017/1‑62 (MQ=255)
                    cAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTg                                      <  1:2424299/69‑1 (MQ=255)
                    cAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTg                                      >  1:5204703/1‑69 (MQ=255)
                     aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAg                                   <  1:2090741/71‑1 (MQ=255)
                      gTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                    <  1:5109469/69‑1 (MQ=255)
                       tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGAt                                 <  1:3056273/71‑1 (MQ=255)
                       tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGAt                                 <  1:1958653/71‑1 (MQ=255)
                       tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGAt                                 <  1:124103/71‑1 (MQ=255)
                       tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGAt                                 <  1:5350044/71‑1 (MQ=255)
                       tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGAt                                 <  1:4896597/71‑1 (MQ=255)
                       tACACGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGAt                                 <  1:2655143/71‑1 (MQ=255)
                         cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGaa                                                               >  1:4654059/1‑39 (MQ=255)
                         cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc                               >  1:2845667/1‑71 (MQ=255)
                           cGGATGTCGTTAAATCGATCGTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGc                             <  1:2644833/71‑1 (MQ=255)
                           cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGc                             <  1:4919816/71‑1 (MQ=255)
                            ggATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCa                            >  1:3008416/1‑71 (MQ=255)
                             gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGc                                                            <  1:2326509/38‑1 (MQ=255)
                             gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGc                                                            <  1:4220195/38‑1 (MQ=255)
                             gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGc                                                            <  1:4403067/38‑1 (MQ=255)
                             gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGc                                                            <  1:230534/38‑1 (MQ=255)
                               tGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAg                                                          >  1:3374263/1‑38 (MQ=255)
                               tGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                    <  1:858369/60‑1 (MQ=255)
                               tGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                    <  1:2089286/60‑1 (MQ=255)
                               tGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg                         <  1:3186101/71‑1 (MQ=255)
                                gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCt                                       <  1:3392760/56‑1 (MQ=255)
                                   gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                     >  1:3284821/1‑71 (MQ=255)
                                   gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                     >  1:3046298/1‑71 (MQ=255)
                                   gTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                     >  1:3761167/1‑71 (MQ=255)
                                    ttAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAcc                          <  1:2714541/65‑1 (MQ=255)
                                    ttAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAcc                          <  1:2482250/65‑1 (MQ=255)
                                    ttAAATCGATCCTTACCGAGGTCGTCAATGCAGTg                                                        <  1:4924703/35‑1 (MQ=255)
                                    ttAAATCGATCCTTACCGAGGTCGTCAATGCAGTg                                                        <  1:3240402/35‑1 (MQ=255)
                                    ttAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                    <  1:2021558/55‑1 (MQ=255)
                                    ttAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTgg                    <  1:4939939/71‑1 (MQ=255)
                                        aTCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGTCAc                >  1:3956597/1‑71 (MQ=255)
                                        aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTgg                    <  1:5102724/67‑1 (MQ=255)
                                        aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTgg                    <  1:1868592/67‑1 (MQ=255)
                                         tCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTg                                      <  1:3455423/48‑1 (MQ=255)
                                         tCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTg                                      <  1:2284747/48‑1 (MQ=255)
                                          cGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                    >  1:243965/1‑49 (MQ=255)
                                          cGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                    >  1:3241311/1‑49 (MQ=255)
                                          cGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                    >  1:431998/1‑49 (MQ=255)
                                          cGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                    <  1:3178624/49‑1 (MQ=255)
                                           gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg              <  1:3120779/70‑1 (MQ=255)
                                           gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg              <  1:1077812/70‑1 (MQ=255)
                                           gATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTa                                           >  1:131723/1‑41 (MQ=255)
                                           gATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg             <  1:2928151/71‑1 (MQ=255)
                                           gATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg             <  1:680092/71‑1 (MQ=255)
                                               cTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                     <  1:4969543/59‑1 (MQ=255)
                                               cTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                     <  1:3456884/59‑1 (MQ=255)
                                               cTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                     <  1:5139341/59‑1 (MQ=255)
                                                 tACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                   >  1:275314/1‑59 (MQ=255)
                                                 tACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTa                                           >  1:3221800/1‑35 (MQ=255)
                                                 tACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACAcc       >  1:5344348/1‑71 (MQ=255)
                                                 tACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACAcc       >  1:5681696/1‑71 (MQ=255)
                                                  aCCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca         <  1:3341946/68‑1 (MQ=255)
                                                  aCCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca         <  1:5170299/68‑1 (MQ=255)
                                                    cGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTa    <  1:5143346/71‑1 (MQ=255)
                                                     gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                   <  1:3464707/55‑1 (MQ=255)
                                                     gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                   <  1:1480544/55‑1 (MQ=255)
                                                     gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTa    <  1:750231/70‑1 (MQ=255)
                                                      aGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAc  >  1:41939/1‑71 (MQ=255)
                                                      aGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAc  >  1:2499402/1‑71 (MQ=255)
                                                             gAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGaa           >  1:3704718/1‑55 (MQ=255)
                                                             gAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGaa           >  1:1101639/1‑55 (MQ=255)
                                                             |                                                               
TCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAC  >  Exported/2225887‑2226011

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: