Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,560,460	A→T	20.4%	pseudogene (146/150 nt)	glpR ←	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,560,460	0	A→T	20.4%	17.2 ‑2.0 ‑0.2	50	pseudogene (146/150 nt)	glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  new base (6/4):  ref base (11/25):  total (20/30)
Fisher's exact test strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.91e-01

CATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATC  >  NC_000913/3560162‑3560751                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                    cATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAATGCAGCGCCGCACCGCCATGC‑TG‑GCGCAGGAt                                                                                                                                                                                                                                                                                                    >  2:358050/1‑301 (MQ=255)       ggtggcggaagtggacACCGGAGAGAAGGATGCGAAAATCTGCGTTGACAACAAACGGGTTGGAAACGGTGAGATGGTTGGTGACAATGCGCAAATTGCGGGGATGGAGCAGTGCGGGCGCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGGTCCGCCACTTTGCGGGCGATGCGCTCTGTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGTTCAGGt                                                                                                                                                                                                                                                                                               <  1:399702/285‑1 (MQ=255)              tAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGGTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACGGGAGGGCATCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTGGcgc                                                                                                                                                                                                                                                                                        >  2:40507/1‑298 (MQ=255)                        ccGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCCCCGCACCGCCCTGC‑TG‑GCGCAGGATCCGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                                                                                               >  1:348981/1‑299 (MQ=255)                          ggCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCCTTGGGGATTTGCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTCACCGCACTGGACGGTAGCGCCGCCCCGCCATGT‑TG‑GCGCCGGTACAGGTTTTGCTCCGCCCGCTcttcg                                                                                                                                                                                                                                                                           >  2:376692/1‑297 (MQ=255)                           gctaaaaTGATGCGAAAAGCGTCTTTGCACATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATCCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                                                                           <  2:348990/295‑1 (MQ=255)                              agaATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGC‑TG‑GCGCAGGCTCAGGTTTTGCTCCGCCAGGTCCttacggt                                                                                                                                                                                                                                                                       >  2:186204/1‑296 (MQ=255)                                                tcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGC‑TG‑GCGCAGGCTCAGGTTTTGCTCCGCCGGCTCATTGAGGTCCCGGCGACTGGTCTggg                                                                                                                                                                                                                                                     >  1:423470/1‑299 (MQ=255)                                                         atgagcGTGTTCGATACGTTGAGATGGTCGTTGACAATGCTCAAGTGGATGTGATTGCGCAGTGCCTGTGCCACCGCTGCTGGTTGTGTGCCGATGTCGATAACCAGCGTCGAGCCATGGGGGAGTTGTTCTGCCCAGTTGCGGGCGATGCGCTAGTTTTCGTCGGTCTGGGTGGCCTGGCGATCGGGCCACGGCGGGTTACCCGAACTGGAAGGCAGCGTCGCACCGCCATGA‑TG‑TCGCAGGAGCAGGTTTGGCTCCGCAAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                                                                            <  2:390340/296‑1 (MQ=255)                                                              cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGA‑TG‑GCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGGGGGCTGACGGAGaaa                                                                                                                                                                                                                                       >  1:421519/1‑301 (MQ=255)                                                              cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGG‑TG‑GCGCAGGGTCCGGTTTTGCCCCGCCAGCTCCTTGAGGTCCCGGCGAATCGCCGGCGGGCTGCCGGGGaaa                                                                                                                                                                                                                                       >  2:74461/1‑301 (MQ=255)                                                                            tggggtgggtgGTGACGATGGGCAAATGGCGGTGATTGGGCAGGGCGTGCGGTAGCGCTTCCGGAGTGGTGCCGATATCGGTAAATAGCGTCTGGCCATTGGGGATGTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                                                                                                                                                                                         <  1:73068/292‑1 (MQ=255)                                                                            tgggaTGGTTGGGGCAAAGGCGCAAAGTGGTGGGATTGAGCAGTGCGTGGTCTACCGGTTCCGGCGTGGGGCCGAGAGCGAAAAACAGCTGCGAGCCATTGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                                                                                                                                                                                         <  1:217352/298‑1 (MQ=255)                                                                            ggtggttgagggttcaaaTGCGGAAATTGTTGTGGTGGAGCAGGGGGGTGGCTACCGTTGATGGGGTGGTGCCTATATCGATAAGAAGGGTCGTGCCAGTGGGGATTTGCTCCGAGACGTGGCGGTGTACGCGCTCTTGTCCTTAGGTCTGGGTGGCCTCGCGATCGCGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                                                                                                                                                                                          <  1:433727/284‑1 (MQ=255)                                                                                         gACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCTTTGCGGTCCCGGCGAATAGTCTGCGGCCTGACGGAGGAATGATCTACCCGCCCTTCGGTCATGACa                                                                                                                                                                                                            >  2:238145/1‑301 (MQ=255)                                                                                              tgcgtaattggatgtgATTGGGGAGTGCGTGAGCTACCGGTTCTGGCGGGGTGCCGAGAGCGATAAACAGGGTCTAGCCGGTGGGGGTTGGCTCCGCCACTTGGCGGGCGATGCGCTCTGTTTCTGCGGTCTGGGTTGCCTTGCGATCGTGCCACGGCGTGTTAAGCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                                                                                       <  1:241798/289‑1 (MQ=255)                                                                                              tgagcaatgtgttgtgATTGGCCAGTGCGGGCGCGACCGCTCGCGGCGGGGGTCCGATTGTGATAAACAGCGTCGAGCCATGGGTGAGTCGCTCCTCCACTGTGCGTGCGATGCGCTCGTGTTCTTCGGGCTGGGTGTCCGTGCGATGGTGCCACTGCGTGTTAATCGAACTGGAAGGCAGCGCCGCACCGCCATCA‑GG‑GCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                                                                                       <  1:187200/289‑1 (MQ=255)                                                                                                ggcAAATTGCTGTGATTGGGAAGGGCGGGCGGGCGGGGTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTGTGCGGGCGATGCGGTCTTTTTCTTCGGTCTGGGTGGCCTTGCGTTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                                                                                      <  1:186203/299‑1 (MQ=255)                                                                                                cgtcatgtttggtgATGGAGCAGTGCGTGTGCTACCGCTGCGGGGGTTGTGCCGATCGCGATAAACAGCGTCGAGCCATGGGGGCTTTCCTCCGCAACTTTGCTGGCGATGCGCTCTTTTTCTGCGGTCTGGGGGGCCTTCCGATCCTGCCACGGCGTGTTAACCGAACTTGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                                                                                      <  1:171448/289‑1 (MQ=255)                                                                                                cgtaaaTTGGCGTGCTTGAGAAGTGCGTGGGCAACGCGTTCCGGGTTGGTGACGATACCGGTAACCAGCGGCTTGCCATTGGGGATTTGCTCCGCGACTGGGCGGGCTATGCGCTCTTTTTCTTGGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCACGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                                                                                      <  1:66576/297‑1 (MQ=255)                                                                                                cggtaaTTGCTGTGATGGGGCAGTGTGTGCCGTACCGCTTCCGGCGTGGGGCCGATATCGATAACCAGCGTCGAGCCATTGGGGATTTGCTACGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                                                                                      <  1:304583/296‑1 (MQ=255)                                                                                                 gcAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGACATGCTCTAACAGCTCTTCGGTAc                                                                                                                                                                                                                 >  2:404175/1‑288 (MQ=255)                                                                                                 gcAAATTGCTGTGATTGAGCAGGGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                                                                                                                                                                 <  1:404166/288‑1 (MQ=255)                                                                                                           tgggttGAGCAGGGGGGGGGGTACGGGTTCCGGCGGGGGGCCGATATCGATAAACAGCGTCGGGCAATTGGGGAGTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGAACTGACATAACCCTGCTGTTTAAcc                                                                                                                                                                                          <  1:358041/297‑1 (MQ=255)                                                                                                           tggcttggtcAGTGCGTTCGGGTCCTCTGCCGGCTTGGTGCCTCTATCGATAACCACGCTCTGGCCATTGGGTATTTGCTCCGCCACTTTGCGGGCTATGCGGTCGTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGAACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                                                                                                                                                                          <  1:162958/292‑1 (MQ=255)                                                                                                              gtttggagggcggggcggtcccGGTTCCGGCGTGGTGCCGAGTTCGATAAACGGCGTCTGGCAATTGGGGATTTGCGCCGCCACTTTGCGGGGGATGCGCTCTGTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                                                                                                                                                                        <  1:376683/280‑1 (MQ=255)                                                                                                                      gtgcgtgcGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGCTTTGCTCCGCCACTTTGCGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACCGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATCA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGCATAGTCCGGGGGCTGACGGCGAAATTCTCTACCAGCTCTTCGGTACTTCCATAACGCTGCTGTTCAACCCGTTCGaaaaa                                                                                                                                                                               >  1:288276/1‑300 (MQ=255)                                                                                                                         cgtgcgAGACCGCCTGCGGGGTGGTGCTGAGATGGCTAAAGACCGTTTAGCCATGGGGGATTGGCGCCGCACGTTTGCGGGCGAGCCGCTCTTTTTCTTCGGTGTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCATGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCTGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                                                                             <  2:10348/300‑1 (MQ=255)                                                                                                                          gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGA‑TG‑GCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATCGTCTGCGTGCTGACGTAGAAATGCTCTCCCAGCTCTTCGGTGCTGACATAACCCTGCTGGTTAACAATTTCGAAAATGCCg                                                                                                                                                                           >  1:134751/1‑301 (MQ=255)                                                                                                                                       ccccgCGTGGTGCCGATATCGAGAAAGAGCGTCGAGCGATGGGGGATTTGCTCCGCCATGTGGCGGGCGATGCGCTCTGTTTCTGCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                                                                              <  1:142919/297‑1 (MQ=255)                                                                                                                                         cggggtggttccgcttggaaTAATCAGCGCGTAGTCATTGGGGTTTTGGTCGGCCACTGTGCGGTCGAGGGGCTGTTGTTGTTAGGTCTTGGTGTCCGTGCGAGCGGGCCACGGTGTGGTAACCGAACTGGAAGGCAGCGCCGCAACTCCATGA‑TT‑GCTCAGGATCAGGTTTTGCCCCGCCAGCTCATCGAGGTCGCTGCGAATAGTCGGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt                                                                                                                                                            <  2:423479/282‑1 (MQ=255)                                                                                                                                                        ttggatgaaggcggtcggggctttGGGGTTTGGCTCCGCCACTTTGCGGGGGAGGCGCTCTGTTTCTCGGGTATGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                                                                                             <  1:40507/279‑1 (MQ=255)                                                                                                                                                         tCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGCATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACATCACCCTGCGGTTTAACCAGTTCGATCATACCGTTGTGCCTTTGTGTTTTTTTCATTTATAAAt                                                                                                                                            >  2:308052/1‑301 (MQ=255)                                                                                                                                                         tCGATAAACAGCGTCTAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTACTGACATAACCATGCTGGTTCACCCTGTGGATAATACCGTTGTGCCGTTGAGTTTGTTTCTTTtatacaa                                                                                                                                            >  2:28934/1‑298 (MQ=255)                                                                                                                                                                 catgtccggtccATTGGGGTTGTGCTCCGCCCCTTGGCGGGCGATGCGCTCTTTTTCTTCGGGTTGGGTGGCCTTGGGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCGATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGaa                                                                                                                                    <  2:288282/291‑1 (MQ=255)                                                                                                                                                                          gCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTCCGGGCTGACGGCGAAATGCCCTACCAGCTGTTCGGTCCTGACATACGCCTGCTGGTTAACCAGTTCGCTAATAACGTTCTTACGTTTTGTTTGGTTCATTTTTCAATCCCTGGAATTAtttt                                                                                                                             >  1:292013/1‑299 (MQ=255)                                                                                                                                                                                       tGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTCGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTATACCAGATCTTCGGTACTGACCTACCCCTGCTGTTTCACCAGTTCGATACTACCGTTGTGACGTTGTGGTTGTTTAATTTTATAATCTCTGGTATTATTTTGGTTTTCTCTCATTg                                                                                                              >  1:330819/1‑301 (MQ=255)                                                                                                                                                                                                                        tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg                                                                             >  2:76426/1‑301 (MQ=255)                                                                                                                                                                                                                        tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTAACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCCGTTCGATAATACCGTTGTGACGTTTTGTTTGTTTAATTTATAAATCCCTGGCTTTCTTTTCGTTTTCGCTCATTGCGCGAATCCACAAAAGCCATCGCTAAACCCCCg                                                                             >  1:300530/1‑301 (MQ=255)                                                                                                                                                                                                                          tgcgggttGGGGGGCTTTGCGTGCGGGCCACGGGGTGTTAACCGAACTGGAAGGCAGAGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGc                                                                           <  1:308044/294‑1 (MQ=255)                                                                                                                                                                                                                                tgggttGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc                                                                      <  1:28934/295‑1 (MQ=255)                                                                                                                                                                                                                                       gCCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCGGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaacaa                                                                                                 <  2:15092/266‑1 (MQ=255)                                                                                                                                                                                                                                               ttagtagggtcaggggggacacGAACCGGGAGGCCGCGCCGCACCGCAATGA‑TG‑GCGCAGGTTCCGGGTTGCACCCGCCAGCTCATCGCGGTCGCGGCGAATTGTCGGCGGGCCGACGGAGAAAGGCGCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTGAACCAGTTCGATAATACCGTTGTTACGGTGTGTTTGTCTCATTTATAAATCCCTGGAATTATGTTCGTTGTCGCGCATTGAGCGAATCAACAAAAGCCATCTCTAAACCCACGGCTAACCCGGCGATGTGTGCTc                                                       <  1:14625/279‑1 (MQ=255)                                                                                                                                                                                                                                               aTCGGGCCCCGCGGTGTGAACCGAACGGGAAGGCAGCGCCGCACCGCAATGA‑TGCGG‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTc                                                       <  1:198951/300‑1 (MQ=255)                                                                                                                                                                                                                                                            gtgtTAACCGAACTGGAAGGCAGCGCCGCACAGCAATGA‑TG‑GCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAGATGCTCTACCAGCTCTTCGGTACTGACGTAACCCTGCTGTTTACCCAGTTCGCTAATACCGGTGTGACGTTGTGTTTGTTTCATGTATACATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGGATCAACAAACGCCCTCGCTAAACCCAAGGCTAACACGGCGCTGTGTGCGcctgtctctgctcc                                          >  2:14625/1‑288 (MQ=255)                                                                                                                                                                                                                                                             tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                                                                                             >  1:360360/1‑200 (MQ=255)                                                                                                                                                                                                                                                             ggtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                                                                                             <  2:360369/199‑1 (MQ=255)                                                                                                                                                                                                                                                                   cagAACTGGTAGGACGCGCCGAACCGCCATGAGTG‑CG‑CAGGACCAGGTGTTGCTACGCCATCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTAACGGAGAAATGCTCTACCAGCTCTTCGTTACTGACATAACCCTGCTGTTTAACAAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATccc                                    <  2:300538/297‑1 (MQ=255)                                                                                                                                                                                                                                                                                                     tggcgcAGGTTGAGGTTTTGTTCGGCTATGTCATTGTGGTGGCGGAGAATAGACTGAGGGCTGACGGAGAAAGGCTATACCAGCTCTTCGGTACTGACATAACCCTGGTGTTTAACCAGTTCGATATATCCGTTGTGACGTTGTGTTTGTTTCATTTATAAAGCCCTAGAAGTATTTTCGTATTCGCGCAGTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCAGCCGGCGACAATCCAGATc  <  1:109178/296‑1 (MQ=255)                                                                                                                                                                                                                                                                                                      gggggaGGATTAGGGTTTGTCCCACCAGCTCATGGAGGTTGCGGCGAATAGTCGGCGGGCTGACGGAGAAATGTTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATc  <  1:238140/295‑1 (MQ=255)                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                    CATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TG‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATC  >  NC_000913/3560162‑3560751

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ 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.