Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,529,218	C→T	100%	pseudogene (1317/2037 nt)	rhsE →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	1,529,218	0	C→T	100.0%	78.2	42	pseudogene (1317/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  new base (16/22):  ref base (0/0):  total (16/26)

CGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGAC  >  NC_000913/1528943‑1529508                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    aggagaaagggcgaGCCGCGCAGCGTGGCGGTGACGCTCCAGGAGGAAGGGGGTGCGGAGGGCAACGGCGTGGGGTTCCCGGCGGAACTGGGGCGGCTTGTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGGGGCTTGCGCAGTGCGGGCTGACGGGGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACgg                                                                                                                                                                                                                                                                               <  1:38068/285‑1 (MQ=11)   ggAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACCGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCCGCGCAGACGGCAATACGGCGTgg                                                                                                                                                                                                                                                                          >  2:350159/1‑301 (MQ=2)                   cgcAGCCGGGGAGGGACGCTCCAGCAGGAAGGGAGTGAGACCGGCCACGGTGTGGGGTTCCCGGCTGAACGGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGGTAGCA‑GTGAAAGCCGGGCTTGGCTTGCGCAGTGCGTGCGGACGGTTGAGCAACTGGCCAG‑ACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                                                                                                                                                                                           <  2:268403/300‑1 (MQ=11)                      cgtcgggCAGGGACGGTCGAGGAGGAAGGGAGTGAGAACGGCAGCGGCGTGGTGGGCCCGGAGGGAGTGGGGCGGCTGCTGGACAGGCGGGAGGAAGAAAGCCGGGCAGACCGAGGGAGAA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGCCGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                                                                                                                                                                       <  1:394084/296‑1 (MQ=11)                            gCAGAGGCGCGCCAGCAGGAAGGGAGGGAGAACGGCCAGGGCGTGGGGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAAGCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGc                                                                                                                                                                                                                                                  <  1:350150/300‑1 (MQ=255)                                   cGCTCCAGCAGGAAGGGAGTGAGAACGGCAACGGCGTGGTGTTCACGGCTGAACTGGTGCGGCTGCTGGACAGGATGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGACGGGCGTGGCTTGCGCAGTGAGGGCTGACGGTGGAGGAACTGGCAAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGGAGCGTTCGTTTGTCACGCTGCGACCACCGGGGCCTCACGCTGGGGCTTATCAGAGTAGACGGCGATACGGCGTAGCGCTGTGAGTGTGATGTATGAGGCCACCAGcg                                                                                                                                                                                                                                          >  2:151583/1‑300 (MQ=11)                                         atcagGAAGGGAGTGAGAACGGCCACGGCGTGGGGTTCCCGGCTGAACTGGTGCGGCTGATGGACAGGCGGGAGGAAGCATTCCGGGCAGACCGCGTGGGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                                                                                                                                                    <  2:266551/299‑1 (MQ=255)                                          gcagGCAGGGAGTGCGAACTGCTACGGCGGGGTGTTCCCGGCTGAACTGCTGCTGCTGCTGGACAGGCTGGAGGAAGAACTCCGGGCAGCCCGCGTGAGCA‑GTGAAGGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGTTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATTGGGCAACCAGCTTAATga                                                                                                                                                                                                                                    <  2:346303/300‑1 (MQ=11)                                          gcagGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAACAGTTCATTTTTACCACTGCGACCACCGGGGCCCGCCGCTGGCGCTTATCAGCGAAGACGGCAATCCGGCGTGGCGCGGGGCGGAGGATGAATGGGGGAACCCGACTAATgaa                                                                                                                                                                                                                                   >  1:277242/1‑300 (MQ=11)                                             ggAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGCCAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTAAAGAAACTGGCCAG‑ACAGGTGGAGCCGGAGTCCGCCCGGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCGTGCGGCTGGCGCTTATCAGCGAAGACGGCAATACGGAGTAGCGCGGGGAGGAAGGTGAATGGGACAAGCCGCTTCATgaggag                                                                                                                                                                                                                                >  2:334075/1‑301 (MQ=11)                                               aatgtggagagaACGGCGTGGGCGTGGGGTTACCGGGTGAACTGGTGCGGCTGATGGACAGGCTGGAGGAAGAAGTCCGGGCAGACCGCGGGAGTA‑GGGCAAGCCGGGCGTGGCTTGCGAAGTGCGTGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTTGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCTCCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCTCGGTGAGTATGATGAATGGGGCAACCAGCTTAATGAGgaga                                                                                                                                                                                                                               <  2:404813/292‑1 (MQ=11)                                                                     cGTGGGGTTCCCGGCTGAAATGGTGCGGCTGCGGGACAGGCTGGAGGAAGAAAGCCGGGCAGACCGCGTGGGCA‑GTGAAAGCCGGGCGTTGCTTGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGACGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGcc                                                                                                                                                                                                        <  2:250330/300‑1 (MQ=255)                                                                      gtggtgTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGCGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCGGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGTAATACACACCGGCGCGAAACGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCCCTGGCGCTTGTAAGCGAACACGCCAATACGGCGTGGCGCGGGGGGTGTGGTGGATGCGGCAACACGCATCATGAGGACAACCCGCCG‑CACATGCACCCCCCGt                                                                                                                                                                                                      >  2:311953/1‑301 (MQ=11)                                                                         gtgtTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTAACCGCGAAGACGGCAATAAGGCGTGGCGCGGGGGGGATGACGACTGGGGCAACCAGCTTACTGAGGAGAACCCGCAT‑ACCCTGCCCCAGCCGTAc                                                                                                                                                                                                    >  2:156593/1‑300 (MQ=14)                                                                                cGGTTGAACTGGTGCGGCGTCGGGACAGCCTGGGGGAAGAAACCCGTGCAGGCCG‑GTGGGGACGTGAAAGCCGGGCGTGGCTTTCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGcc                                                                                                                                                                                            <  2:277248/301‑1 (MQ=14)                                                                                                ggtgtggggcgggCGGGGGGATGAAACCCGGTCAGAGGGGGGGGGGA‑GTAAACGCCTGGGGGGGCGGGCGCGGTGTGGGCGGACTGTGGCGCAACTGGACAG‑AAAGGGGGAGCCGGAAGATACACCGGCGCGAAAAGTTCATTGTTACCACGGCGACCACCGGGGCCTGCCGCGGGCGCTTATCAGCGAAGACGGCAAGACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGAACCAGCCGTACCGTCTGCCAGGGCAGCAGCatga                                                                                                                                                                             <  1:151582/289‑1 (MQ=14)                                                                                                gctgcgGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGGTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatga                                                                                                                                                                             <  2:52466/300‑1 (MQ=255)                                                                                                       atgtgtataaagggcagaaTCGGGGCAGACCGGGGGAGCA‑GGGAAAGCCGGGCGTGGCTTGCGAAGGGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGAGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag                                                                                                                                                                      <  1:103364/283‑1 (MQ=14)                                                                                                           ggCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCCACCAGCTTAATGAGGAGCACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCggg                                                                                                                                                                 >  2:69811/1‑301 (MQ=255)                                                                                                           ggCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAC‑CACCTGCCCCAGCCCTCCCGTCTGCCAGGGCAGCAACACGGTGGGGAGGCggg                                                                                                                                                                 >  1:143481/1‑301 (MQ=255)                                                                                                               ggaggaAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTCTCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCGTAATGAGGAGGAGCCGCATACACTGACACC‑GGCGTACCGTCTGCCAggggcacagcaaggaggagagacgggg                                                                                                                                                                >  2:393405/1‑274 (MQ=11)                                                                                                                        acgcgggggagagcgggggggCA‑GGGAAAGCCGGGGGTGGCGGGTGGAGTGCGGGCTGACGGGGGGGCAACTGGCCAG‑ACAGGTGGGGGCGGAATACGCACCGGCGCGAAAAGTTCATTTTTACCACGGCGACCACCGGGGCCTGCCGCTGGCGCGTATCAGCGAAGACGGCAATACGGCGGGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAc                                                                                                                                                    <  1:398810/281‑1 (MQ=14)                                                                                                                        aaTCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGAGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTACTGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCAAGGGCAGCAGCAGGATgaggagctgtcccttacacacat                                                                                                                                                     >  2:103364/1‑283 (MQ=14)                                                                                                                         aTCCGGGCAGATCGGGGGAGCA‑GTGAAAGCCGGGTGGGGCGTGCGGAGTGGGGGCGGACGGTGGAGCAACTGGGCAG‑ACAGGTAGAGCCGGAATCCACAGCGGCGCGAAAAGTTCATTGTTACCACAGCGCCCACCGGGGCCTGCCGCGTGTGGTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCTGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATaa                                                                                                                                                     <  2:45699/299‑1 (MQ=14)                                                                                                                           ggcgggcagtccgggggggcgtggAAGGCCGGGGGTGGCTGGCGCGGGGCGGGCTGCGGGGGGGGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGGTCATTTGTACCGATGCGACCACCGGGGCCTGGCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAc                                                                                                                                                    <  1:135089/276‑1 (MQ=14)                                                                                                                                           ggggggggAAGGCGGGGCGGGGGTGGCGGAGGGCGGGGTGACGGGGGAGCAACGGGCCAG‑ACGGGTGGGGCCGGAATACAAACCGGCGCGAAAAGGTCATTGTGACCACTGCGACCACCGGGGCCGGCCGCTGGCGCTTATCAGCGAAGACGGCAAGACGGCGGGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGa                                                                                                                                   <  1:406286/293‑1 (MQ=14)                                                                                                                                           ggggaggcaaaGGCGGGGGGGGGTGGGGCAGTGTGGGCGGACGGGGGAGCAACTGGCCAG‑GAAGGTGGAGCCGGAATACACACCGGCGCGAAAATGTCATTTTTACCACGGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTGTAACCGTCACCGGCACTACGa                                                                                                                                   <  1:428184/292‑1 (MQ=14)                                                                                                                                                   agagtggggaggggCTTTGGCAGTGCGGGCAGACGGGGGA‑CGCCAGGTCAGTACACGTGGAGCCGTAAGACACAGCGGCGCGAACAGTTCAGTGTTACCACTGGGACCAGCGGGGCCTGCCGCGGGCGCGTACCAGGGAAGACGGCAATACTGCGTGGGGCGGTGAGTATGTTGAATGGGGCAACCATCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGc                                                                                                                           <  1:197562/288‑1 (MQ=14)                                                                                                                                                                  ttGCGCAGGGCGGGCTGGGGGGGGAGCAACTGGCCAG‑ACGGGTGGAGCCGGAAGACACCGCGGCGCGAAAAGTTCATTTTTACCACGGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCAc                                                                                                           <  1:156592/301‑1 (MQ=255)                                                                                                                                                                  tggggcggggggggGTGAGGGGGGGGCAACTGGCCAG‑GAAGGGGGAGCCGGAATAAACACCGGCGCGAAAAGTTCATTTTGACCACTGCGACCACCGGGGCCTGCCGCGGGCGCGTATCAGGGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCAc                                                                                                           <  1:392648/290‑1 (MQ=14)                                                                                                                                                                   gcgggagggggggCGGACGGGGGGGGAACGGGCGAG‑GAAGGGGGGGCCGGAATGCACACCGGCGGGAAAAGTTCATTGCTACCGGCGGGACGACGGGGGCGGGGCGCTGGCGCTGAGCAGCGAAGACGGCAATACGGCGTGGCGGGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGGCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCAc                                                                                                           <  1:232415/290‑1 (MQ=14)                                                                                                                                                                       tgaggggggCTGACGGGGGGGCAACTGGGCAG‑GCAGGGGGAGCCGTAAGACCAACCGGGGAGAAAAGTTGATTGTGACCACTGGGACCACCGGGGCCGGCCGGTGGCGCTGATCAGCGAAGACGGCAATACGGCGGGGCGCGGGGAGTAGGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGGAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCg                                                                                                       <  1:403881/294‑1 (MQ=17)                                                                                                                                                                       cAGTGCGGGCTGACGGTGGTGCAACTGGCCAG‑TCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTTCGACCACCGTGGCCTGCCGCTGGCGCTTATCAGCGa                                                                                                                                                                                                                                                                                             <  2:282609/119‑1 (MQ=11)                                                                                                                                                                       cAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCCCCGGGGCCTGCCGCTGGCGCTTATCAACGa                                                                                                                                                                                                                                                                                             >  1:282603/1‑119 (MQ=255)                                                                                                                                                                         gTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAACGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTGTGATGAATGGGGCAAACCGCTTCATGAGGAGCACCCGCAT‑AACCTGCACCAGCCGTACCGTATTCCAGGGGAGGAGCATGATCCGGATTCGGGGCGGTACTATAACCGTCACAGGACGTAGGATCAGTTGCCGGGGTGGtatagtactcagga                                                                                                     >  2:246592/1‑290 (MQ=14)                                                                                                                                                                                                             tGGAGCCGGAATAGACACCGGCGCGAAAAGGTCATGTGGACCACTGCGACCACCGGGGCCGGCCGTGGGCGGTGAGCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGATGATGAATGGGGCAACCAGCTTAATTAGGAGAAACCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATc                                                                  <  2:297956/300‑1 (MQ=17)                                                                                                                                                                                                               gAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCAACGGGGCCTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGGGGAGTATGATGAATGGGGGAACCAGGTTAATGAGGAGAACCAGCAG‑AGCATACAGCAGCCGAACAGGCAGCCAGGGCAGCAGCGTGATACGGCGTCGGGACGGTAATATAACCGTCCCCGGCACTACTGACGGTGGCGGGGGCGGTCGAACAGCCCTGGCCCAGGTGGGTTGGGAGGgggcgtgcagatgtgtagga                                                               >  2:385528/1‑284 (MQ=9)                                                                                                                                                                                                                    ggAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCGTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGAAGTATGATGAATGGGGCAACCAGCTTAACGAGGAGAAACGCCAT‑CACCTGAAGCAGCCGTACCGTGTGCCAGGGCAGCAGACTGATGAGGAGTCGGGGCTATACTATACCCGTCACCGGCACTACGATCCGTTGCACGGGCGGTATATGACCACGGACCCGCTTGGGTTggcagtaggaagggagatgaaacagttgc                                                            >  2:331228/1‑271 (MQ=255)                                                                                                                                                                                                                            cacCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCCTGATGAGGAGTCGGGGCTGTACTATAACCCTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCCCCCCGAGCCCGGTTGGGTTGAGAGGGGGGTGGAAGATGTATCAGTATCCGtttgatc                                                   >  1:297704/1‑295 (MQ=17)                                                                                                                                                                                                                                  agggaaaattgaatgtttACCCCGGGGACCCCCGGGGCCGTCCGCTGGCGCTTAGCAGCGAAGACGGCAATAAGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCATCTTACGAAGGAGAACCCGGAG‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATAAGGAGTGGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATAc                                             <  1:331220/285‑1 (MQ=255)                                                                                                                                                                                                                                                            gACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATGGGGTTGAGCGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATCCAAGTGATGGACCCAATGGGGTTAGATg                  >  1:322090/1‑301 (MQ=255)                                                                                                                                                                                                                                                                            cGCGGGCGCTAAGCAGCGTAGACGTCCATACGGCGTGGCTCGGGTGTTAGGACGAATGGGGCACCCAGCTTAATGCGGAGAACCCGCAT‑CACCTGCAACAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATATCCGGCACCGGCACTACGATCCGTTGCAGGGGCTGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGAc  <  2:297710/301‑1 (MQ=21)                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    CGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCA‑GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAG‑ACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAT‑CACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGAC  >  NC_000913/1528943‑1529508

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