Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,218	C→T	pseudogene (1317/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,218	0	C	T	97.8%	102.5 / NA	46	pseudogene (1317/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/21);  minor base G (1/0);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGG  >  NC_000913/1528977‑1529514                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         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)                                                                                                                                                                                                 agaaaaacccctctttcgcACTGCCACCACCGGGGCCGTCCGCTGCCGCTTATCAGCGAAGCCGGCAATCAGGCGTGGGGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCGGAACCAGCCGGACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAg                                                <  2:361646/283‑1 (MQ=17)                                                                                                                                                                                                       cctcttttttACCCCTGCGACCACCGGGGACTGCCGGTGGCGCTTTTCAGCGAAGGCGGCAATAGGGGGGGGGGGGGGGAGTATGAGGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCAGGGCAGAAGCAGGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAg                                          <  2:460931/296‑1 (MQ=17)                                                                                                                                                                                                        cggatcttgagtccTGCGACCACGGGGGACGGCCGGTGGAGCTTAACAGGGGAGGCGGCAGTACGGGGTAGCGGGGGGAGTAAGACGAATGGGGCAACGAGGGTAATGAGGAGAACCCGAATCACCTGCACCAGCCGTACCGTCTGCCGGGGCAGCAGAATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCCGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGa                                         <  2:131437/288‑1 (MQ=255)                                                                                                                                                                                                          cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCCGCTTAATGAGGAGAACACGCATCACCTGCACCAACCGTCACGTATGCCAGGGAGGCAGCCTGATGAGGAGTCGGGACTGTACGATAACCGTAACCCGCCATACGATCCGATGCCGGGGCGGTCGACCACCCCGGAACCGTGTGGGTTGGGAGGTTGATGGAGTATGTGACATTAACCGTTCACTCCCATTCAAGTGGTTGAcc                                       >  2:347327/1‑301 (MQ=14)                                                                                                                                                                                                                 atcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAACGTCACCGGCACTACGCCCCGTTCCAGGGGCGGTAGATCACCCCAGACCCGCGTGGGGGGAGAGGAGGATGGACTATGTATTAGTAGCCGTGGAAACCGCTAACAGTGAAAGACACAACggg                                >  2:464738/3‑301 (MQ=17)                                                                                                                                                                                                                  tcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                                                      >  1:76181/2‑166 (MQ=255)                                                                                                                                                                                                                  tcACTGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                                                      <  2:76166/165‑1 (MQ=255)                                                                                                                                                                                                                       gCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCGGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGTTATTTCAACCCGGACCCGACTGGGGTGGGGGGGGGATGGTATATGTCTCAGTACACGGTGAAGCCCAAAACAGTGaatggacctatgggggtgg                           >  2:163367/1‑282 (MQ=14)                                                                                                                                                                                                                                      cgCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGTCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAatat           <  1:287812/300‑1 (MQ=255)                                                                                                                                                                                                                                               cgcTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCGGTGGCAGGGGCGGTATATGACCACGGACCCGATTGGGTGGAGCGGTGGATGGAATGTGTATCCGGGTACGTTGAATGCCATACCAGTGAGAGACCCACAGTGGTGAGATGCGAGTGGTAATATGACATAAgg  >  2:92422/1‑301 (MQ=17)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         CGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGG  >  NC_000913/1528977‑1529514

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