Predicted mutation
evidence seq id position mutation annotation gene description
RA 178 928 G→A intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1789280GA95.8% 73.0 / ‑3.9 24intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (14/9);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GGTAAGATTTAACATATCCGGGGATTTGAAGCCGATAAATCCTGATAAATATCCATGAACGTAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGA  >  178/898‑1025
                              |                                                                                                 
gggacagataaaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACt                             >  2:674865/11‑101 (MQ=255)
 ggacagctaaaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCGATGAACGCAAAAACAGGATACGACCTGCCGGCTGCGgtactggtactg                            >  2:1631107/10‑101 (MQ=255)
 ggacagataaaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTg                            <  2:659166/92‑1 (MQ=255)
   acagataaaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGAt                          >  1:1197181/8‑101 (MQ=255)
    cagataaaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGAtt                         <  2:1396358/95‑1 (MQ=255)
     agataaaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTg                        >  2:154656/6‑101 (MQ=255)
     agataaaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTg                        <  1:1027431/96‑1 (MQ=255)
     agataaaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTg                        >  2:4284/6‑101 (MQ=255)
       tttAACATATCCGGGGATTTGAAGCCGATAAATCCTGATAAATATCCATGAACGTAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTgcc                      >  1:1586216/1‑101 (MQ=255)
          aaCATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTgccgtc                   <  2:672952/101‑3 (MQ=255)
            cATATCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTgccgtcgg                 >  2:1411185/1‑97 (MQ=255)
             atatCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTgccgtcggt                <  1:1558566/101‑6 (MQ=255)
             atatCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGGCGGCTGCCGTTCTGGCACTGATTgccgtcggt                <  2:1163612/101‑6 (MQ=255)
              tatCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTg               <  2:1313376/101‑1 (MQ=255)
              tatCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTg               <  2:1074918/101‑1 (MQ=255)
                  cagggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGcc           >  1:1122916/10‑101 (MQ=255)
                   agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt          <  2:973764/93‑1 (MQ=255)
                   agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt          >  2:124465/9‑101 (MQ=255)
                   agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt          >  1:630300/9‑101 (MQ=255)
                   agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt          >  1:754578/9‑101 (MQ=255)
                     ggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGc        >  2:241240/7‑101 (MQ=255)
                         cagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTgcgcctgcgcct    >  2:152024/3‑101 (MQ=255)
                         aagAACCCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTgcgcctgcgcct    <  2:322058/99‑1 (MQ=255)
                          agAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGcgcctgcgcctg   >  2:773086/2‑101 (MQ=255)
                          agAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGcgcctgcgcctg   >  2:1449259/2‑101 (MQ=255)
                           gAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGa  <  2:402908/101‑1 (MQ=255)
                              |                                                                                                 
GGTAAGATTTAACATATCCGGGGATTTGAAGCCGATAAATCCTGATAAATATCCATGAACGTAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGA  >  178/898‑1025

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