Predicted mutation
evidence seq id position mutation annotation gene description
RA 39 13,354 A→G intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*3913,3540AG96.7% 195.0 / ‑2.1 60intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (30/28);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

AGTGCCGGTAATACCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTATGGGCTCAGTT  >  39/13259‑13450
                                                                                               |                                                                                                
aGTGCCGGTAATACCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCa                                                                                             <  1:1258999/101‑1 (MQ=255)
 gTGCCGGTAATACCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAg                                                                                            >  2:1608259/1‑101 (MQ=255)
   tccGGTAATACCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCa                                                                                          <  1:707222/100‑1 (MQ=255)
      ggTAATACCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCaaaa                                                                                       >  2:99493/1‑101 (MQ=255)
       gTAATACCGGGGCACTGATGGGGCTGGCAAAATTACTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCGTCAGCAAAAg                                                                                      <  2:1090143/101‑1 (MQ=255)
           tACCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCa                                                                                  <  1:1617097/101‑1 (MQ=255)
           tACCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCa                                                                                  >  2:128386/1‑101 (MQ=255)
            aCCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCaa                                                                                 <  1:1545884/101‑1 (MQ=255)
              cGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCaaaa                                                                               <  1:1066339/101‑1 (MQ=255)
                   cACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACggtg                                                                          >  2:643324/1‑101 (MQ=255)
                    aCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACggtgg                                                                         >  1:1003944/1‑101 (MQ=255)
                    aCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACggtgg                                                                         <  1:844760/101‑1 (MQ=255)
                        aTGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCt                                                                     >  1:236700/1‑101 (MQ=255)
                           gggCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                  >  1:1006846/1‑101 (MQ=255)
                           gggCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                  <  2:524554/101‑1 (MQ=255)
                             gCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACt                                                                <  2:697481/101‑1 (MQ=255)
                              cTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACtt                                                               <  1:274640/101‑1 (MQ=255)
                              cTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACtt                                                               >  2:1201933/1‑101 (MQ=255)
                              cTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACtt                                                               >  2:453163/1‑101 (MQ=255)
                               tGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTa                                                              >  2:1574812/1‑101 (MQ=255)
                                 gCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAgg                                                            <  1:848304/101‑1 (MQ=255)
                                  cAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAggg                                                           <  1:1473687/101‑1 (MQ=255)
                                    aaaTTATTACTCAAGCCCGTGGAGGGGATTGTGCGTCCTGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTGCCTGACTTAGGGGc                                                         <  2:268032/101‑1 (MQ=255)
                                    aaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGc                                                         >  2:1247354/1‑101 (MQ=255)
                                    aaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGc                                                         >  1:1381111/1‑101 (MQ=255)
                                     aaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGcc                                                        >  1:720167/1‑101 (MQ=255)
                                        tattaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAAc                                                     <  1:1311015/101‑1 (MQ=255)
                                         attaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACg                                                    >  1:1542206/1‑101 (MQ=255)
                                           taCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTc                                                  <  1:183532/101‑1 (MQ=255)
                                            tcTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCg                                                 <  2:572490/100‑1 (MQ=255)
                                            aCTCAAGCCCCTGGAGGGGGTTGAGCGGCCGGCGCTGGGGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCg                                                 >  2:1689714/1‑101 (MQ=255)
                                            aCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCg                                                 <  1:655504/101‑1 (MQ=255)
                                            aCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCg                                                 <  2:1377905/101‑1 (MQ=255)
                                             cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGa                                                >  2:94009/1‑101 (MQ=255)
                                               cAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAtt                                              <  1:1608265/101‑1 (MQ=255)
                                                aaGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtg                                             >  1:402191/1‑101 (MQ=255)
                                                   ccccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGCCCTTGACTTCGGGGCCAACGTCGATTGTGa                                          >  2:1675568/1‑101 (MQ=255)
                                                    cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGAc                                         >  2:1270174/1‑101 (MQ=255)
                                                     ccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACa                                        <  2:364379/101‑1 (MQ=255)
                                                     ccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACa                                        >  1:1652763/1‑101 (MQ=255)
                                                     ccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACa                                        >  2:1504156/1‑101 (MQ=255)
                                                     ccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACa                                        <  1:1045454/101‑1 (MQ=255)
                                                       tGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGt                                      >  1:1013651/1‑101 (MQ=255)
                                                         gAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTAc                                    <  2:1652757/101‑1 (MQ=255)
                                                          aGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACa                                   <  2:1574783/101‑1 (MQ=255)
                                                          aGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACa                                   >  2:491321/1‑101 (MQ=255)
                                                             ggATTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATg                                >  1:818826/1‑101 (MQ=255)
                                                               aTTGAGCGTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGCGACAGAACAATGtt                              >  2:1142438/1‑101 (MQ=255)
                                                                 tGAGCGTCCGGCGCTGGTGACGGTATTACCGCATGAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTgg                            >  1:1701103/1‑101 (MQ=255)
                                                                      gTCCGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCaa                       >  1:1214713/1‑101 (MQ=255)
                                                                         cGGCGCTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAAttt                    <  2:1461312/101‑1 (MQ=255)
                                                                           gcgcTGGTGACGCTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGc                  <  1:608624/101‑1 (MQ=255)
                                                                              cTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCAt               >  1:1094325/1‑101 (MQ=255)
                                                                              cTGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCAt               <  1:453168/101‑1 (MQ=255)
                                                                               tGGTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCAtt              <  2:1330612/101‑1 (MQ=255)
                                                                                 gTGACGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTAt            <  2:1003939/101‑1 (MQ=255)
                                                                                 gTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTAt            >  2:1660518/1‑101 (MQ=255)
                                                                                    aCGGTATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTATggg         >  1:997594/1‑101 (MQ=255)
                                                                                        tATTACCGCATCATCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGGCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTATGGGCTCa     <  1:1403909/101‑1 (MQ=255)
                                                                                        tATTACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTATGGGCTCa     >  1:567664/1‑101 (MQ=255)
                                                                                           tACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTATGGGCTCAGtt  >  1:1547691/1‑101 (MQ=255)
                                                                                           tACCGCATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTATGGGCTCAGtt  >  2:1068913/1‑101 (MQ=255)
                                                                                               |                                                                                                
AGTGCCGGTAATACCGGGGCACTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGTACAATGTTGGTGCAATTTGCCATTATGGGCTCAGTT  >  39/13259‑13450

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