Predicted mutation
evidence seq id position mutation annotation gene description
RA 44 8,024 A→G intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*448,0240AG98.4% 215.6 / ‑4.6 64intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (32/31);  total (33/31)
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.07e-01

ACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACG  >  44/7928‑8121
                                                                                                |                                                                                                 
acaATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                                                                               >  2:1246304/1‑101 (MQ=255)
acaATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                                                                               <  2:1359626/101‑1 (MQ=255)
       cGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCa                                                                                        <  2:478614/101‑1 (MQ=255)
         gTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACg                                                                                      >  1:994120/1‑101 (MQ=255)
         gTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACg                                                                                      >  1:1483251/1‑101 (MQ=255)
          tCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACgg                                                                                     <  1:351014/101‑1 (MQ=255)
          tCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACgg                                                                                     <  2:445531/101‑1 (MQ=255)
            aGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGAc                                                                                   >  2:1526294/1‑101 (MQ=255)
            aGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGAc                                                                                   >  1:835691/1‑101 (MQ=255)
               ggTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTc                                                                                >  1:1023560/1‑101 (MQ=255)
                 ttCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCAt                                                                              >  1:220264/1‑101 (MQ=255)
                   ccTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCTAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                                                                            >  2:1371754/1‑101 (MQ=255)
                   ccTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                                                                            <  1:670814/101‑1 (MQ=255)
                   acTGAAGTGCTTCGAGGTATTCATTGGAGTACAGGGTAAATCCGCGTGGAACTCCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATcn                                                                            <  1:825709/100‑2 (MQ=255)
                          tGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACCCCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGAt                                                                     <  2:994115/101‑1 (MQ=255)
                            cTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                                                                   <  1:1248793/101‑1 (MQ=255)
                            cTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                                                                   >  2:1023568/1‑101 (MQ=255)
                               cGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGcc                                                                >  1:678367/1‑101 (MQ=255)
                               cGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGcc                                                                >  2:381305/1‑101 (MQ=255)
                                gAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCg                                                               >  1:1672540/1‑101 (MQ=255)
                                gAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCg                                                               >  2:308692/1‑101 (MQ=255)
                                 aGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGa                                                              >  2:1404642/1‑101 (MQ=255)
                                    aTTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAAcc                                                           >  1:1065987/1‑101 (MQ=255)
                                    aTTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAAcc                                                           >  1:1330002/1‑101 (MQ=255)
                                       tCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGCTAAGCCGGACCGAt                                                        >  2:668250/1‑101 (MQ=255)
                                       tCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAt                                                        <  2:1422787/101‑1 (MQ=255)
                                           tGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCg                                                    >  1:254444/1‑101 (MQ=255)
                                           tGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCg                                                    >  1:1103290/1‑101 (MQ=255)
                                              aGTTCATGGTAAAACCGCGTGGACCCCCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGc                                                 <  2:724706/101‑1 (MQ=255)
                                                   aGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGAt                                            >  2:1306056/1‑101 (MQ=255)
                                                    gggTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATa                                           <  2:1449108/101‑1 (MQ=255)
                                                     ggTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATAc                                          <  1:515230/101‑1 (MQ=255)
                                                        aaaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGtt                                       <  1:1023573/101‑1 (MQ=255)
                                                        aaaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGtt                                       >  2:948663/1‑101 (MQ=255)
                                                         aaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTc                                      <  2:1415426/101‑1 (MQ=255)
                                                         aaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTc                                      >  1:1719007/1‑101 (MQ=255)
                                                          aaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCt                                     <  1:1246309/101‑1 (MQ=255)
                                                          aaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCt                                     <  1:400296/101‑1 (MQ=255)
                                                          aaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCt                                     <  2:1346287/101‑1 (MQ=255)
                                                           aCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTg                                    <  2:716349/101‑1 (MQ=255)
                                                            ccGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGc                                   <  2:1051526/101‑1 (MQ=255)
                                                             cgcgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCg                                  >  2:1468708/1‑101 (MQ=255)
                                                                      cacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCa                         <  2:1071645/101‑1 (MQ=255)
                                                                       acCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCAc                        >  1:912108/1‑101 (MQ=255)
                                                                       acCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCAc                        <  2:39098/101‑1 (MQ=255)
                                                                       acCTTCTGCCAGATTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCTGAACCGATTTCGTGCTTGCTACGTTGTGCgcagacttgcc                        >  2:840158/1‑91 (MQ=255)
                                                                         cTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCg                      <  1:1498544/101‑1 (MQ=255)
                                                                         cTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCg                      <  1:866759/101‑1 (MQ=255)
                                                                         cTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCg                      >  1:911584/1‑101 (MQ=255)
                                                                             tGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCa                  <  1:501225/101‑1 (MQ=255)
                                                                               ccAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCaga                >  2:819205/1‑101 (MQ=255)
                                                                                 aGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAgaga              >  1:1014637/1‑101 (MQ=255)
                                                                                  ggTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTTCGGTGGCTTCACCGATCAGAGaa             >  1:255637/1‑101 (MQ=255)
                                                                                  ggTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGaa             >  1:1166874/1‑101 (MQ=255)
                                                                                  ggTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGaa             >  2:983368/1‑101 (MQ=255)
                                                                                    ttACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAAcc           <  1:899440/101‑1 (MQ=255)
                                                                                    ttACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAAcc           <  1:181802/101‑1 (MQ=255)
                                                                                    ttACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAAcc           <  1:624137/101‑1 (MQ=255)
                                                                                       cGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTa        <  2:371248/101‑1 (MQ=255)
                                                                                         gCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAAt      <  2:1329996/101‑1 (MQ=255)
                                                                                         gCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAAt      <  1:1428403/101‑1 (MQ=255)
                                                                                          ccACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAAtt     >  1:1691241/1‑101 (MQ=255)
                                                                                          ccACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAAtt     <  2:592906/101‑1 (MQ=255)
                                                                                           cACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTa    >  2:1436645/1‑101 (MQ=255)
                                                                                             cGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTacg  >  1:979103/1‑101 (MQ=255)
                                                                                                |                                                                                                 
ACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACG  >  44/7928‑8121

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: