Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 5,217 C→T intergenic (–/+125917) – / ← unknown –/tRNA‑Ala‑GGC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*15,2170CT93.3% 65.4 / ‑2.9 30intergenic (–/+125917)–/unknown–/tRNA‑Ala‑GGC
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (18/10);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

CCGGTTTTAGCCCCACGCATTAAACACTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCACGCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAG  >  1/4975‑5438
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                             
ccGGTTTTAGCCCCACGCATTAAACACTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAAc                                                                                                                                                                                                                       <  1:493246/251‑1 (MQ=255)
               cGCATTAAACACTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAgg                                                                                                                                                                                                        >  1:735912/1‑251 (MQ=255)
                 cATTAAACACTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGAt                                                                                                                                                                                                      >  1:732022/1‑251 (MQ=255)
                  aTTAAACACTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGGCACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAg                                                                                                                                                                                                                     >  1:645318/1‑235 (MQ=255)
                        cacTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTa                                                                                                                                                                                               >  1:790635/1‑251 (MQ=255)
                          cTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTAt                                                                                                                                                                                              >  1:11409/1‑250 (MQ=255)
                          cTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTAt                                                                                                                                                                                              >  1:237662/1‑250 (MQ=255)
                          cTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTAt                                                                                                                                                                                              >  1:93824/1‑250 (MQ=255)
                                    gTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTa                                                                                                                                                                                               <  1:504748/239‑1 (MQ=255)
                                         cgGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGctgaa                                                                                                                                                                              >  1:635018/1‑247 (MQ=255)
                                         cgGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGctgaa                                                                                                                                                                              >  1:631828/1‑247 (MQ=255)
                                           gTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGctgaac                                                                                                                                                                             <  1:475827/250‑6 (MQ=255)
                                               cTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCCAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCAt                                                                                                                                                                        >  1:514412/1‑251 (MQ=255)
                                                 cccTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCAtca                                                                                                                                                                      >  1:555863/1‑249 (MQ=255)
                                                  ccTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCAtcag                                                                                                                                                                     <  1:134948/251‑4 (MQ=255)
                                                   cTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCTGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGt                                                                                                                                                                                     <  1:434363/234‑1 (MQ=255)
                                                   cTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGt                                                                                                                                                                                     >  1:592600/1‑234 (MQ=255)
                                                        gATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTg                                                                                                                                                                                         >  1:591187/1‑225 (MQ=255)
                                                         aTCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCAtcagatttct                                                                                                                                                               >  1:846594/1‑241 (MQ=21)
                                                          tCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCAtcagatttctcg                                                                                                                                                             <  1:443709/251‑12 (MQ=21)
                                                                   gtggCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCAtcagatttctcgttctcaag                                                                                                                                                     <  1:105708/250‑20 (MQ=18)
                                                                     ggCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAAc                                                                                                                                                                                                                       >  1:792451/1‑182 (MQ=255)
                                                                     ggCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCAtc                                                                                                                                                                       <  1:737475/230‑2 (MQ=255)
                                                                          tAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAgg                                                                                                                                                                                      >  1:156129/1‑210 (MQ=255)
                                                                                                        cggcggcggTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGctgaact                                                                                                                                                                            <  1:518583/190‑7 (MQ=255)
                                                                                                          gcggcggTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAAt                                                                                                            <  1:675978/250‑1 (MQ=11)
                                                                                                                                    gtAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAg                                                                                                                                                                                                                     >  1:524268/1‑121 (MQ=255)
                                                                                                                                          aaaTAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGt                                                                                                                                                                                                >  1:666057/1‑136 (MQ=255)
                                                                                                                                                                                                                     gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAg  >  1:897128/1‑251 (MQ=18)
                                                                                                                                                                                                                             ctTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAAtc                                                                                                                                                                  >  1:94868/1‑83 (MQ=25)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                             
CCGGTTTTAGCCCCACGCATTAAACACTTAGACGACGTGCGCGGTCACTCCCTCATGATCCACGCTGGTGGCGATAACCACTCCGATCACCCGGCCCCACTTGGCGGCGGCGGTCCGCGTATGGCGTGTGGTGTAATTAAATAATCCGCTTTCCGAAACAAAAAAGTGATCTGCACCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCACGCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAG  >  1/4975‑5438

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