Predicted mutation
evidence seq id position mutation annotation gene description
RA 10 21,791 T→C intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1021,7910TC98.0% 167.5 / ‑4.2 49intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (13/35);  total (14/35)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

GTGACTGATGGTGGAAATAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCACTGCTCATAG  >  10/21693‑21888
                                                                                                  |                                                                                                 
gTGACTGATGGTGGAAATAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCccacca                                                                                                 >  1:850603/1‑101 (MQ=255)
    cTGATGGTGGAAATAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAAcc                                                                                             <  2:532315/101‑1 (MQ=255)
      gATGGTGGAAATAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGt                                                                                           <  1:87840/101‑1 (MQ=255)
         ggtggAAATAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGAc                                                                                        <  2:1509297/101‑1 (MQ=255)
           tggAAATAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGt                                                                                      >  1:129267/1‑101 (MQ=255)
              aaaTAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTa                                                                                   <  1:236394/101‑1 (MQ=255)
              aaaTAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTa                                                                                   <  1:838604/101‑1 (MQ=255)
              aaaTAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTa                                                                                   >  2:1158702/1‑101 (MQ=255)
              aaaTAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTa                                                                                   <  2:1402817/101‑1 (MQ=255)
               aaTAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTAc                                                                                  <  1:184187/101‑1 (MQ=255)
                  aaaCGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTAcccc                                                                               >  1:1493017/1‑100 (MQ=255)
                    aCGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCAcg                                                                             <  1:332570/101‑1 (MQ=255)
                              gTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCgtagt                                                                   <  1:122463/101‑1 (MQ=255)
                                cgcgGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTAc                                                                 >  1:1156661/1‑101 (MQ=255)
                                 gcgGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACg                                                                <  2:284383/101‑1 (MQ=255)
                                 gcgGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACg                                                                <  2:1455029/101‑1 (MQ=255)
                                  cgGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACgg                                                               <  1:168958/101‑1 (MQ=255)
                                   ggCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGc                                                              >  1:1309259/1‑101 (MQ=255)
                                    gCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCa                                                             <  1:905628/101‑1 (MQ=255)
                                        gcgcGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACttt                                                         <  1:583356/101‑1 (MQ=255)
                                          gcgTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTAc                                                       <  1:1382550/101‑1 (MQ=255)
                                            gTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACaa                                                     <  1:117696/101‑1 (MQ=255)
                                                   cgcATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCaa                                              <  1:1084604/101‑1 (MQ=255)
                                                   cgcATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCaa                                              <  1:1591900/101‑1 (MQ=255)
                                                    gcATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCaaa                                             <  1:1710726/101‑1 (MQ=255)
                                                     cATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAg                                            <  2:558680/101‑1 (MQ=255)
                                                     cATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAg                                            <  1:1186263/101‑1 (MQ=255)
                                                       tCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCa                                          <  1:893841/101‑1 (MQ=255)
                                                        ccTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAg                                         <  1:506817/101‑1 (MQ=255)
                                                         cTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAgg                                        <  1:1007416/101‑1 (MQ=255)
                                                         cTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAgg                                        <  2:1381688/101‑1 (MQ=255)
                                                         cTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAgg                                        <  1:64332/101‑1 (MQ=255)
                                                           gATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCt                                      <  1:1130134/101‑1 (MQ=255)
                                                            aTACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCtt                                     >  2:1230350/1‑101 (MQ=255)
                                                             tACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCttt                                    <  1:581658/101‑1 (MQ=255)
                                                             tACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCttt                                    <  1:643769/101‑1 (MQ=255)
                                                              aCCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCCCCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTg                                   <  2:108818/101‑1 (MQ=255)
                                                               attcTTACGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGc                                  >  2:1212988/4‑101 (MQ=255)
                                                                  cTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGcgtc                               >  1:755267/1‑101 (MQ=255)
                                                                     aCGCTCGGTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGcgtccgt                            <  2:814945/101‑1 (MQ=255)
                                                                           ggTAATGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGAt                      >  1:396408/1‑101 (MQ=255)
                                                                              aaTGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAAt                   >  1:1450607/1‑101 (MQ=255)
                                                                                tGTCGCGACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATcc                 >  1:383500/1‑101 (MQ=255)
                                                                                   cgcgACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGc              <  1:1222288/101‑1 (MQ=255)
                                                                                      gACCAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCAc           >  2:1408381/1‑101 (MQ=255)
                                                                                        ccAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCACTg         <  1:545859/101‑1 (MQ=255)
                                                                                         cAAAGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCACTGc        >  1:1003557/1‑101 (MQ=255)
                                                                                          aaaGCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCACTGCt       >  2:1004405/1‑101 (MQ=255)
                                                                                             gCCCACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCACTGCTCAt    <  2:1079793/101‑1 (MQ=255)
                                                                                              cccACCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCACTGCTCATa   <  1:616637/101‑1 (MQ=255)
                                                                                               ccaccaAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCACTGCTCATAg  <  1:1007218/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
GTGACTGATGGTGGAAATAAACGTCGTTTTGTCGCGGCTGGCGCGTTCAAGCGCATCCTGATACCGCTTACGCTCGGTAATGTCGCGACCAAAGCCCATCAAACCGTGACGTTTACCCACGCGGTCGTAGTACGGCACTTTACAAATTTCAAAGCAGGCTTTGCGTCCGTCGGGATAATCCAGCCACTGCTCATAG  >  10/21693‑21888

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