Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,629,247	C→T	I11I (ATC→ATT)	ydfZ →	putative selenoprotein YdfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,629,247	0	C	T	100.0%	178.6 / NA	51	I11I (ATC→ATT)	ydfZ	putative selenoprotein YdfZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/28);  total (23/28)

ACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCC‑AGATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTA  >  NC_000913/1628951‑1629534                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                 aCGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCC‑AGATCCGCCAGGAAAACA‑GCGACTGG‑TGTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAACTTGATCAAAAACAAAAAATATTCCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGTAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATtacc                                                                                                                                                                                                                                                                                               <  2:448887/300‑1 (MQ=255)   gTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCC‑AGATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATtaccac                                                                                                                                                                                                                                                                                             >  1:67315/1‑300 (MQ=255)                  aTTTGCAGGAAATCAGCGAATCCGTGCGCC‑AGATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAATAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTa                                                                                                                                                                                                                                                                              <  2:152757/299‑1 (MQ=255)                  aTTTGCAGGAAATCAGCGAATCCGTGCGCC‑AGATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTa                                                                                                                                                                                                                                                                              <  2:152887/299‑1 (MQ=255)                    ttGCAGGAAATCAGCGAATCCGTGCGCC‑AGATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATgg                                                                                                                                                                                                                                                                           <  2:451643/300‑1 (MQ=255)                    ttGCAGGAAATCAGCGAATCCGTGCGCC‑AGATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATgg                                                                                                                                                                                                                                                                           <  1:94030/300‑1 (MQ=255)                              tCAGCGAATCCGTGCGCC‑AGATCCGCCATGAAAACA‑GCTACTGG‑TGTAGCGAAGGGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACATGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCAc                                                                                                                                                                                                                                                                 <  2:202048/300‑1 (MQ=255)                               agccGAATCCGTGCGCCNAGATCCGCCAGGAAAACANGCGACTGGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCNGTGTTATGGTTAGCgg                                                                                                                                                                                                                                                                    <  1:564369/297‑1 (MQ=255)                                         gTGCGCC‑AGATCCGCCATGAAAACA‑GTGACTGG‑TTTATCGAAGAGTAATTCATTTCCTCCCATCCCATCCGGTGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGGGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCacac                                                                                                                                                                                                                                                        <  2:141165/298‑1 (MQ=255)                                         gTGCGCC‑AGATCCGCCATGAAAACA‑GCGACTGG‑TTTAGTGAAGAGTGATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGGTGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGATTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCTGCACCGGTCacac                                                                                                                                                                                                                                                        <  2:140977/298‑1 (MQ=255)                                          tGCGCC‑AGACCCGCTAGGAAAACA‑GCGACTGG‑TTTAGCTAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGTAGAAAATACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCacac                                                                                                                                                                                                                                                        <  2:140423/297‑1 (MQ=255)                                                  gATCCGCCATGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATc                                                                                                                                                                                                                                              <  1:253334/300‑1 (MQ=255)                                                  gATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATNGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATc                                                                                                                                                                                                                                              <  2:533069/300‑1 (MQ=255)                                                  gATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATc                                                                                                                                                                                                                                              <  1:522575/300‑1 (MQ=255)                                                                         tGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTc                                                                                                                                                                                                                        <  2:621597/300‑1 (MQ=255)                                                                           g‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCt                                                                                                                                                                                                                       >  1:170395/1‑299 (MQ=255)                                                                                                   ccTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGNCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTAGATTGATACTGAAGGTCTGACCGCTGAGCAAATCcgccgc                                                                                                                                                                                                 >  2:483420/1‑298 (MQ=255)                                                                                                    ctctCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGa                                                                                                                                                                                              >  1:621286/1‑300 (MQ=255)                                                                                                            ccatccGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGAACTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGAAAAAAcagt                                                                                                                                                                                      >  2:181178/1‑297 (MQ=255)                                                                                                            ccatccGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGAAAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCAATGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAAACGt                                                                                                                                                                                      >  2:181034/1‑300 (MQ=255)                                                                                                             catccGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCA‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTa                                                                                                                                                                                     >  2:53507/1‑300 (MQ=255)                                                                                                                         gagTCTTTTCCCCCTCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGt                                                                                                                                                                         <  2:245540/300‑1 (MQ=255)                                                                                                                         gagTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGt                                                                                                                                                                         <  2:622313/300‑1 (MQ=255)                                                                                                                                ttCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAag                                                                                                                                                                  <  1:637480/300‑1 (MQ=255)                                                                                                                                 tCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAaga                                                                                                                                                                 >  2:345026/1‑300 (MQ=255)                                                                                                                                       gCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCCAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACt                                                                                                                                                           >  1:454365/1‑300 (MQ=255)                                                                                                                                       gCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACt                                                                                                                                                           >  1:454371/1‑300 (MQ=255)                                                                                                                                        ccTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACAGGTCACACTGGCAAGATCTTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACt                                                                                                                                                           >  2:168067/1‑299 (MQ=255)                                                                                                                                               gCTCATGCATGCATCAAAAACGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGAAAGGAGTATTTATATTGCGCCCGTTAAGTGGGCTTCGAATGTAAATCAGAAAGGAGAAAACACCTATGACGAACTACGATCGTAACCGTAACGCAATTACCAATGGCAGAC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAAGTCTGACCGCTGAGCAAATCCGCCGAGGAACAAACGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGcacaaa                                                                                                                                                   >  2:382959/1‑297 (MQ=255)                                                                                                                                               gCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGcaccac                                                                                                                                                   >  2:412074/1‑300 (MQ=255)                                                                                                                                               gCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGCCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGcaccac                                                                                                                                                   >  1:210450/1‑300 (MQ=255)                                                                                                                                               gCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAAACGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGcaccac                                                                                                                                                   >  2:616725/1‑300 (MQ=255)                                                                                                                                                cTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTANGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACt                                                                                                                                                  >  2:485394/1‑300 (MQ=255)                                                                                                                                                cTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACt                                                                                                                                                  >  2:403940/1‑300 (MQ=255)                                                                                                                                                     gcatgcatCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGAcc                                                                                                                                             <  1:460778/300‑1 (MQ=255)                                                                                                                                                     gcatgcagCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATAGTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGAcc                                                                                                                                             <  2:473808/300‑1 (MQ=255)                                                                                                                                                                              tGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGc                                                                                                                    <  2:586381/300‑1 (MQ=255)                                                                                                                                                                               gATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACTTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGAACTACGATCGTAACCGTAACGCAATTACCACTTGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGGGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCg                                                                                                                   <  2:189464/300‑1 (MQ=255)                                                                                                                                                                                   aaaaacaaaaaATATTCCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAACCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCg                                                                                                                   <  2:95720/296‑1 (MQ=255)                                                                                                                                                                                                             ggAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCAGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTAGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttt                                                                                     >  2:439712/1‑300 (MQ=255)                                                                                                                                                                                                                gTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTAGGCATGAATTAAGCGTGTGAATGCCGCCGATGGCGGCATTGTtttttt                                                                                  >  2:48978/1‑300 (MQ=255)                                                                                                                                                                                                                  aTTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTAc                                                                                <  1:53507/300‑1 (MQ=255)                                                                                                                                                                                                                                      gTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                             <  1:108958/299‑1 (MQ=255)                                                                                                                                                                                                                                      gTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATGGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAAAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGGATATTTTGCCa                                                             <  1:108945/299‑1 (MQ=255)                                                                                                                                                                                                                                      gTGGGCTGCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                             <  1:108782/299‑1 (MQ=255)                                                                                                                                                                                                                                                                       cacCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTAGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACaa                           >  1:160354/1‑300 (MQ=255)                                                                                                                                                                                                                                                                           tATGACGACCTACGATCNTAACCGTAACGCAATTACNACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGATCTGAGTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACTGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGNTAATAAAGACAAGTAc                       <  1:576208/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                  acCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                 >  1:238316/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                  acCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTAGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                 >  1:449834/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                     tacaaTCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCa             <  2:67315/296‑1 (MQ=255)                                                                                                                                                                                                                                                                                     tACGATCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGCCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTNGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGt               <  2:210450/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                          tCGTAACCGTAACGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTAGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCtt        >  2:2946/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                ccgtaacGCAATTACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTa  >  1:556520/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                 ACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCC‑AGATCCGCCAGGAAAACA‑GCGACTGG‑TTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCC‑GTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTA  >  NC_000913/1628951‑1629534

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑