Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,175	T→C	pseudogene (1274/2037 nt)	rhsE →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,529,175	0	T	C	100.0%	82.1 / NA	36	pseudogene (1274/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/16);  total (20/16)

CGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGT  >  NC_000913/1528977‑1529447                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                 cGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTCCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCCGACCGCGTGAGCAGTGACAGCCGGGCGTGGCTTGCGCAGTTCGGGCTGACGGTGGAGCAACTGCCCAGCCACGTGGAGCCGGCATACACC‑CCGGCGCGAAAAGCCCAACTTTATCACTGCCAACAACGGGTACTGCCGCTGACGCCGCTCAGCGAATACGCCACTACGGCGGGGCggggggcgattgcgtaattgggcaaacggcg                                                                                                                                                                            >  2:239934/1‑276 (MQ=11)     ccAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGATGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTCGCATGCGCAGTGCGGGCTGGCGGAGGAGCAACTGGCCCGACCGATGGAACCGGAATACACA‑CCGGCGCGACAAGCTAATCTTTATCACTGCGACCACCGGGGACTGACTATGGCGCTTATAAGCGAAAACGGCAATACGCCAGGGCGCGGGCACTCTGATACATggggtcagcagccaaa                                                                                                                                                                         >  1:394383/1‑288 (MQ=11)                            gcgaggacGTGGTGTGCCAGGCTGAACTGGTGCGGCTGCTGGACAGTGAGGGGGAAGAAATCCTGGCAGACCGCGTGAGCGGTGAAAGCCGGGCGTGGCTTGCGCAGGGCGGGCTGACGGTGGCGGAAATGGCCAGACAGGTGGAGCCGGAATACACA‑CAGGCGCGAAAAGCGCGTCTTTATCACTGCGAACACCGGGGACTGCCGCTGGCTCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGTACCCGCATCACCTGCa                                                                                                                                                 <  1:57334/294‑1 (MQ=255)                                   cGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCg                                                                                                                                          >  1:131468/1‑301 (MQ=14)                                                             ggCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCCCATCTTTATCAATGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAACCGCCCATACGGCGTGGCGCGGGGAGTATGATGAACGGGGCAACCAGCTCCAAGAGAAGAACCCGCACGACCTGCACAAGCCGTACCTTCAGCCAGGggcgcagcacg                                                                                                                 >  2:332684/1‑298 (MQ=11)                                                               gagctgGAAAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGGAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgatg                                                                                                              <  1:483070/299‑1 (MQ=9)                                                                 gctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTCTCACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAAGACGGCGGGGCGCGGGGGGTATGATGAAGGGGGCAAACAGCTTACTGCGGAGAACGCGCATCACCCGCACCAGCCCTCACGCGTGCCAGGgcaaacgcatcagt                                                                                                              >  2:482953/1‑288 (MQ=11)                                                                 gctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgatg                                                                                                              >  1:73076/1‑299 (MQ=14)                                                                     gACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCAGGATgaggagg                                                                                                        >  1:351050/1‑300 (MQ=12)                                                                        aGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCACCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTggg                                                                                                     >  1:283072/1‑301 (MQ=12)                                                                         ggCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGCGCGGGGGGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCCCCAGCCGTACCGCCTGCCCGGGCGGCCACCAGACGAGGAGTCggg                                                                                                    >  2:526160/1‑301 (MQ=11)                                                                           ggggagaaggaagcccGGGCAGACCGGGGGGGCGGGGAAAGCCGGGGGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                  <  2:350961/287‑1 (MQ=11)                                                                             ggaggaAGAAATCCGGGCAGACCGCGTGAGCAGTGACAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTg                                                                                                <  1:459897/301‑1 (MQ=255)                                                                                 gaagaaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACt                                                                                            <  1:394217/301‑1 (MQ=18)                                                                                             gCAGACCGCGTGAGAAGTGAAAGCCGGGCGTGGCTGGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCAcc                                                                                <  1:332770/301‑1 (MQ=255)                                                                                             ccgtttcccagcgggagggAAAGCCGGGCGGGGGTTGCGCAGTGCTGGCTGACGGTGGAGCAACGGGCCAGCCAGTTGGAGCCGGAATACACC‑CCGGCGCGAAAAGCTAATCGGTATCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAGGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCAcc                                                                                <  2:214833/283‑1 (MQ=11)                                                                                               aGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCCCCAGCCGTACCGTCTGCCAggg                                                                                                                           >  1:554499/1‑256 (MQ=12)                                                                                                          gCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGGGTATGATGAATGGGGAAACCAGCTTAATGAGGAGCACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCCGCCGCATGATGAGGAGGTGGGGCCGTACTCTTAGCGGCAccagataggaagcgc                                                                   >  1:214892/1‑288 (MQ=11)                                                                                                               gAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCCCTACGATCCGTTGc                                                              >  1:313212/1‑301 (MQ=255)                                                                                                               gAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACgggatcgga                                                               >  1:559625/1‑292 (MQ=255)                                                                                                               gAAAGCCGGGCGTGGCTGGCGGCGGGCGGGCGGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACCCA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACg                                                                       <  2:559493/292‑1 (MQ=11)                                                                                                                   ggagggCGGGGATTGCGCAGTGCGGGCTGACGGTGGAGAAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGGGAAAAGCTCATCTTTATCACTGAGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAggg                                                          <  1:239999/298‑1 (MQ=11)                                                                                                                      gggCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACGGTCACCGGCACTCCGATCCGTTGCAGGGGCg                                                       >  1:452930/1‑301 (MQ=255)                                                                                                                              ggtGCGCAGTGGCGGCCGACGGTGGAGGAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCa                                               <  1:141079/299‑1 (MQ=255)                                                                                                                                        gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGCTGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGAccc                                     >  1:206753/1‑301 (MQ=255)                                                                                                                                          gggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCgg                                   >  1:565989/1‑300 (MQ=255)                                                                                                                                               gACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTggg                              <  1:526286/301‑1 (MQ=255)                                                                                                                                                   gTGGAGCAACTGGCCAGACAGGTGGGACCGGAATACACA‑CCGGCGCGAAAAGCGCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTga                          <  1:86910/301‑1 (MQ=255)                                                                                                                                                         cAACTGGCCAGACAGGTGGGGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag                                                                                                         <  2:86181/216‑1 (MQ=255)                                                                                                                                                         cAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag                                                                                                         >  1:86196/1‑216 (MQ=16)                                                                                                                                                              ggCCAGACAGGTGGGGCCGGAAGACCCC‑CCGGCGGGAAAAGCTCATCTTTAACACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtggat                 <  2:73063/299‑1 (MQ=255)                                                                                                                                                              ggCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTAATATAACCGTCACCGGCACTACGACCCGTTGCAGGGGGGGGATATCACCCCGGCCCCGATTGGGTTGTGAGGtgtaggg               >  2:223269/1‑296 (MQ=14)                                                                                                                                                               gCCAGACAGGTGGAGCCGGAAT‑CACACCCAGGGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACGGCCGCTGGGGCGTATCAGGGAAGACTGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtggatgga              <  2:452819/301‑1 (MQ=14)                                                                                                                                                                 cagcaaGGGGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTATTCTTTATCACTGCGACCACGGGGGACGGCCGCTGGCGCTTATCAGCGAGGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAat            <  2:565855/296‑1 (MQ=14)                                                                                                                                                                       ggTGGAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATAACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGAGGAGTGTTCGGGGCTGTACTATAAACGTCACCGGCACCACGATCCGTTGCAGGGGCGGTATGTCACACCGGACCCGATTGGGGTGAGAGGtggggggattggcttt      >  2:287729/1‑288 (MQ=14)                                                                                                                                                                           gAGCCGGAATACACA‑CCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTTGGTTGAGAGGTGGATGGAATATGTATCAgg  >  1:528965/1‑300 (MQ=255)                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                 CGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGT  >  NC_000913/1528977‑1529447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑