Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 862,789 (CATAT)2→1 intergenic (‑13132/‑59425) unknown ← / → unknown tRNA‑Met‑CAT/tRNA‑Leu‑GAG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1862,7840C.86.7% 38.7 / 1.1 15intergenic (‑13127/‑59434)unknown/unknowntRNA‑Met‑CAT/tRNA‑Leu‑GAG
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (3/10);  total (3/12)
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.31e-01
*1862,7850A.86.7% 38.2 / 1.0 15intergenic (‑13128/‑59433)unknown/unknowntRNA‑Met‑CAT/tRNA‑Leu‑GAG
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (3/10);  total (3/12)
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.31e-01
*1862,7860T.86.7% 38.0 / 1.2 15intergenic (‑13129/‑59432)unknown/unknowntRNA‑Met‑CAT/tRNA‑Leu‑GAG
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (3/10);  total (3/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*1862,7870A.100.0% 53.6 / NA 15intergenic (‑13130/‑59431)unknown/unknowntRNA‑Met‑CAT/tRNA‑Leu‑GAG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (3/12);  total (3/12)
*1862,7880T.100.0% 53.6 / NA 15intergenic (‑13131/‑59430)unknown/unknowntRNA‑Met‑CAT/tRNA‑Leu‑GAG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (3/12);  total (3/12)

CTTTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAATTCGGAATGCGT‑AAACTATTTCGGAAAATCAGGGTCAAACGGGCATTCTTGCTTCAGTCGCCATTTATATTAAAGCCATATCATATATAGTTACACGTAATTCGGACTTTTAAGGTTGCTGCAGCGAGGCAATGATTCATCTATTTTCAACCT‑AAGGAATATCATGTTAAAAAAATACTTACTTGTTTTACCGATTTCTATTTTAACCGCCTGCGCGGGTTATTATCCTTATCCGCAATTTGAGCCCAGTTATCC‑GTCCTTTTCCGTGCAGCATAATGA  >  1/862671‑862986
                                                                                                                  |||||                                                                                                                                                                                                        
cttTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAAACGGGCTTTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATAATTAATCTATTTCCACCCTGAAGAAATATCATGTAAAAAAAATCCTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGtt                                                             <  1:189979/250‑1 (MQ=255)
      cGATCCGTATCTGAGCGAGAAGTAAAAAAAAGTCGAAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAACCGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCAT‑‑‑ATATGT‑‑TTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTTTTTTCAACCTGAAGGAATATCATGTAAAAAAAATACTTATTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatcc                                                      <  1:386685/251‑1 (MQ=255)
       gATCCGTATCTGAGCGAGAAGTAAGAAAAAGGCGGAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAAACGGGCTTTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTAAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatcc                                                      <  1:43973/250‑1 (MQ=255)
        aTCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCt                                                                                                                                             <  1:460939/162‑1 (MQ=255)
        aTCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatc                                                       >  1:444596/1‑248 (MQ=255)
        aTCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCtgaagg                                                                                                                             <  1:391927/178‑6 (MQ=255)
        aTCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTCCCGATTTCTATTTTAACCGCCTGTGCCGGtt                                                             <  1:580545/242‑1 (MQ=255)
            gTATCTGAGCGAGAAGTAAGAAAAAGTCGGAGAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCg                                                                                         >  1:42256/1‑210 (MQ=255)
             tATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCg                                               >  1:711569/1‑251 (MQ=255)
                                 aaaagtcggaaagtgaaaaaCTATTTCGGAAAACCAGGGTAAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCAT‑‑‑ATATGT‑‑TTATACGTAATTCAGGCTTTTAAGTTTGCTGCAGCGGGAAAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTTCTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATcc                           <  1:223172/234‑1 (MQ=255)
                                         gaaagtgtaaaaCTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAat                                                                                                                          <  1:424091/139‑1 (MQ=255)
                                             gtgtaaaaCTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAAc                                                                           <  1:771184/186‑1 (MQ=255)
                                             gtgtaaaaCTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCC‑ggtcc                     <  1:396556/239‑5 (MQ=255)
                                                         ttCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGTTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTCCCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa  <  1:710923/251‑1 (MQ=255)
                                                             gaaaaacggggTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGC‑‑‑‑‑CATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGAACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTTTTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa  <  1:143318/239‑1 (MQ=255)
                                                                                                                  |||||                                                                                                                                                                                                        
CTTTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAATTCGGAATGCGT‑AAACTATTTCGGAAAATCAGGGTCAAACGGGCATTCTTGCTTCAGTCGCCATTTATATTAAAGCCATATCATATATAGTTACACGTAATTCGGACTTTTAAGGTTGCTGCAGCGAGGCAATGATTCATCTATTTTCAACCT‑AAGGAATATCATGTTAAAAAAATACTTACTTGTTTTACCGATTTCTATTTTAACCGCCTGCGCGGGTTATTATCCTTATCCGCAATTTGAGCCCAGTTATCC‑GTCCTTTTCCGTGCAGCATAATGA  >  1/862671‑862986

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑