Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,706	T→G	100%	K173Q (AAA→CAA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	4,116,706	0	T→G	100.0%	127.7	40	K173Q (AAA→CAA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  new base (19/21):  ref base (0/0):  total (19/21)

ATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATA  >  NC_000913/4116408‑4116992                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                               taaaaGCATCCAGTGTCCAAGGTGGTCTTCGCCATCCCTTCGTTCACGCACCACTGACCAAAGAGCCCGCCCTGCTGGTCACCGGCGTTCCCGGGGATTGGAATACGCGTGCCGCCTTGGCCGCCAATGTCAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGcc                                                                                                                                                                                                                                                                                              <  2:176402/299‑1 (MQ=255) aTAAAGCGGCCAGTGCCATAGGTTGGCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCGGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGcc                                                                                                                                                                                                                                                                                              <  2:324834/301‑1 (MQ=255)                           ctgtgccatccatgacgtCACGCACAACGGACTAAACAGCGCGGCATGTTGGTCACGCGCGATACCGGGGATTGGAAGACGCGTGCCGCCTTTGCCGACAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCa                                                                                                                                                                                                                                                                    <  2:28489/284‑1 (MQ=255)                                                     aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAGTCGGTCACATGGACACGGCCCTGATGCATTTGCCGTATAAGCCACGTAGCAACCCTACCAAACAGCACTTCACCACGACGTGCAc                                                                                                                                                                                                                                           >  1:183686/1‑300 (MQ=255)                                                        gACCAAACAGCGCGGCCTGCTGGTCAACGGCGATCCCGGAGATTGGAATACGCGTGCCGCCGTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCt                                                                                                                                                                                                                                        <  2:10263/300‑1 (MQ=255)                                                                         ttgtgGTCACCGGGGATGCAGGGGATGGGAAGAAGCGCTGGGGCGGTGCCGGCAGTGTGTGTCGGACCGGATACAGCGAAAGAAAGAAGCACTGCTGGCAGCATCTCGCGCGGAATATACAGCACGTCCAGCATTTTGGCGTCCCAGTTCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTGATAGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACAt                                                                                                                                                                                                                      <  1:5454/298‑1 (MQ=255)                                                                                   cGGGGATCCCGTGGATTGGAATACGCGTGCCGCCTTGGCCGCCAATGTTAGTCTGACCGGATACTTCGGAAGAACGACGCACTGCTGGGAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGAt                                                                                                                                                                                                            <  1:419548/301‑1 (MQ=255)                                                                                               agaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTACCCACGACGTGCACGCCCGCGAGAGCCTTCTACATGGTTGAGGACCCCCTTCACttt                                                                                                                                                                                                >  2:359350/1‑301 (MQ=255)                                                                                                          cgcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCAATTCCAGCATTTTGTCGTCCCAGTCCAGCGCATGGATGTTGAACAACATGGTACGAGAGGCTTTGGGGTAATCGGTCACATGGACACGCCCCTTAGTCATTTGCCAGCTAAGCTACGGCACAACCGTACCAAACAGCAATACACCCCGACGTGCACCCTCGCGAGCACCTTCCCCACGGTCCCGGCTCCCCCTACCTTTGGTCCCAGaaa                                                                                                                                                                                     >  1:73271/1‑301 (MQ=255)                                                                                                                               aTGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCCGTCCAGGGTCTGGATGTTGAACAACATGGTACGAGAGTCGTTGGTGTAATCCGTCACATGGACACGGCACTGAGTCATTTGCCAGCTAAGCCACGTATACACCGTAACAAACAGCACTTCACCACGACGTGCCCGCCCGCCATCGTATTCCTCCTGGTCGAGGAGCCACTCACCTCTGGAGCCAGAAACTGACGGGTCAACCACCaaca                                                                                                                                                                >  1:12770/1‑298 (MQ=255)                                                                                                                               aTGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGCGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAACCAGCAATTCACCACGACGTGCACGCt                                                                                                                                                                                                                                        >  2:18400/1‑229 (MQ=255)                                                                                                                               aTGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCt                                                                                                                                                                                                                                        <  1:18400/229‑1 (MQ=255)                                                                                                                                tGTTAGTCTGACCGTAGACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAcc                                                                                                                                                               <  2:38455/301‑1 (MQ=255)                                                                                                                                  ttcgtgctactggtaaactacgaagtaCGACTCACTTTGGGCAGCAGTTCTGGTGGAATATCCCGCACTCCCAGCATTGTTTATGCACAGTCCAGCGTATGGATGTTGAACAACATTGGACGAGAGGCGTTTGTGTAATCGGTCACATGGACACGGCCCTGAGTCAGTGGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAAGTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAcc                                                                                                                                                               <  2:204111/278‑1 (MQ=255)                                                                                                                                          aCCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcaccac                                                                                                                                                                                                                                                    >  1:101475/1‑206 (MQ=255)                                                                                                                                          aCCGTATACTTCGGAAGAAAGACGCACTTCTGGCAGCAGCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGGCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcaccac                                                                                                                                                                                                                                                    <  2:101475/206‑1 (MQ=255)                                                                                                                                          aCCGTATACTTCCGAAGAACGACGCGCCTCGTGCATCATACCGCGCGGAAGCTTCAGCACTTGCAGCATTTTTGCGTGCCTGGCCAGGAGCTGGAGGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACCTGGACACGGCCCTGATTCATTTGCCAGACAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTGTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctg                                                                                                                                                     <  2:138740/301‑1 (MQ=255)                                                                                                                                                  cTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTTCACGCTCGCCACAGCCTTCCACATGGTCGAGGATCCACTTCCCTTTGGTGCCAGAAAAGTACGGTTCCATCCCCAGCCCGGTTTTGCTTCGGAcata                                                                                                                                             >  1:27153/1‑301 (MQ=255)                                                                                                                                                     cGGCAGCACGACGCACCTCTTGCAGCATCTCTCGCGGACTAGCCAGCACTGCCCGCAGTTTGTCGTCCCAGTCCAGTGTATGGAGGTGGAACAACATGGTACGAGATGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTCTCAACCGTACCAAACCGCAATTCACCACGACGTGCACGCTCGCGAGAGCTTTCCACATGTTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtc                                                                                                                                          <  2:183687/301‑1 (MQ=255)                                                                                                                                                       gaagaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCCCATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGCCAATCCCCAGACCGGTATTCCTGCGGATGTAAtctt                                                                                                                                        >  1:99982/1‑301 (MQ=255)                                                                                                                                                             ctacgCACTTCTGGCAGAATCTCGCGCGGAATATCCATCACGTCCAGCATTTGGGCGTCCCAGTCCAGGGTATGGAGGGTGAACAACAGGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAAc                                                                                                                                  <  1:28838/299‑1 (MQ=255)                                                                                                                                                                  cACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCt                                                                                                                                                                                                                                        <  1:21962/194‑1 (MQ=255)                                                                                                                                                                           aagcaGCTCGCGCGGAATAGCCAGCACTTCCAGCAGTTTGTCGTCCCAGGCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGGAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTaaa                                                                                                                    <  2:27153/300‑1 (MQ=255)                                                                                                                                                                                                                                gggTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTAACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGCAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGAAGATTTATGCGGTACGCCGGCACTGCCAGCCAATGGCGTTatagat                                                               >  2:43965/1‑301 (MQ=255)                                                                                                                                                                                                                                           ttGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCACTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCAACATGGTCGAGGATCCACTTCACTTTGGTGACAGAAAAGTACGGGTCAATCACAAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGATTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTctagatctgtcgctaca                                                    >  2:94371/1‑290 (MQ=255)                                                                                                                                                                                                                                           ttGAACAACATGGTACGAGAGGCGTGGGTGTAATCGTTACCATGGACACGGCCCTGAGTCATGTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTAACCACGACGTGCACGCTCGCGAGAGCCTTCCACATCGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATCATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTatagat                                                               <  1:94371/290‑1 (MQ=255)                                                                                                                                                                                                                                                caacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCCCGTTTTAAATGCTCGCCGATTTCTGCGGTACGACGGCACTGCCAGACCATGGATTTATAGATAGGATTTCCGGttt                                                 >  1:374062/1‑299 (MQ=255)                                                                                                                                                                                                                                                  acaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTAGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCTGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTTCGGCTATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCACACAATGGCGTTATAGATAGTCTTGACGGTTTCttt                                             >  2:434974/1‑301 (MQ=255)                                                                                                                                                                                                                                                     tGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTAGAGGATCCACTTCACTTTGGTGCCAGCAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGGTTTCTGCGGTACGACGGCACTGCCTGACAATGGAGTTATAGATAGGCTTGCCGGTTTCttttt                                           >  2:103724/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 cGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCCCCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGCGGATCCACTTCACTTTGGTGCCAGACAAGTACGGGTCAATCACCCGCCCGGTATTTCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTTCGGTCCGCCGGCCCTGCCATACAATGGCGTTATATATAGGCTTGCCGGTTTCTTTTTCCCAGACAATCGt                              >  2:106958/1‑301 (MQ=255)                                                                                                                                                                                                                                                                     ggTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCAAGAGCCTTACACATGTTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAACCTTCTAAACCGGCACGTTTTCACTGCTCGCAGATTTCTGCGGTACGCCGTCACTGCCAGGCAATGGCGTTATAGAGCAGCTTACCGGGTCCTTTTTCCCTGCCACTAGTGGtt                          >  2:61941/1‑301 (MQ=255)                                                                                                                                                                                                                                                                     ggTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGGATCAACCGTACCAAACAGCAATTCACCACGAAGTGCACGCTCGCGAGAGCCTTCCACATGGTCGCGGATCCACTTCACTTTGGTGCCAGACAAGTACGGGTCAATCACCAGACCGGTATTGGTGCGGATATAATATTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGCCGATGGAGTTATGGTAAGGCTTGCCGGTTTCTTATTccccgcgaatgatggtt                          >  2:331657/1‑287 (MQ=255)                                                                                                                                                                                                                                                                               gTCACATGGACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGGAAAGTACGGGTCAATCAACAGCCCGGTATTGATGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTAGCAGATTCCTGAGGTACGAAGGAACTGCCAGACACTGACGTTATAGATAGGCTTGCCGGTTTGTTTTTCCCATACAACAATGTGTTCAGACTGGt                >  2:113168/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                     tGGACACGGCCCTGAGTCATTTGCCAGAGAAGCCACGTATCAACCAGACCAAACAGCAATACACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACGTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCTGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTaa          <  1:113168/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                      ggACACGGCCCTGATTCATTTGCCAAATAAGCCACGTATCAAACGTACCAAAAAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCGTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTaa          <  1:43965/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                       gACACGGCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATCTAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGTTTTCTGCGGTACGACGGCACTGCCAGACACTGGCGTTCTAGATCGGCTTGcct                                                     >  1:224444/1‑255 (MQ=255)                                                                                                                                                                                                                                                                                       gACACGGCCCTGAGTCATTTGCCAGAGAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCg                                                     <  2:224449/256‑1 (MQ=255)                                                                                                                                                                                                                                                                                       gACACGGACCCCAGTCATTTGCCAGATACTCCACGGATCAACCGTACCCAACACCACTTCACCACTACGTGCACGTTCGCGAGAGCCTTCCACAGGGTCGCGGTTCCAATTCACTTTGGTGCCAGAAAAGTACGGGTCAATCCCCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAAt         <  2:99982/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                        acacGGCCCGGAGTCATGTGCCAGTTAAGCACCGTATCAACCGTACCAAACAGCAATCCACCACGCCGTTCCCTCTCGCGAGAGCCGTCAACATGGGCGAGGATCCACTTCACTTTGGTGCCAGAAAACTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGGTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc       <  1:61941/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                             gCCCTGAGTCATTTGCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGCTTTCTGCGGTACGCCGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGGGGTTTCACGCTGGTTCGTAATACCGATa  >  1:194940/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                               ATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATA  >  NC_000913/4116408‑4116992

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.