Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,532	C→T	I70I (ATC→ATT)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,532	0	C	T	100.0%	238.7 / NA	73	I70I (ATC→ATT)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (43/30);  total (43/30)

CAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAA‑AACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGT  >  NC_000913/4617285‑4617778                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        cAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCATGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCt                                                                                                                                                                                                                                                      >  1:222799/1‑251 (MQ=255)  aaGGTGAATTTTATGGGCGACAAGCAAGGAGAATGAAATGAATGATCTGAAAGCAAGACGCCGGCGGGCGCTGAAATGGATGGCCCTGACCCCCCTGAAGGACAACGACACCGACGAGAAAGTAAT‑GCCTCGGGGTCATCAGGCCAA‑ACCCCGGGCCGGAAATACCGCCGCTATCTGTATCGATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTa                                                                                                                                                                                                                                                     <  1:614593/251‑1 (MQ=255)    ggTGAATTTTATTGGCGACAAGCCAGGAGATTGAAAGGACGGATCGAAAAGAAAGCAGCCTGCGTGCACTGAATGTGATGGACCTGACCCCCCTGAATGACGACGACCCCGACGAGAAAGTAATTGCT‑CTGTGTCTTAAGGCAAA‑AACCCCGTTCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACg                                                                                                                                                                                                                                                   <  1:332094/251‑1 (MQ=255)         aTTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAc                                                                                                                                                                                                                                              >  1:494934/1‑251 (MQ=255)                     aCAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCcac                                                                                                                                                                                                                                   >  1:108032/1‑250 (MQ=255)                      cAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCACAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCacac                                                                                                                                                                                                                                 >  2:417511/1‑251 (MQ=255)                       aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                                                                                                                                                >  1:316608/1‑251 (MQ=255)                       aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                                                                                                                                                >  1:120103/1‑251 (MQ=255)                          ccAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGt                                                                                                                                                                                                                              >  2:154579/1‑250 (MQ=255)                           cAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATACTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTa                                                                                                                                                                                                                             >  1:530941/1‑250 (MQ=255)                            aGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGACTGACGACAACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGATCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGCCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTcacaacaggtaac                                                                                                                                                                                                                           >  2:523596/1‑239 (MQ=255)                            aGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAac                                                                                                                                                                                                                           >  1:227865/1‑251 (MQ=255)                                   tgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTACCgatgac                                                                                                                                                                                                                     >  1:506848/1‑250 (MQ=255)                                     aaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGCGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGAtgacatc                                                                                                                                                                                                                  >  1:224896/1‑251 (MQ=255)                                            tgatctgaacgcacaaacccTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGCAGACACCGCCGAGACACTAATTGCC‑CTGTGTCATCAGGCCAACAAC‑CCGGCCGGCAATACCGCCGCTACCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCCAACCACTGAAAGCGCAGGGACCCCCGGAAATACGTATTGCTACGGTACCCAACTTCCCACACGGCAACGAagccatcgaaat                                                                                                                                                                                                             >  2:269142/19‑246 (MQ=255)                                             aacaagAAAGCAAGGAGCTTGGGTCCACTGAAATTGAGGGTCCTCCCCCCCGTGAATGCCAACGCCCCCGAAGAGAAATTAATTTCC‑CTTTGTCATCAGTCAAA‑ACCCCCGGCCGGCATTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATgacatcgaaatc                                                                                                                                                                                                            <  2:565572/244‑1 (MQ=255)                                                       aaagcaGCCTGAGTGCACTGAAATTGATGGACCTGACCACCATGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGAAGGGCACCCCGGAAATCCGTATTGCTACGGTAAACAACTTCACACACGGTACCGATGACATAGAAATCGCGCTGGCAGa                                                                                                                                                                                                 >  2:443218/2‑250 (MQ=255)                                                        aagcaGCCTGGGTGCACTGAAATTGATGGACCTGACCCCCCTGAAGGACGACGACACCGACGAGAAAGTATTTGCT‑CTGTGTAATCAGGCCAA‑AACCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                                                                                                                                                                               <  1:109799/251‑1 (MQ=255)                                                                 gCGTGCACTGAACTTGATGGACCTGACAACCCTGAATGACGACGACACCGACGAGAACGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCAGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTCTCCCGATTGCTCGCAAACCACTGAACGAGCAGGGCACCCCGGCAATCCGTATTGCTACGGTAACCAACTTCCCAAACGGTAACGATGACATCGAAATCGCGCTGGCAGAACCCCGTGCgg                                                                                                                                                                                      >  2:37992/1‑251 (MQ=255)                                                                      cacgaaaaTTGAGGGCCCTGACCCCCCTAAATGACGCCCACCCCTACGAAAAAGTAATGGCT‑CTGTTCCACCAGGCCAA‑AACCCCGGCCTGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTGCCCACACGGTAACTATGACAGCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                                                                                                                                 <  2:82453/246‑1 (MQ=255)                                                                      cACTGAACTTGATGGACCTGACCACCCTGACTGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                                                                                                                                 >  2:65502/1‑251 (MQ=255)                                                                           aaatggagggcccTGACCCCCCTGAATGACGACGCCACCAACGAGAAAGTAATTGCT‑TTGTGTCATCTGGCCAA‑CACCCCGGTCGGCAGTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAGTTCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTATCCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTa                                                                                                                                                                            <  1:91159/240‑1 (MQ=255)                                                                           aaaTTGATGGCCCTCACCACCCTGAATCACGACGCCACCGACGAAAAAGTATTTGCT‑CTGTGTCATCAGGCCAA‑ACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTa                                                                                                                                                                            <  2:76264/251‑1 (MQ=255)                                                                            aaTTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGAT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTAc                                                                                                                                                                           >  2:243875/1‑251 (MQ=255)                                                                                   ggccctgccccccctaaagacggacgacCCGGACGAAAAAGTAATTGCT‑CTGTGTCATCAGACAAA‑ACCCCCGGTCGACATTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccg                                                                                                                                                                    <  2:181774/229‑3 (MQ=255)                                                                                          accaccCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACAGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGt                                                                                                                                                             >  1:561909/1‑251 (MQ=255)                                                                                                    aggacgccgacACCGACGAAAAATTAATTGCT‑TTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACgtggtg                                                                                                                                                   <  1:545218/249‑1 (MQ=255)                                                                                                      gacgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGtt                                                                                                                                                 >  1:57641/1‑251 (MQ=255)                                                                                                            gacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGt                                                                                                                                            >  1:182200/1‑250 (MQ=255)                                                                                                            gacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGt                                                                                                                                            >  1:327681/1‑250 (MQ=255)                                                                                                                  gacgaGAAAGTCATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGACATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATAGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTAAGGTGCCGTTGAAGTTGACGTGGTGTTCCCGTACcgcgc                                                                                                                                     >  2:516144/1‑251 (MQ=255)                                                                                                                    cgaGAAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCACTCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCgcgcgc                                                                                                                                   >  2:574632/1‑251 (MQ=255)                                                                                                                     gagaAAGTAATTGCT‑CTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTACCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCt                                                                                                                                  >  2:463881/1‑251 (MQ=255)                                                                                                                      aagaaaggaattgctcTGTGACATCAGGCCAA‑AACCCCGGTCGGCAATACCGACGCTATCTGTATCTATCCTCGCTATATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCACACACGGTAACGATGACATCGAAATCGCGATGGCAGAAACCCGCGAGACAATCGCCTACGATGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCt                                                                                                                                  >  2:342111/16‑251 (MQ=255)                                                                                                                              aattgctcTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGAACCCCGGCAATCCGTATTGCTACGGTAACCAACTTCCCACACGGAAACGATGACATCGCAATCGCGCTGGCAGACACCCGGGCGGCAATCGCCTACGGTGCCGATGAGGTCGACGTGGTGTTCCCGTACCGCGCGCTGATGGCgg                                                                                                                          >  2:514123/8‑251 (MQ=255)                                                                                                                                ttgctcTGTGTCATCAGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGt                                                                                                                        >  1:581302/6‑251 (MQ=255)                                                                                                                                    gctgttccatcaGCCCAA‑ACCCCCGGGCGACAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCGGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACg                                                                                                                    <  2:17534/244‑1 (MQ=255)                                                                                                                                    cccgtggtctacagcacaaacccccGGTCGGAAAATCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCGACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAc                                                                                                                     <  2:518666/228‑1 (MQ=255)                                                                                                                                       gtgctatcaGCCCAA‑ACCCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTAGCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGc                                                                                                                 <  2:143696/246‑1 (MQ=255)                                                                                                                                         ttcatcaGGCCAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTCTTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCTTACCGCGCGCTGATGGCggggccatcgcac                                                                                                               >  2:214177/2‑241 (MQ=255)                                                                                                                                            atcaGGCCAA‑ACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggtt                                                                                                            <  1:283307/251‑1 (MQ=255)                                                                                                                                             tcaGGCCAA‑ACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGAGGAAGTGGACGTGGTGTTCCCGGACCGCGCGCTGATGGCGGGTAACGAGCAggtt                                                                                                            <  2:444435/250‑1 (MQ=255)                                                                                                                                                  ccAA‑AACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGCGCAGGTTGGtttt                                                                                                      >  1:440625/1‑251 (MQ=255)                                                                                                                                                     aaaacccaggTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGa                                                                                                    <  1:37992/243‑1 (MQ=255)                                                                                                                                                         ccccgttCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGGATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTg                                                                                                <  2:364108/245‑1 (MQ=255)                                                                                                                                                         ccccGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGAtgacat                                                                                                                                                                                                                   >  1:227591/3‑136 (MQ=255)                                                                                                                                                          cccGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGAGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGGCCTgg                                                                                               >  2:578897/2‑251 (MQ=255)                                                                                                                                                            cGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCCGGTTGGGTTTGCCCTgggg                                                                                             >  1:83400/1‑249 (MQ=255)                                                                                                                                                             ggTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGa                                                                                            <  2:327681/251‑1 (MQ=255)                                                                                                                                                              gTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGTAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGCAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCCGGTTGGTTTTGACCTgggga                                                                                            >  2:306995/1‑247 (MQ=255)                                                                                                                                                                 ggaaaTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGc                                                                                        <  1:229055/248‑1 (MQ=255)                                                                                                                                                                  ccAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCGCGCAAAACACTTAAAGAGCATGGCACCCCGGAAAGCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGAGGACATCGAAATCGCGCTGGCAGAATCCCGTGCGGCAAGCGCCTACGGTGCCGATGAAGTTTACGGGGGGTTCCCGTTCCGCGCGCTGATTTCGGGTAACGAGCAGGTCGGTTTTGACCTGGTGaaa                                                                                          <  2:596364/247‑1 (MQ=255)                                                                                                                                                                         cgccgcTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGCAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGCAAGCCTGTat                                                                                  >  2:109041/1‑248 (MQ=255)                                                                                                                                                                          gccgcTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAaga                                                                               <  2:514669/251‑1 (MQ=255)                                                                                                                                                                                    gTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTgcgc                                                                      <  2:329654/250‑1 (MQ=255)                                                                                                                                                                                    gTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTGCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTgc                                                                        <  1:211869/248‑1 (MQ=255)                                                                                                                                                                                    gTATCAATCCGCGCTTTATCCCGATGGCTCGCAAAACACTGAAAGAGCAGGGCACACCGGAAATCCGTATTGCTACGGTAACCAACTGCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTgcg                                                                       <  1:55042/249‑1 (MQ=255)                                                                                                                                                                                          atcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCg                                                               <  1:36154/251‑1 (MQ=255)                                                                                                                                                                                                 cTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGCACGAGCGGGTTGGTTTTGACCTGGTGTAAGCCTGGAAAGAGGCTTGCGCGGCAGCG‑AATgga                                                        >  2:25739/1‑248 (MQ=255)                                                                                                                                                                                                     aTCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AATGTActgct                                                   <  2:279751/251‑1 (MQ=255)                                                                                                                                                                                                      tCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGGAGGG‑GATGTAatgctg                                                  >  2:616023/1‑251 (MQ=255)                                                                                                                                                                                                              gCTCGCAACACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGCAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTccc                                                                                                                                              >  2:317523/1‑152 (MQ=255)                                                                                                                                                                                                              gCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTccc                                                                                                                                              <  1:317523/152‑1 (MQ=255)                                                                                                                                                                                                                  gcaaaccactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCAGGTTGGGTTTTACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AAGGTACTGGTGGAAGTGctcatc                                      >  2:111065/9‑251 (MQ=255)                                                                                                                                                                                                                            gAAATAGCAGGGCACCCCGGAAATCCGGACTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCTTGCGGCAATCGCCTACGTTGCCGATGAAGTTGACGTGGTGTTCTCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGATCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCATCG‑AATGTACTGCTGAAAGTGATCAGCGAAACCGGCg                            <  2:7900/251‑1 (MQ=255)                                                                                                                                                                                                                            gAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AATGTACTGCTGAAAGTGATCATCGAAACCgg                              <  2:58940/249‑1 (MQ=255)                                                                                                                                                                                                                              aaGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AATGTACTGCTGAAAGGGATCATCGAAACCGGCGa                           >  1:233019/1‑250 (MQ=255)                                                                                                                                                                                                                              aaGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTACCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGCAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGTGTGGTTTTGACCTGGTGAAAGACTGTAACGAGGCTTGCGCGTCAGCGGAATG‑ACTGCTGAAATTGATCATCTAAAACGGCgg                           >  2:590751/1‑249 (MQ=255)                                                                                                                                                                                                                                      ggCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGAAAGCCTGTCAAGAGGCTTGCGCGGCAGCG‑AATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGaa                     >  1:83328/1‑248 (MQ=255)                                                                                                                                                                                                                                            ccGGAAATCCGTATTGCTCAGGTAACCAACTTCCCCCACGGTAACGATGACCTCGAAATCGCGCTGGCAGCCACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGCCCTGGTGAAAGACTGTAAAGAGGCTTGGGCGGCAGCG‑AATGTAATGGCGAAAGGGATCATCGAAACCGGCGAACTGGAAGAAGAAgc            >  2:58600/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aaTCCGGATTGCTACGGTAACCAACTTCCCACACGGTATCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCTCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGa       <  1:125960/251‑1 (MQ=255)                                                                                                                                                                                                                                                    ccGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATcc    <  1:43996/251‑1 (MQ=255)                                                                                                                                                                                                                                                      gTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCCGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AATGTACTGCTGAAAGGGATCATCGAAACCGGCGAAATGAAAGACGAAGCGCTGATCCGt  >  1:406626/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        CAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAA‑AACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCG‑AATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGT  >  NC_000913/4617285‑4617778

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