Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,203	C→A	pseudogene (1302/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,203	0	C	A	95.5%	96.3 / NA	44	pseudogene (1302/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (23/19);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAG  >  NC_000913/1528977‑1529489                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                               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)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                               CGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACA‑CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAG  >  NC_000913/1528977‑1529489

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