Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	175	7,220	C→A	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	175	7,220	0	C	A	98.2%	189.0 / ‑5.7	56	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (31/24);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAA  >  175/7125‑7314                                                                                                |                                                                                               cttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGaa                                                                                           <  2:930206/101‑1 (MQ=255)  ttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAg                                                                                          >  1:945284/1‑101 (MQ=255)  ttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAg                                                                                          <  1:1184291/101‑1 (MQ=255)  ttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAg                                                                                          >  1:640029/1‑101 (MQ=255)  ttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAg                                                                                          >  2:863208/1‑101 (MQ=255)    atatGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCa                                                                                        <  1:714992/101‑1 (MQ=255)       tggtggGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAGCGTGGGAACGAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGaa                                                                                     <  2:1251254/101‑1 (MQ=255)       tggtggGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaa                                                                                     >  1:190202/1‑101 (MQ=255)       tggtggGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaa                                                                                     >  1:161853/1‑101 (MQ=255)              cGATAAATTAGAAGATATGCAGGCAGCGGCTGGGGCGAACGTGGGAACAAAAGTGCTGGTGCGCACGGGCAAACCTATTACGACTGAAGCAGAAAACgcgg                                                                              >  2:539293/1‑101 (MQ=255)              cGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcgg                                                                              <  1:1571382/101‑1 (MQ=255)               gATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggc                                                                             <  1:683231/101‑1 (MQ=255)               gATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggc                                                                             >  2:42311/1‑101 (MQ=255)                aTAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggcg                                                                            >  2:92683/1‑101 (MQ=255)                aTAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACgcggcg                                                                            >  2:1416139/1‑101 (MQ=255)                   aaTTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGaaaaaacggggggg                                                                         >  2:1071223/1‑91 (MQ=255)                   aaTTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAt                                                                         >  1:1674573/1‑101 (MQ=255)                         aagaTATGCAGGCAGCGGCTGCGGCGAACGGGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGgtg                                                                   >  2:210384/1‑101 (MQ=255)                            atatGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTa                                                                <  1:1302445/101‑1 (MQ=255)                              atGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa                                                              >  2:729895/1‑101 (MQ=255)                              atGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGGGTTaaa                                                              >  2:1626014/1‑101 (MQ=255)                               tGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAAt                                                             <  2:190202/101‑1 (MQ=255)                                gcaggcagCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATa                                                            <  1:1071228/101‑1 (MQ=255)                                 caggcagCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAg                                                           >  2:1698265/1‑101 (MQ=255)                                     cagCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTg                                                       <  2:73807/101‑1 (MQ=255)                                      agCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTgg                                                      >  1:1495763/1‑101 (MQ=255)                                       gCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGc                                                     >  1:351451/1‑101 (MQ=255)                                        cGGCTGCGGCCAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTCCGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCa                                                    <  1:1486687/101‑1 (MQ=255)                                           cTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGAc                                                 >  2:522857/1‑101 (MQ=255)                                            tGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGAcc                                                >  1:149035/1‑101 (MQ=255)                                                gcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGcc                                            <  2:946262/101‑1 (MQ=255)                                                  gAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGc                                          >  1:485559/1‑101 (MQ=255)                                                   aaCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCa                                         <  1:202809/101‑1 (MQ=255)                                                     cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAg                                       >  2:839870/1‑101 (MQ=255)                                                     cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAg                                       >  1:476412/1‑101 (MQ=255)                                                     cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAg                                       >  2:327231/1‑101 (MQ=255)                                                      gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGc                                      >  1:355035/1‑101 (MQ=255)                                                         ggACCAAAAGTGCTGGTGCGTAGGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGAt                                   <  2:1224629/101‑1 (MQ=255)                                                           aaCAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaa                                 >  1:1299724/1‑101 (MQ=255)                                                             cAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaaaa                               <  1:282054/101‑1 (MQ=255)                                                              aaaaGTGCTGGTGCGTACTGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaaaaa                              >  1:1258976/1‑101 (MQ=255)                                                                   tGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAagcag                         <  1:1416145/101‑1 (MQ=255)                                                                     cgggTGCCTACGGCCAACCCTATTCCGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATCGCCTGGCAGCCCTGCCGCAAGCGATAAAAAAgcagca                       <  2:765959/99‑1 (MQ=255)                                                                     cTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAgcagca                       <  2:671027/101‑1 (MQ=255)                                                                      tGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCaa                      <  2:215879/101‑1 (MQ=255)                                                                        gTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCaaaa                    <  1:488503/101‑1 (MQ=255)                                                                           cGTACGGGCAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAAcc                 <  2:1458285/101‑1 (MQ=255)                                                                            gTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCg                <  2:485554/101‑1 (MQ=255)                                                                               cGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGCCCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCa             <  1:1389230/101‑1 (MQ=255)                                                                               cGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCa             >  2:106385/1‑101 (MQ=255)                                                                                gggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCAt            <  2:266164/101‑1 (MQ=255)                                                                                 ggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATa           >  1:999203/1‑101 (MQ=255)                                                                                  gTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATaa          >  2:591434/1‑101 (MQ=255)                                                                                  gTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATaa          >  1:614471/1‑101 (MQ=255)                                                                                    aaaCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATg        <  2:716645/101‑1 (MQ=255)                                                                                    aaaCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATg        >  1:876732/1‑101 (MQ=255)                                                                                    aaaCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATg        <  1:839875/101‑1 (MQ=255)                                                                                      aCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGAt      <  2:476407/101‑1 (MQ=255)                                                                                        cTATAACGACTGAAGCAGAAACCGCGGCGGTTGGGGTGTTAAATAGCCTGCCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTa    <  1:42311/101‑1 (MQ=255)                                                                                         tATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTaa   <  2:479007/101‑1 (MQ=255)                                                                                          aTTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTaaa  >  1:246633/1‑101 (MQ=255)                                                                                                |                                                                                               CTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAA  >  175/7125‑7314

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑