Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,586,710	Δ1 bp	9.3%	coding (114/1743 nt)	ggt ←	gamma‑glutamyltranspeptidase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,586,710	0	T→.	9.3%	6.5 0.0 ‑0.4	43	coding (114/1743 nt)	ggt	gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (20/18):  total (22/21)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.01e-01

GATTTTTTGCTGTCCGGGTTGCCCTGATCATCGAGGAACATATCGCGGGTCGCTTTGGCGGGTGCCATTTCGCGGAAATCGATAGCCGTGGTATTGCCATTTTTCGAGCGGATTAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAGAGGAGGATACGCGATGAAACGACTTCTGATTCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAAT  >  NC_000913/3586411‑3587005                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                        gaTTTTTTGCTGTCCGGGTTGCCCTGATCATCGAGGAACATATCGCGGGTCGCTTTGGCGGGTGCCATTTCGCGGAAATCGATAGCCGTGGTATTGCCATTTTTCGAGCGGATTAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTCCCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCtt                                                                                                                                                                                                                                                                                                        >  1:149934/1‑301 (MQ=255)     ttttGGTGGCCGGGTTGCCCTGATCAGCGAGGAACATATCGCGGGTCGCTTTGGCGGGTGCCATTTCGCGGAAATCGATAGCCGTGGTATTGCCATTTTTCGAGCGGATTAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTtcctc                                                                                                                                                                                                                                                                                                    <  2:149927/301‑1 (MQ=255)                    tGCCCTGATCATCGAGGAACATATCGCGGGTCGCTTTGGCGGGTGCCATTTCGCGGAAATCGATAGCCGTGGTATTGCCATTTTTCGAGCGGATTAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGAc                                                                                                                                                                                                                                                                                     <  1:506393/301‑1 (MQ=255)                          gATCATCGAGGAACATATCGCGGGTCGCTTTGGCGGGTGCCATTTCGCGGAAATCGATAGCCGTGGTATTGCCATTTTTCGAGCGGATTAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCACCCTGAGTGCCAGTGGCGTCCACAGACGCTACCCTTCCCTGTTCCCCGCGGACCGGTGGGCAGACACCTTCCCCCCCAcccacccacccgggccg                                                                                                                                                                                                                                                                               >  2:293073/1‑286 (MQ=255)                                                    gCTTTGGCGGGTGCCATTTCGCGGAAATCGATAGCCGTGGTATTGCCATTTTTCGAGCGGATTAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGTCGTCCACCGACGCTACCATTCCCTGTTTCGCGCGTACCTGGTGGACGACATCTTCCTCCACCCCATACGAGACGGGCGGCGCCGAGGGCGCGGCGGCGgcgccca                                                                                                                                                                                                                                                     >  2:466400/1‑298 (MQ=255)                                                                     ttgccggtaacacagtgcctttgTATTGCCATTTTTCGGGCGGATTAACATACAACCACCCCCGTCCCGTGCCCCTTCCTGCGGATGCGTTCCCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCACTCCTTGTGAATATCCACCCCCACCTGCGTGGCAGTGGCGTCCACAGACGCTACCATTCCCGGTTGCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGTGCGGCTCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTgag                                                                                                                                                                                                                                       <  1:720955/276‑1 (MQ=255)                                                                           gAAATCGATAGCCGTGGTATTGCCATTTTTCGAGCGGATCAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACAGCCAGCGCGTAGCTCACCGCCACGGCGGCATCAACGGCATTCACGCCCTCCTTGAGCATACCCACCCCCACCTGAGTGGCAGTGCCGTCCACCGACACTACCATTCCCTGTTTCGCTCCTACCGGGTGGGCGACATCTTCCTCGACACCATAAGAGACGGGAGGGCCAGGAGGCTCGGCGGCGGCGCGACACCAACTTccacaccacagagc                                                                                                                                                                                                                               >  2:364423/1‑288 (MQ=255)                                                                                cGATAGCCGTGGTATTGCCATTTTTCGAGCGGATTAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGCCGCTACCATTCCCTGTTTCGCGCGTACCGGTTTGCAGCCATCTTACTCCACACCATACCCGACCGGCGGCGCAGGAGGCCCGGCGGCGGCGCGAAcaccgtaccgagagaagcgcacccacgg                                                                                                                                                                                                                         >  2:517364/1‑276 (MQ=255)                                                                                                                          aaCCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCCGGAGGCGCGGCGGCGGCGCTAACACAACTTCCTGAGAGCAGCGCCGCACTGGCCCCCCGGCGTAAAAACGTCCGTTTTATCATACTTAtt                                                                                                                                                                                 >  1:596498/1‑299 (MQ=255)                                                                                                                          aaCCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGCGAGCAGCGCAGCAAGGGCACGCCGGCGTAAGCACGTCGGTTTTATCATCGTTAtt                                                                                                                                                                                 >  1:504581/1‑299 (MQ=255)                                                                                                                                 cctgCCCAGATTCCCTGCCGGCGGTGGCGTTACCCCCAGCGCGGAGCCCACCGCCACGTCTGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTGCCTCCAAACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAgaga                                                                                                                                                                        <  1:293079/298‑1 (MQ=255)                                                                                                                                     cccAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCCTTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCCGCGGCGGCGCCACAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTCAAAACGTCGGCTTTATCATCGCTCTTCCCCACACATTa                                                                                                                                                                     >  2:352199/1‑300 (MQ=255)                                                                                                                                           ttCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCCGAGCAGCCATGGCCCCCCGGCGTAAAAACGACGGTTTTCGCCTCGTCTTTCCCCCGAGCTTAAGGGGCa                                                                                                                                                              >  2:155868/1‑301 (MQ=255)                                                                                                                                                     gcGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCAGCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCt                                                                                                                                                     <  1:466414/300‑1 (MQ=255)                                                                                                                                                          cgcgtccccGCCCGCGCGTAGCCCACCGCCCCGGCGGCCGCAACGGCATTCCCGTCCTCCTTGAGAATATCCACCCCCACCGGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCt                                                                                                                                               <  2:504558/294‑1 (MQ=255)                                                                                                                                                             gTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCCTACGACACGTGCGGCGCAGGACGCGCGACGCAGGAGCTAAAACAACTTACTTCGAGCGGAGCCGCACATGCCCCCCGGAGTCACAACAGCGGTTTTATCCACGTTATTCTCCAAGGATTAAGGTGCAACACCAGCTAACCCTggg                                                                                                                                            >  2:589277/1‑300 (MQ=255)                                                                                                                                                               tACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGGGAATATCCACCCCCACCTGAGTGGCAGGGGCGTCCACAGACGCTACCATTCCCTGGTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGtat                                                                                                                                          <  2:596473/301‑1 (MQ=255)                                                                                                                                                                    cccgcgcGTACCCCAACGCCACGGCCTCATCAACGGTATTCCCGCCCGCCTAGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGGACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCT‑‑AACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAActc                                                                                                                                   <  1:517387/298‑1 (MQ=255)                                                                                                                                                                    ccAGCGCGTCGCCCACCGCCACGGCGGCATACACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCGGAGTTGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGGACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAAc                                                                                                                                     <  1:155875/301‑1 (MQ=255)                                                                                                                                                                              gCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATC‑TCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAACACACCTTCCTGCGAGCAGAGCAGCCATGGCCCCCCGGCGTACACACGTCGGTTTTATCCTCGTTATTATCCAGAGCTTCAGGGGCACCCCCCGCTAAGCCTGGCATCTAACTCTGACTTaaa                                                                                                                          >  2:337638/1‑300 (MQ=255)                                                                                                                                                                               cccACCGCCCCGGCGGCAGCAACGGCATTCCCGCCCTCCTGGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCAGTCCCGGTTTCGCGCGTACCGGGTGGAAGACATC‑TCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAatc                                                                                                                         <  1:337644/301‑1 (MQ=255)                                                                                                                                                                                ccACCGCCACGGCGGCCTCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCTCCTGAGTGGCAGGTGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGGACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAatc                                                                                                                         <  1:491059/301‑1 (MQ=255)                                                                                                                                                                                                                 cctcctTGAGACAATCCGCCCCCCCCTCGGTGTCATGTGCGTACCCAGACCCTACCATTCCCTGGTCCGCGCGTACCGTGTGGAAGACTGCTTCCTCCACACCATACGAGACGGGCTGCGCAGGGGGCTCGGCGGCGGCTCTAAAACAACTACCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGAATAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAgg                                                                                        <  1:364432/301‑1 (MQ=255)                                                                                                                                                                                                                 ccgcctTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTGTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAgg                                                                                        <  1:472858/298‑1 (MQ=255)                                                                                                                                                                                                                   tcctTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATC‑TCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCGATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCACAGCTAAGCCTGGTATATCACCCTGAATTAATCATCGTTGTGTCGGTAAGCGAGTAAACTAAAAGGaga                                                                                     >  2:361459/1‑299 (MQ=255)                                                                                                                                                                                                                    cctTGAGAATATCCACCCCCACCTGCGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGGAGACATCTTCCTCCACACCATACGAGACGGTCGGCGCAGGAGGCGCGGCGGCGGCGCTCAACCAACTTCCTGAGAGCAGCGCCTCAATGGCCACCCGGCGTAAAAACGCCGGTTTTCTCATCGTTATTCTCCAGCGATTACGGGGCAACCCCAGCTAAGCCTGGTATATCACTTTGACTTAACCTTCGGTTTGTCGGTACTCGATTAAACTCAAAGGata                                                                                     >  2:344036/1‑301 (MQ=255)                                                                                                                                                                                                                    cctTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATAAAACTCTGAATTAATCATCGTTTTGTCGGGCAGCGAGTAAACTTAACGGata                                                                                     >  2:282073/1‑301 (MQ=255)                                                                                                                                                                                                                         cgaagagccccccctcccctgggggcagtggcGTCCACAGACGCTACCATTCCCGGTTTCGCGCGTACCGGGTGGAAGACATCGTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATAAAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTc                                                                                <  1:282079/278‑1 (MQ=255)                                                                                                                                                                                                                         agaCTATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTatca                                                                                                                                                                                         <  1:529866/196‑1 (MQ=255)                                                                                                                                                                                                                         agaATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTatca                                                                                                                                                                                         >  2:529843/1‑196 (MQ=255)                                                                                                                                                                                                                         agaATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGGAAAAACGTCGGTTTTATCATCGTTATTCTCCAAAGATTAAGGGACAACCCCAGCTAAGCGTGGTATATAACTCTGAATTAATGCTCGATTTTTCGGGACGCGAGTCAACTTAAAGGCTATCtt                                                                                 >  2:656381/1‑300 (MQ=255)                                                                                                                                                                                                                              cgccACCCCCACCTGCGTGGCAGGGGCGTCCACAGACGCTACCATTCCCTGTTGCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGTAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAACCGTCGTTTTTATCATCGTTATTCTCCAGAGATTAAGGTGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAgagg                                                                           <  1:142905/299‑1 (MQ=255)                                                                                                                                                                                                                              attctcccccGCCCGAGTGGCCGGGGCGTCCACAGACGCTACCATTCCCTGGTTCGCGCGTACCGCGGGGAAGACACCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTCCCTGAGAGCAGCGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTGATCATCGTTATTCTCCAGAGATTAAGGGGCACCCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAgagg                                                                           <  1:153761/296‑1 (MQ=255)                                                                                                                                                                                                                              atCCGCCCCCACCTGAGTGGCAGTGGCGTACACAGACGATACCATTCCCGGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAgagg                                                                           <  1:245662/301‑1 (MQ=255)                                                                                                                                                                                                                                  cccaccgtcggatgcgccgcggcgctcccaGACGCCACGCCTCCCTTGTTCGCGCGTCCCGTGTGGAAGACATC‑TCCGCCACACCATAAGAGGCGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAGAGGAGGAt                                                                      <  1:538852/273‑1 (MQ=255)                                                                                                                                                                                                                                         ccTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACCCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAGAGGCGGATACGCGa                                                                >  2:350901/1‑301 (MQ=255)                                                                                                                                                                                                                                                ggcagtggcGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACCCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAGAGGAGGATACGCGATGAAACg                                                         >  1:563064/1‑301 (MQ=255)                                                                                                                                                                                                                                                  caggggcGGCAACAGACGGCACAAGTCCCGGTGGCGCGCGTACGGGGGGGAAGACATCTTGCGACAAAACTGGCGGAACGGGCGGCGGAGGAGGCGCGGAGGCGGCGCGAAAGCAACTACCTGAGAGCAGAGCAGCAATGGCAACACGGAGTAAAAACGTCGGTTTTATCATCGTTAGTCTCCAGAGATTAAGGAGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAGAGGAGGATACGCGATGAAACGa                                                        <  1:352207/300‑1 (MQ=255)                                                                                                                                                                                                                                                           ccACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCCATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAAACCCAGCTAAGCCTGGTATATAACTCTGAATTAATAATCGTTTTGTCGGGAAGCGGGTAAACTTAAAGGATTTTTTCAGGGGAGGATACGCAATGAAACggattcgggttt                                              >  2:188199/1‑290 (MQ=255)                                                                                                                                                                                                                                                                gACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCCGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTCAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTTTCGGGAAGCGAGTAAACTTAACGGATATATTCAGAGGAGGATACGCGATGAAACGACTTCTGATTCTTACg                                         >  2:707477/1‑301 (MQ=255)                                                                                                                                                                                                                                                                    cTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGGAAACTTAAAGGATATCTTCAGTGGAGGATACGCGATGAACCGACTTCTGATTCTTCCggaa                                      >  1:131380/1‑298 (MQ=255)                                                                                                                                                                                                                                                                         agtCCCGGTTCCGCGCGTACCTGGTGGAAGACATCTTCCTCCACACCACACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGTCGCTAAAACACCTTACTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTa                                                                                                                                                    <  2:432888/183‑1 (MQ=255)                                                                                                                                                                                                                                                                         aTTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAAACGTGCGGCGCAGGAGGCACGGCGGCGGCGCTAAAACGACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGACAACCCCAGCtc                                                                                                                                                    >  1:432899/1‑184 (MQ=255)                                                                                                                                                                                                                                                                                                       aCATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAGAGGAGGATACGCGATGAAACGACTTCTGATTCTTACGGCACTCCTGCCGGTTGTCGGCTTTGACAAGCCCATTaaa  >  2:402848/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                        GATTTTTTGCTGTCCGGGTTGCCCTGATCATCGAGGAACATATCGCGGGTCGCTTTGGCGGGTGCCATTTCGCGGAAATCGATAGCCGTGGTATTGCCATTTTTCGAGCGGATTAACATAAAACCACCACCGCCCAGATTCCCTGCCTGCGGATGCGTTACCGCCAGCGCGTAGCCCACCGCCACGGCGGCATCAACGGCATTCCCGCCCTCCTTGAGAATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCTACCATTCCCTGTTTCGCGCGTACCGGGTGGAAGACATCTTCCTCCACACCATACGAGACGGGCGGCGCAGGAGGCGCGGCGGCGGCGCTAAAACAACTTCCTGAGAGCAGAGCAGCAATGGCCACCCGGCGTAAAAACGTCGGTTTTATCATCGTTATTCTCCAGAGATTAAGGGGCAACCCCAGCTAAGCCTGGTATATAACTCTGAATTAATCATCGTTTTGTCGGGAAGCGAGTAAACTTAAAGGATATCTTCAGAGGAGGATACGCGATGAAACGACTTCTGATTCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAAT  >  NC_000913/3586411‑3587005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.