Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,185	C→T	pseudogene (1284/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,185	0	C	T	97.3%	78.0 / NA	38	pseudogene (1284/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/15);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACC  >  NC_000913/1528977‑1529477                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                     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)                                                                                                                                                                                                       cctcttttttACCCCTGCGACCACCGGGGACTGCCGGTGGCGCTTTTCAGCGAAGGCGGCAATAGGGGGGGGGGGGGGGAGTATGAGGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCAGGGCAGAAGCAGGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAg     <  2:460931/296‑1 (MQ=17)                                                                                                                                                                                                          cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCCGCTTAATGAGGAGAACACGCATCACCTGCACCAACCGTCACGTATGCCAGGGAGGCAGCCTGATGAGGAGTCGGGACTGTACGATAACCGTAACCCGCCATACGATCCGATGCCGGGGCGGTCGACCACCCCGGAACCGTGTGGGTTGGGAGGTTGATGGAGTATGTGACATTAACCGTTCACTCCCATTCAAGTGGTTGAcc  >  2:347327/1‑301 (MQ=14)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                     CGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACC  >  NC_000913/1528977‑1529477

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