Predicted mutation
evidence seq id position mutation annotation gene description
RA 40 30,807 2 bp→GA intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*4030,8070CG100.0% 115.6 / NA 37intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (20/17);  total (20/17)
*4030,8080CA97.3% 113.2 / NA 37intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (20/16);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCGTGTCATACAAGGTTTCGGCGGGGGCGGCCCCCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGTGAACGGCAAAAC  >  40/30712‑30908
                                                                                               ||                                                                                                    
gaataccccgaagataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGAccgc                                                                                                  >  1:1076218‑M2/91‑101 (MQ=255)
gaataccccgaagataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGAccgc                                                                                                  >  2:1590623‑M2/91‑101 (MQ=255)
 aataccccgaagataatcgggttggaaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGAccgcc                                                                                                 <  2:199144‑M2/12‑1 (MQ=255)
  ataccccgaagataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGAccgccg                                                                                                >  2:1590810‑M2/89‑101 (MQ=255)
  ataccccgaagataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgagatctgattaaagtgagcgccgcaacgGCGGCGAccgccg                                                                                                >  2:1709757‑M2/89‑101 (MQ=255)
           agataatcgggttgggaagcacgagccgcactattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATc                                                                                       <  1:1078725‑M2/22‑1 (MQ=255)
           agataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgTCGGCGACCGCCGTGGTTTATc                                                                                       >  1:465430‑M2/80‑101 (MQ=255)
                atcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCAt                                                                                  >  1:1689069‑M2/75‑101 (MQ=255)
                 tggggtgggaaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaccgGCGGCGACCGCCGTGGTTTATCCTCAtt                                                                                 <  1:607838‑M2/28‑1 (MQ=255)
                     gttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTAc                                                                             >  1:991082‑M2/70‑101 (MQ=255)
                         ggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGCCCGCCGTGGTTTATCCTCATTCTACGCtt                                                                         <  2:1689063‑M2/36‑1 (MQ=255)
                              cacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGc                                                                    >  2:1177343‑M2/61‑101 (MQ=255)
                              cacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGc                                                                    >  1:1088065‑M2/61‑101 (MQ=255)
                                    ccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCNTTCTANGCTTGCGGCAAGCGt                                                              >  2:1069908‑M2/55‑101 (MQ=255)
                                     cgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGtt                                                             >  2:1513357‑M2/54‑101 (MQ=255)
                                       cacgattttagtttgccccgtcgagatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTcc                                                           <  1:1709763‑M2/50‑1 (MQ=255)
                                        acgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCg                                                          <  1:1710387‑M2/51‑1 (MQ=255)
                                        acgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCg                                                          <  1:1106823‑M2/51‑1 (MQ=255)
                                             tttcgtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGCCGACCGCCGTGGTATATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGc                                                     >  2:462125‑M2/46‑101 (MQ=255)
                                             tttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGc                                                     >  2:1626794‑M2/46‑101 (MQ=255)
                                                 gtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTaca                                                 <  2:1459561‑M2/60‑1 (MQ=255)
                                                 gtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGCTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTaca                                                 >  2:747359‑M2/42‑101 (MQ=255)
                                                    tgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACAccc                                              >  1:918702‑M2/39‑101 (MQ=255)
                                                    tgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACAccc                                              <  2:245906‑M2/63‑1 (MQ=255)
                                                          gtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCa                                        >  1:946472‑M2/33‑101 (MQ=255)
                                                           tcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTCCCGGCAGCTACACCCGCGCCag                                       >  1:231755‑M2/32‑101 (MQ=255)
                                                             gcgatctgattaaagtgaccgccgcagcgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCagag                                     <  2:232834‑M2/72‑1 (MQ=255)
                                                               gatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGAt                                   >  1:1715335‑M2/28‑101 (MQ=255)
                                                               gatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGAt                                   >  1:1371656‑M2/28‑101 (MQ=255)
                                                                 tctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATaa                                 <  2:1353903‑M2/76‑1 (MQ=255)
                                                                 tctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATaa                                 <  2:1354578‑M2/76‑1 (MQ=255)
                                                                  ctgattaaagtgagcgccgcaacgTCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt                                <  2:465425‑M2/77‑1 (MQ=255)
                                                                  ctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt                                <  2:1076213‑M2/77‑1 (MQ=255)
                                                                  ctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt                                <  1:1349251‑M2/77‑1 (MQ=255)
                                                                  ctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt                                >  1:1296569‑M2/25‑101 (MQ=255)
                                                                  ctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt                                <  1:1702168‑M2/77‑1 (MQ=255)
                                                                         aagtgagngcngcaangGCGGNGACCGCCGTGGTTTATNCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCTCCAGAGATAATGCCCCTg                         <  2:1203060‑M2/84‑1 (MQ=255)
                                                                           gtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGac                       >  2:973106‑M2/16‑101 (MQ=255)
                                                                              agcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACt                    <  2:991077‑M2/89‑1 (MQ=255)
                                                                                  ccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCAGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAg                >  2:769146‑M2/9‑101 (MQ=255)
                                                                                    gcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGt              >  1:763240‑M2/7‑101 (MQ=255)
                                                                                         gGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTCCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGTGAACg         <  2:231755‑M2/100‑1 (MQ=255)
                                                                                                cccGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGTGAACGGCAAAAc  <  1:519009/101‑1 (MQ=255)
                                                                                               ||                                                                                                    
TCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCGTGTCATACAAGGTTTCGGCGGGGGCGGCCCCCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGTGAACGGCAAAAC  >  40/30712‑30908

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