Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,616	T→C	A98A (GCT→GCC)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,616	0	T	C	100.0%	169.6 / NA	51	A98A (GCT→GCC)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/28);  total (23/28)

CCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGA  >  NC_000913/4617378‑4617861                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         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)           cgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTcc                                                                                                                                                                                                                                    >  1:121559/1‑251 (MQ=255)           cgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTcc                                                                                                                                                                                                                                    >  1:236254/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)                                                                                         atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCgg                                                                                                                                                         <  1:178662/250‑1 (MQ=255)                                                                                           ctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGCCCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGccg                                                                                                                                                      >  1:437563/1‑249 (MQ=255)                                                                                                       tATCCCGATTGCTCGCCAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTActtc                                                                                                                                         >  1:412365/1‑249 (MQ=255)                                                                                                            cGATTTCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAAGCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggt                                                                                                                                                                                                  <  2:423539/190‑1 (MQ=255)                                                                                                            cGATTGCTCGCAAAACCCTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGGGCAgga                                                                                                                                                                                                  >  1:423539/1‑189 (MQ=255)                                                                                                                 gCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGat                                                                                                                              <  1:243380/251‑1 (MQ=255)                                                                                                                     gCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGGACTGCTG‑AAAGTGatcatc                                                                                                                          >  1:263191/1‑251 (MQ=255)                                                                                                                      cAAAACACGGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCGACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCg                                                                                                                         <  2:263191/251‑1 (MQ=255)                                                                                                                                       aGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGaaa                                                                                                                      <  1:448362/237‑1 (MQ=255)                                                                                                                                       aGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCCGCGGA‑TGTACTGCTG‑AAAGTGATCATCGAAAACGGCGAACAGaa                                                                                                         >  1:329459/1‑250 (MQ=255)                                                                                                                                             ccccGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGaa                                                                                                 <  1:16768/251‑1 (MQ=255)                                                                                                                                               ccGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAgc                                                                                               <  1:267087/251‑1 (MQ=255)                                                                                                                                                   aaaTCCTTATTGCTACGGTAACCAACTGCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTg                                                                                           <  1:242681/251‑1 (MQ=255)                                                                                                                                                     aTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGAt                                                                                         >  1:312813/1‑251 (MQ=255)                                                                                                                                                     aTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGGGGCCGGGAA‑TGGACTGCTG‑AAAGGGATCATCGAAACCGGCGGACTGAAA‑GAAGAAGCGCTGa                                                                                          >  1:265984/1‑250 (MQ=255)                                                                                                                                                       ccgtagtgCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATcc                                                                                       <  2:412365/244‑1 (MQ=255)                                                                                                                                                       ccGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATcc                                                                                       <  1:303330/251‑1 (MQ=255)                                                                                                                                                           attgCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAAGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCCGCGAA‑TGTACTGCTG‑GAAGTGATACTCGAAACCGGCGAACTGAAA‑GACGAAACGCTGAACCGTa                                                                                    >  1:296078/4‑250 (MQ=255)                                                                                                                                                             tgCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTAGGGTGCCGATGACGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGGACAGAGGCTTGCGCGGGAGCGAA‑TGTACAGCTG‑AAAGGGATTAACGAAACCGGCGAGCTGCAA‑GACGAAgggccgcgccgtcaa                                                                                  >  1:207439/2‑236 (MQ=255)                                                                                                                                                                              tCCCACACGGTAACGATGACAGCGAAAACGCGATGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCa                                                                  <  2:171144/249‑1 (MQ=255)                                                                                                                                                                               cccccacGGTAACGATTACATCGAAATCGCGCTGGCAGAACTCCGTGCGGCAATCGCCTACGGTGCCGATGAAGGTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATc                                                                <  2:124532/246‑1 (MQ=255)                                                                                                                                                                                       gTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGt                                                       <  1:451368/251‑1 (MQ=255)                                                                                                                                                                                       gTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGGACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCgggg                                                       >  1:375251/1‑250 (MQ=255)                                                                                                                                                                                         aaCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGc                                                     <  1:481760/251‑1 (MQ=255)                                                                                                                                                                                                acatcgatatcgCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGTCGGCGGTTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGGAAAGATGCTTGCGCGGCAGCGGA‑TGTACGGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GTCGACGCGCTGAGGCGTAAAGTGTCTGAAATCACCATCAAAGCGGGTGCGGACTtc                                              <  1:418260/251‑1 (MQ=255)                                                                                                                                                                                                    cgaaatcgCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTtcatca                                          <  1:25387/247‑1 (MQ=255)                                                                                                                                                                                                             gcTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTa                                 <  1:215748/251‑1 (MQ=255)                                                                                                                                                                                                                 gCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCgg                             <  1:10781/251‑1 (MQ=255)                                                                                                                                                                                                                   aGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGGGAAATG‑ACTTCTGGAAAG‑GATCATCGAAAACGGCGAAATGAAAAGA‑GAAACGCTGGTCCGGAAAGCGGCCGGAAACTCCACCAAAACGGGGGCGGACTTTATCCAAAcccccacccggt                           >  2:85379/1‑241 (MQ=255)                                                                                                                                                                                                                    gAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGGGATCATCGAAACCGGCGAACTGAAA‑GACGAAGAGCTGATCCGTAAAGCGTCTTAAATCTCCATCAAAGCGGGGGCGGACTTCCTCAAAACCTCTACCGGTa                           >  1:93269/1‑250 (MQ=255)                                                                                                                                                                                                                        cccGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGGA‑TGTACTGCTG‑AAAGTGGTCATCGAAACCGGCGGACTGAAA‑GACGCAGCGCTGATTCGTAAAACGTCTGCAATCTCCCTCAAAGAGGGTGCGGACTTCCTCCAAACCTCTACCGGGaaaatg                      >  2:213898/1‑248 (MQ=255)                                                                                                                                                                                                                                 cAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAAc             <  1:254690/251‑1 (MQ=255)                                                                                                                                                                                                                                       ccTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGa               <  2:108365/243‑1 (MQ=255)                                                                                                                                                                                                                                       ccTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGGGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGaagatcggg       >  1:108365/1‑244 (MQ=255)                                                                                                                                                                                                                                            ggTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCCGCGCA‑TGTACTGCTG‑AAAGGGATCATCGAAACCGGCGGACTGAAA‑GACGCAGCGCTGATCCGCAAAGCGTCTAAAATCTCCATCAAAGCGGGTGCGGACCTCATCAAAACCTCTACCGGTAAAGCGGGTGTGAAAGCGAAGCCggg  >  1:277692/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         CCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAA‑TGTACTGCTG‑AAAGTGATCATCGAAACCGGCGAACTGAAA‑GACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGA  >  NC_000913/4617378‑4617861

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