Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,574	C→A	D84E (GAC→GAA)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,574	0	C	A	100.0%	160.2 / NA	49	D84E (GAC→GAA)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/25);  total (24/25)

CTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTC  >  NC_000913/4617327‑4617817                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    atGATCTAAAACCAAGCAGCCTGTGTGCACTGAAATGGA‑GGACCCTCACCCCCCTGAATGACGACGACCCCGCCGAAAAAGTAATTGCT‑CTGTTTCAGCAGGCCAAACCCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATgacatcgaaatc                                                                                                                                                                                                                                                  <  1:285173/250‑1 (MQ=255)        gAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGCCC‑TGCCCCCCCGAAATGACGCCAACCCCGACGAAAAAGTAATTGCT‑CTGTTTCACCAGACAAAACCCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATcgc                                                                                                                                                                                                                                                <  1:259405/246‑1 (MQ=255)                     cTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGc                                                                                                                                                                                                                              >  1:139142/1‑251 (MQ=255)                       gCGTGCACTAAAATTTATGGACC‑TGCCCCCCCTGAATGCCGACGACCCCGACGATAAATTAATTTCT‑CTGTGTCTTTGGGCCAAACCCCCGGTC‑GCCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACTATGACATCGAAATCGCGCTGGCAGAAACCCGTGCgg                                                                                                                                                                                                                            <  1:396346/251‑1 (MQ=255)                         gTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CCGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                                                                                                                                                                                                           >  2:303330/1‑250 (MQ=255)                          tGCACTGAGATAGCTGGACC‑TGACCACCCTGAATTCCGACGACACCGACGAGAACGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAAAACCCTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGCCAGAAACCCGTGCGGcc                                                                                                                                                                                                                          >  2:267937/1‑249 (MQ=255)                                     tGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACggg                                                                                                                                                                                                              >  1:252660/1‑250 (MQ=255)                                                     cccTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTACGGCAATCGCCTACGGTGCCGATGAAGTTGAc                                                                                                                                                                                               >  1:341694/1‑251 (MQ=255)                                                       cTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCACAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCACTCGCCTCCGGTGCCGATGAAGTTGACGt                                                                                                                                                                                             >  1:301789/1‑251 (MQ=255)                                                           aggacgacgacACCGACGAAAAAGTAAT‑GCCTCTGTGTCATC‑GGGCAAAACCCCCGTCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACgtggtg                                                                                                                                                                                         <  1:463292/249‑1 (MQ=255)                                                           aTGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAAccc                                                                                                                                                                                                                                  >  2:6805/1‑210 (MQ=255)                                                               cgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTcc                                                                                                                                                                                     >  1:236254/1‑251 (MQ=255)                                                               cgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTcc                                                                                                                                                                                     >  1:121559/1‑251 (MQ=255)                                                                             gagaaagtatgtgctttttgtCATCACGCAAAACCCCCGGCC‑GACAACACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCt                                                                                                                                                                        <  2:160787/233‑1 (MQ=255)                                                                                     aattgctcTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCgg                                                                                                                                                                >  1:83253/8‑251 (MQ=255)                                                                                                 tcatcaGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGACGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAagagacgg                                                                                                                                                    >  2:254690/1‑247 (MQ=255)                                                                                                   atcaGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggtg                                                                                                                                                  >  2:25387/1‑250 (MQ=255)                                                                                                   acaaggccaaatccgcGGTT‑GGCAATACCGCCGCTAACTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACTCCGGAAGTCCGTATGGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACGCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTTGTGTGCCCGTTCTGCGCGCTGATGGCTGGTAACGAGCAggtt                                                                                                                                                  <  2:280939/236‑1 (MQ=255)                                                                                                         ccaaacccccGGTC‑GGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGtttt                                                                                                                                            <  1:20784/243‑1 (MQ=255)                                                                                                           aaacccccGGTG‑GGCAACTCCGCCGCTATCTGTATCTATCCTCGCTTTAGCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGa                                                                                                                                          <  2:139142/245‑1 (MQ=255)                                                                                                              aaccccgtgaggaaaaaCCGCCGCTACCTGTATCTATCCGCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTTCCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAcc                                                                                                                                        <  1:267937/235‑1 (MQ=255)                                                                                                                   gggcgggcactgccgccgGTATCTTTATCTATCCTCGCGTGATCCCGAGTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAGTCCGTATTGCTACGGTAACCAACTGCCCTCACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTg                                                                                                                                   <  2:291173/239‑1 (MQ=255)                                                                                                                      c‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAg                                                                                                                               <  1:214695/251‑1 (MQ=255)                                                                                                                                         cTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGGGCAGGGTGGGTTTTACCTGGTGAAAGCCTGGAAAAAGGCTTggg                                                                                                             >  1:129367/1‑249 (MQ=255)                                                                                                                                             atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccg                                                                                                                                                                                                          <  1:324549/154‑3 (MQ=255)                                                                                                                                             atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccg                                                                                                                                                                                                          >  2:324549/1‑152 (MQ=255)                                                                                                                                             atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCgg                                                                                                          <  1:178662/250‑1 (MQ=255)                                                                                                                                               ctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGCCCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGccg                                                                                                       >  1:437563/1‑249 (MQ=255)                                                                                                                                                           tATCCCGATTGCTCGCCAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTActtc                                                                                           >  1:412365/1‑249 (MQ=255)                                                                                                                                                                cGATTTCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAAGCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggt                                                                                                                                                   <  2:423539/190‑1 (MQ=255)                                                                                                                                                                cGATTGCTCGCAAAACCCTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGGGCAgga                                                                                                                                                   >  1:423539/1‑189 (MQ=255)                                                                                                                                                                     gCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGat                                                                                 <  1:243380/251‑1 (MQ=255)                                                                                                                                                                         gCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGGACTGCTGAAAGTGatcatc                                                                             >  1:263191/1‑251 (MQ=255)                                                                                                                                                                          cAAAACACGGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCGACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCg                                                                            <  2:263191/251‑1 (MQ=255)                                                                                                                                                                                           aGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGaaa                                                                         <  1:448362/237‑1 (MQ=255)                                                                                                                                                                                           aGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCCGCGGATGTACTGCTGAAAGTGATCATCGAAAACGGCGAACAGaa                                                            >  1:329459/1‑250 (MQ=255)                                                                                                                                                                                                 ccccGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGaa                                                     <  1:16768/251‑1 (MQ=255)                                                                                                                                                                                                   ccGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAgc                                                   <  1:267087/251‑1 (MQ=255)                                                                                                                                                                                                       aaaTCCTTATTGCTACGGTAACCAACTGCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTg                                               <  1:242681/251‑1 (MQ=255)                                                                                                                                                                                                         aTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGAt                                             >  1:312813/1‑251 (MQ=255)                                                                                                                                                                                                         aTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGGGGCCGGGAATGGACTGCTGAAAGGGATCATCGAAACCGGCGGACTGAAAGAAGAAGCGCTGa                                              >  1:265984/1‑250 (MQ=255)                                                                                                                                                                                                           ccgtagtgCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATcc                                           <  2:412365/244‑1 (MQ=255)                                                                                                                                                                                                           ccGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATcc                                           <  1:303330/251‑1 (MQ=255)                                                                                                                                                                                                               attgCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAAGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCCGCGAATGTACTGCTGGAAGTGATACTCGAAACCGGCGAACTGAAAGACGAAACGCTGAACCGTa                                        >  1:296078/4‑250 (MQ=255)                                                                                                                                                                                                                 tgCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTAGGGTGCCGATGACGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGGACAGAGGCTTGCGCGGGAGCGAATGTACAGCTGAAAGGGATTAACGAAACCGGCGAGCTGCAAGACGAAgggccgcgccgtcaa                                      >  1:207439/2‑236 (MQ=255)                                                                                                                                                                                                                                  tCCCACACGGTAACGATGACAGCGAAAACGCGATGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCa                      <  2:171144/249‑1 (MQ=255)                                                                                                                                                                                                                                   cccccacGGTAACGATTACATCGAAATCGCGCTGGCAGAACTCCGTGCGGCAATCGCCTACGGTGCCGATGAAGGTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATc                    <  2:124532/246‑1 (MQ=255)                                                                                                                                                                                                                                           gTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGt           <  1:451368/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGGACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCgggg           >  1:375251/1‑250 (MQ=255)                                                                                                                                                                                                                                             aaCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGc         <  1:481760/251‑1 (MQ=255)                                                                                                                                                                                                                                                    acatcgatatcgCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGTCGGCGGTTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGGAAAGATGCTTGCGCGGCAGCGGATGTACGGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGTCGACGCGCTGAGGCGTAAAGTGTCTGAAATCACCATCAAAGCGGGTGCGGACTtc  <  1:418260/251‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    CTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTC  >  NC_000913/4617327‑4617817

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