Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 430,954 G→A E4E (GAG→GAA lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606430,9540GA98.0% 185.8 / ‑3.6 52E4E (GAG→GAAlonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (29/21);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

ATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCC  >  REL606/430880‑431026
                                                                          |                                                                        
aTGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGc                                                                         <  1:445834/76‑1 (MQ=255)
    cGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTc                                                                     >  2:1406521/1‑76 (MQ=255)
      tGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTg                                                                   >  2:410983/1‑76 (MQ=255)
       gAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGa                                                                  >  2:705206/1‑76 (MQ=255)
       gAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATACTGAACGTTCTGa                                                                  >  1:1315535/1‑76 (MQ=255)
        aaGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGaa                                                                 >  1:483254/1‑76 (MQ=255)
         aGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAAc                                                                <  1:1502626/76‑1 (MQ=255)
           cacaGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGc                                                              >  2:495758/1‑76 (MQ=255)
               gTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTg                                                          <  1:992715/76‑1 (MQ=255)
                  gtgtCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGaaa                                                       >  1:74465/1‑76 (MQ=255)
                           gATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTAt                                              <  1:1334481/76‑1 (MQ=255)
                           gATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTAt                                              >  1:977268/1‑76 (MQ=255)
                            aTTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTAtt                                             >  1:846938/1‑76 (MQ=255)
                                ccTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCg                                         >  2:566933/1‑76 (MQ=255)
                                ccTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCg                                         <  1:604910/76‑1 (MQ=255)
                                 cTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGc                                        >  1:1397576/1‑76 (MQ=255)
                                  tGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCt                                       <  2:921391/76‑1 (MQ=255)
                                     cGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTgcg                                    >  2:82678/1‑76 (MQ=255)
                                     cGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTgcg                                    <  1:706428/76‑1 (MQ=255)
                                      ggAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTgcgc                                   <  1:498339/76‑1 (MQ=255)
                                      ggAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTgcgc                                   <  1:470445/76‑1 (MQ=255)
                                        aaaTTAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGa                                 <  2:1359808/76‑1 (MQ=255)
                                           ttAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGAtgt                              <  1:849422/76‑1 (MQ=255)
                                           ttAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGAtgt                              >  1:592789/1‑76 (MQ=255)
                                            tAAACTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGAtgtg                             >  2:725324/1‑76 (MQ=255)
                                            tAAACTAAGAGAGAGCTCNATGAATCCTGAACGTTCTGATCGCATTAAAATCCCCGTATTGCCGCTGCGCGAtgtg                             >  1:1466302/1‑76 (MQ=255)
                                             aaaCTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATgtgg                            >  2:603933/1‑76 (MQ=255)
                                             aaaCTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATgtgg                            >  2:831740/1‑76 (MQ=255)
                                             aaaCTAAGACAGAGCTCTATGAANCCTGANCGTTCTGAACGCATNGAGATCACCGTATTGCTGCTGCGCGATgtgg                            >  2:206534/1‑76 (MQ=255)
                                               aCTAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATgtggtg                          >  2:1407023/1‑76 (MQ=255)
                                                 tAAGAGAGAGCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTggtggt                        >  2:937908/1‑76 (MQ=255)
                                                 tAAGAGAGAGCTCTATGAATCCTGAACGTTCCGAACGCATNNAAATCCCCGTNNTACCGCTGCGCGATGTggtggt                        >  1:1372229/1‑76 (MQ=255)
                                                   agagagagCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGttt                      <  1:763247/76‑1 (MQ=255)
                                                   agagagagCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGACGTGGTGGttt                      >  2:951938/1‑76 (MQ=255)
                                                     agagagCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTAt                    <  1:748497/76‑1 (MQ=255)
                                                        gagCGCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCg                 >  1:86266/1‑76 (MQ=255)
                                                          gCTCTATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGca               >  1:545512/1‑76 (MQ=255)
                                                           ctctATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGcac              >  2:1154465/1‑76 (MQ=255)
                                                            tctATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGAGATGTGGTGGTTTATCCGcaca             <  1:9110/76‑1 (MQ=255)
                                                              tATGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATg           >  1:1100314/1‑76 (MQ=255)
                                                               aTGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATgg          <  2:953297/76‑1 (MQ=255)
                                                                tGAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGt         <  2:168336/76‑1 (MQ=255)
                                                                 gAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTc        <  1:495451/76‑1 (MQ=255)
                                                                 gAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTc        <  2:910365/76‑1 (MQ=255)
                                                                 gAATCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTc        >  1:698765/1‑76 (MQ=255)
                                                                 gAATCCTGAACGTNCTGAACGCATTGAAATGCCCGTATTGCCGCTTCGCGATGTGGTGGTTTATCCGCACATGGTc        <  2:1397575/76‑1 (MQ=255)
                                                                   aTCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCAt      <  1:494603/76‑1 (MQ=255)
                                                                   aTCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTAAt      >  2:611872/1‑76 (MQ=255)
                                                                    tCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATc     <  1:659956/76‑1 (MQ=255)
                                                                    tCCTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATc     >  1:964649/1‑76 (MQ=255)
                                                                      cTGATCCTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATccc   >  1:84223/1‑76 (MQ=255)
                                                                      cTGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGNGGTGGTTTATCCGCACATGGTCATccc   <  1:1525796/76‑1 (MQ=255)
                                                                       tGAACGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATcccc  >  1:598281/1‑76 (MQ=255)
                                                                          |                                                                        
ATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCC  >  REL606/430880‑431026

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