Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,247,298 A→G D238D (GAT→GAC truB ← tRNA pseudouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,247,2980AG91.4% 204.4 / 8.3 58D238D (GAT→GACtruBtRNA pseudouridine synthase B
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (31/22);  total (34/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 = 8.89e-01

CAGAAGACGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAGGGT  >  REL606/3247228‑3247369
                                                                      |                                                                       
cAGAAGACGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCaa                                                                     <  2:254554/75‑1 (MQ=255)
  gaagtCGTTAACGGAAGATTCACCCCCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGt                                                                   <  1:1190531/75‑1 (MQ=255)
   aagaCGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTa                                                                  <  2:695308/75‑1 (MQ=255)
   aagaCGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTa                                                                  <  2:631535/75‑1 (MQ=255)
   aagaCGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTa                                                                  <  2:230969/75‑1 (MQ=255)
     gaCGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAAc                                                                <  2:582591/75‑1 (MQ=255)
       cGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTc                                                              <  1:634568/75‑1 (MQ=255)
          tAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCcgc                                                           >  1:1079553/1‑75 (MQ=255)
          tAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCcgc                                                           <  2:590451/75‑1 (MQ=255)
                  gATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGaa                                                   <  2:1066653/75‑1 (MQ=255)
                  gATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGaa                                                   >  1:388237/1‑75 (MQ=255)
                   aTTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAat                                                  >  2:923088/1‑75 (MQ=255)
                   aTTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAat                                                  <  2:545315/75‑1 (MQ=255)
                     tCACCGCCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAatat                                                <  2:115187/75‑1 (MQ=255)
                      caccacCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATc                                               >  2:478479/1‑75 (MQ=255)
                       accaccGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATcc                                              >  1:1350044/1‑75 (MQ=255)
                        ccaccGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCct                                             <  2:516673/75‑1 (MQ=255)
                         caccGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCctg                                            >  2:1285882/1‑75 (MQ=255)
                         caccGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAACATcccg                                            >  2:396746/1‑73 (MQ=255)
                          accGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCctgc                                           <  1:3261/75‑1 (MQ=255)
                             ggATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATTGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGt                                        >  1:679031/1‑75 (MQ=255)
                             ggATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGt                                        >  1:1469817/1‑75 (MQ=255)
                              gATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGtt                                       <  1:1418211/75‑1 (MQ=255)
                                tAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCa                                     <  2:1378971/75‑1 (MQ=255)
                                tAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCa                                     <  2:373530/75‑1 (MQ=255)
                                tAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCa                                     >  1:739969/1‑75 (MQ=255)
                                 aGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAg                                    >  2:221993/1‑75 (MQ=255)
                                   tCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCt                                  >  2:1466258/1‑75 (MQ=255)
                                        aGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTc                             >  2:1244223/1‑75 (MQ=255)
                                           tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac                          >  1:772107/1‑75 (MQ=255)
                                           tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac                          >  1:1112390/1‑75 (MQ=255)
                                           tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac                          >  2:1578446/1‑75 (MQ=255)
                                           tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac                          <  1:408652/75‑1 (MQ=255)
                                           tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac                          <  1:408637/75‑1 (MQ=255)
                                           tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac                          >  1:206029/1‑75 (MQ=255)
                                               cTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGt                      <  2:1585565/75‑1 (MQ=255)
                                               cTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGt                      >  1:766074/1‑75 (MQ=255)
                                                 gTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTc                    <  2:1586580/75‑1 (MQ=255)
                                                 gTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTc                    <  1:32538/75‑1 (MQ=255)
                                                  tCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCa                   >  2:1391892/1‑75 (MQ=255)
                                                   ccATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCAc                  <  1:1511923/75‑1 (MQ=255)
                                                   ccATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCAc                  >  2:416591/1‑75 (MQ=255)
                                                    cATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACg                 >  1:1590772/1‑75 (MQ=255)
                                                     aTTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGc                <  2:109749/75‑1 (MQ=255)
                                                     aTTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGc                >  1:926592/1‑75 (MQ=255)
                                                      ttGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCa               >  2:665289/1‑75 (MQ=255)
                                                      ttGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCa               >  2:434581/1‑75 (MQ=255)
                                                       tGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAg              >  2:1123195/1‑75 (MQ=255)
                                                          cATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTg           >  1:1052921/1‑75 (MQ=255)
                                                            tCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCt         >  2:1060669/1‑75 (MQ=255)
                                                              aGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTcc       >  1:1133810/1‑75 (MQ=255)
                                                               gTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCa      >  1:1523320/1‑75 (MQ=255)
                                                                tAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAg     >  1:878604/1‑75 (MQ=255)
                                                                tAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAg     >  1:1198574/1‑75 (MQ=255)
                                                                 aaTGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAgg    >  1:775640/1‑75 (MQ=255)
                                                                  aTGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAggg   >  1:1246130/1‑75 (MQ=255)
                                                                  aTGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAggg   <  1:1027142/75‑1 (MQ=255)
                                                                   tGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAGGGt  >  1:1474935/1‑75 (MQ=255)
                                                                      |                                                                       
CAGAAGACGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAGGGT  >  REL606/3247228‑3247369

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