Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,318,527 +AC coding (1393/1446 nt) tldD ← predicted peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,318,5241.C85.7% 167.8 / 23.4 49coding (1396/1446 nt)tldDpredicted peptidase
Reads supporting (aligned to +/- strand):  ref base . (5/2);  new base C (26/16);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01
*REL6063,318,5242.A85.7% 157.8 / 21.7 49coding (1396/1446 nt)tldDpredicted peptidase
Reads supporting (aligned to +/- strand):  ref base . (5/2);  new base A (26/16);  total (31/18)
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.83e-01

TATCCTTCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGG  >  REL606/3318449‑3318588
                                                                              ||                                                                    
taTCCTTCTGAGAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGc‑‑c                                                                        >  1:53552/1‑76 (MQ=255)
taTCCTTCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGc‑‑c                                                                        >  1:284207/1‑76 (MQ=255)
     ttCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acac                                                               <  1:210273/76‑3 (MQ=255)
      tCTGACAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg                                                              >  2:1351015/1‑73 (MQ=255)
      tCTGAAAATATTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg                                                              >  1:1927666/1‑73 (MQ=255)
      tCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg                                                              >  2:403725/1‑73 (MQ=255)
      tCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg                                                              >  2:1351011/1‑73 (MQ=255)
      tCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg                                                              >  1:1310224/1‑73 (MQ=255)
               aCTNAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGc                                                     >  1:1838405/1‑76 (MQ=255)
                  gAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCaaa                                                  >  2:714562/1‑76 (MQ=255)
                   aaaTTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAAc                                                 <  1:704175/76‑1 (MQ=255)
                     attattACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACtt                                               >  1:1401622/1‑76 (MQ=255)
                       tgttaAGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAAACGGCAAACTTTg                                             <  2:85731/74‑1 (MQ=255)
                       tattaCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTg                                             <  2:1866457/76‑1 (MQ=255)
                       tattaCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTg                                             >  1:216467/1‑76 (MQ=255)
                       tagtaCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTg                                             <  2:1345736/73‑1 (MQ=255)
                          taCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGccc                                          >  2:1841521/1‑76 (MQ=255)
                           aCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCct                                         >  2:701884/1‑76 (MQ=255)
                            cGCAGTACCGCCAACAGTCAGGTGATGNACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCctt                                        >  2:1778/1‑76 (MQ=255)
                              cAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCcttct                                      >  1:411025/1‑76 (MQ=255)
                              cAGNACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCcttct                                      >  1:1834714/1‑76 (MQ=255)
                               aGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCcttctt                                     >  2:804133/1‑76 (MQ=255)
                               aGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCcttctt                                     >  2:1311553/1‑76 (MQ=255)
                                gTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCttt                                    <  1:1069199/76‑1 (MQ=255)
                                gTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCttt                                    <  2:1491489/76‑1 (MQ=255)
                                gTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑ACANGCCGACCGGCAAACTTTGCCCTTCttt                                    <  2:2619/76‑1 (MQ=255)
                                 tACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTa                                   >  2:673420/1‑76 (MQ=255)
                                   ccGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCGTTAcc                                 >  2:78680/1‑76 (MQ=255)
                                    cGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCg                                >  2:1684363/1‑76 (MQ=255)
                                     gCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGc                               <  2:1087375/76‑1 (MQ=255)
                                         aCAGTCAGGTTATCAACGGTCAGCATAGGCTGGCTAA‑‑G‑‑AC‑‑GCCGACCGCCAAACTTTGCCCTACTTTAGCGTAGac                           >  1:2041942/1‑76 (MQ=255)
                                           aGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGacac                         >  2:1849255/1‑76 (MQ=255)
                                              cAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCa                      >  1:1213665/1‑76 (MQ=255)
                                              cAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCa                      >  1:1477615/1‑76 (MQ=255)
                                                ggTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACg                    <  1:460068/76‑1 (MQ=255)
                                                ggTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACg                    <  2:530792/76‑1 (MQ=255)
                                                 gTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGc                   <  1:1311554/76‑1 (MQ=255)
                                                 gTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGc                   <  1:804134/76‑1 (MQ=255)
                                                  ttATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGcc                  >  2:1504264/1‑76 (MQ=255)
                                                    aTCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGt                >  1:659542/1‑76 (MQ=255)
                                                       aaCTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTAt             <  2:776341/76‑1 (MQ=255)
                                                       aaCTTTCAGCGTTGGCTCGC‑‑C‑‑CAGAG‑‑GCCAACCGGCAAANTTTGCACTTCTTTACNGCAGACCCCCACGCCGTTAt             >  2:2019605/1‑76 (MQ=255)
                                                        aCTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATc            >  2:801577/1‑76 (MQ=255)
                                                         cTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATcc           >  1:1054326/1‑76 (MQ=255)
                                                         cTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATcc           <  1:209357/76‑1 (MQ=255)
                                                         cTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATca           <  2:848952/76‑2 (MQ=255)
                                                          tttCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCa          <  2:876869/76‑1 (MQ=255)
                                                           ttCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAg         >  1:691041/1‑76 (MQ=255)
                                                            tCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGt        >  2:1172161/1‑76 (MQ=255)
                                                            tCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGt        >  1:580537/1‑76 (MQ=255)
                                                              aGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGttt      >  2:1834165/1‑76 (MQ=255)
                                                               gCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTc     <  1:1103189/76‑1 (MQ=255)
                                                                  ttGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACGCCCACGCCGNTATCCAGTTTCAgg  <  1:9528/76‑1 (MQ=255)
                                                                  ttGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAgg  >  1:1226596/1‑76 (MQ=255)
                                                                              ||                                                                    
TATCCTTCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGG  >  REL606/3318449‑3318588

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