Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,742	T→C	H380H (CAT→CAC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,742	0	T	C	100.0%	267.4 / NA	83	H380H (CAT→CAC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/44);  total (39/44)

TCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAA  >  NC_000913/4620495‑4620987                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      tCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGt                                                                                                                                                                                                                                                    >  1:197037/1‑251 (MQ=255)        cAACTTCGTTGACTTCGCATCTCCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTCCCGGAGCTGATGTCTCTGCGGCGCGATGACGACCTCCTGATCCTCACTGCTGACCACGGTTGGGATCCGACCTGGCCCGGTACTGACCACACGCGTGACCACATTCCGGGACTGGTTTACGGCCCGACAATCACACCGGGCTAACTGGGCCACCGTGAAACCt                                                                                                                                                                                                                                             >  2:275274/1‑251 (MQ=255)             tCGTTGACTTCGACTCTTCCGGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTcgc                                                                                                                                                                                                                                         <  2:457224/250‑1 (MQ=255)               gTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGgg                                                                                                                                                                                                                                      >  2:512661‑M2/1‑249 (MQ=255)               gTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGATGCGGGTCTGGAGCTGTACGACCGCCGTCTGACGGAGCTGTTGTCTCTGCTGTGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGATAGTAAAACCGGGCTCACTGGGTCACCGTGAACCCTTcgctgg                                                                                                                                                                                                                                      >  1:503917/1‑248 (MQ=255)                    cTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatg                                                                                                                                                                                                                                 <  1:274302‑M2/251‑8 (MQ=255)                     ttCGCCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGCAACCTTCGCGgggggtgc                                                                                                                                                                                                                                <  2:338484‑M2/251‑9 (MQ=255)                        gACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCTCGggtgatgctgc                                                                                                                                                                                                                             >  2:360831‑M2/1‑240 (MQ=255)                              tCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaactt                                                                                                                                                                                                                       <  1:140255‑M2/251‑18 (MQ=255)                              gccTGGGTCCGCCGTCGCGACGTCGCCGGTTATGCTTCTGTTCTGGAGCGGTGCGACCGCCGTCTGCCGGAGCTGATGCCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTTACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCTGTACTGGTGTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACAGTGAAACCTTCGCgggttatgctgacaactt                                                                                                                                                                                                                       <  2:503917/250‑17 (MQ=255)                                  ggggCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCAGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactg                                                                                                                                                                                                                   >  1:594595‑M2/1‑230 (MQ=255)                                   gggCCCCCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCCCTGCTGACCACGGTTGTGATTCGCCCTGGACCGGTACTGCCCACACGCGTGAACACCGTCCGGTACTGGTTTACGGTGCGACAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactga                                                                                                                                                                                                                  >  1:383542‑M2/1‑229 (MQ=255)                                    ggCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgctgctgccaacttactgat                                                                                                                                                                                                                 <  2:392225‑M2/251‑24 (MQ=255)                                     gCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGACACCTTCGCGggggatgctgccaccttactgatt                                                                                                                                                                                                                <  1:276073‑M2/251‑25 (MQ=255)                                             cgcgACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggggatgctgccaacttactgatttagtgtat                                                                                                                                                                                                        <  2:280509‑M2/251‑33 (MQ=255)                                                  cGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatgg                                                                                                                                                                                                   <  1:477435‑M2/251‑38 (MQ=255)                                                  cGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGAAGACGACATCATGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGGCACCGTGAAACCTTCGCGggtgatgctgccaacttactgaattagtgtatgatgg                                                                                                                                                                                                   <  1:282347‑M2/251‑38 (MQ=255)                                                      gCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtt                                                                                                                                                                                               >  2:49901‑M2/1‑210 (MQ=255)                                                            tATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttg                                                                                                                                                                                           <  1:546654‑M2/249‑46 (MQ=255)                                                               gctgcgggtCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggt                                                                                                                                                                                       >  2:197037‑M2/4‑201 (MQ=255)                                                                    gggTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctcca                                                                                                                                                                                 >  2:390006‑M2/1‑196 (MQ=255)                                                                     ggTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccag                                                                                                                                                                                >  1:235811‑M2/1‑195 (MQ=255)                                                                           gagctgtTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgcgccagtggctt                                                                                                                                                                          >  1:493425‑M2/4‑189 (MQ=255)                                                                            agctgTTCGACCGCCGTCTGCCCGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGATCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttc                                                                                                                                                                         <  1:429554‑M2/249‑64 (MQ=255)                                                                             gctGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttct                                                                                                                                                                        >  1:76531‑M2/2‑187 (MQ=255)                                                                                    gACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGCAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggggatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttct                                                                                                                                                                  <  1:538169‑M2/250‑71 (MQ=255)                                                                                              tGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtcc                                                                                                                                                       <  1:297408‑M2/251‑82 (MQ=255)                                                                                               gCCGGAGCTGATTTCTCTGCTGCGCGATGACGACATCCTGATCTTCACTGCTGACCTCGGTTGAGATCCGACCTGGACCGGTACTGACCACAAGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgttttggaggtgctccagtggcttctgtttctatcagctgt                                                                                                                                                         >  2:129475‑M2/1‑169 (MQ=255)                                                                                                 cGGAGCTGATGTCTCGGCTGCGCGATGACGGCATCCTGAGCCTCACTGGTGACCACGGTTGTGAGCCGACCTTGACCGGTACTGACCACACGCGTGAACCCATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtgtttgaggtgctccagttgcttctttttctatcagctgtccctc                                                                                                                                                    >  2:389566‑M2/1‑167 (MQ=255)                                                                                                        gATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcc                                                                                                                                                   <  1:374873‑M2/245‑86 (MQ=255)                                                                                                        gATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcc                                                                                                                                                   >  2:374873‑M2/1‑160 (MQ=255)                                                                                                          ggtCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgagggtgttttggaggtgcgccagtggcttctgtttctatcagctctcccgcctgttcagc                                                                                                                                           >  2:409528‑M2/2‑158 (MQ=255)                                                                                                               ctgctgCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgcgccagtggcttctgtttctatcagctgtccctcctgttcagctac                                                                                                                                        >  2:546654‑M2/1‑153 (MQ=255)                                                                                                                 gctgcGCGATGACGACATCCTGCTCCTCACTGCTGACCGCGGTTGTTATCCGACCTGGTCCGGTACTGTCCACACGCGTGAACTCATTCCGGTACTGTTTTACGGCCCGAAAGTCAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttacttatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctccggttcagctact                                                                                                                                       >  2:184498‑M2/1‑151 (MQ=255)                                                                                                                   tgcGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacgg                                                                                                                                  <  1:260126‑M2/251‑103 (MQ=255)                                                                                                                        gATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTAATGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtgg                                                                                                                             <  2:419059‑M2/251‑108 (MQ=255)                                                                                                                               acATCATGATCCTCCCTGCTGAACACGGTTGTGCTCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgctttagtgtatggtggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtg                                                                                                                       <  1:168843‑M2/250‑114 (MQ=255)                                                                                                                                cATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaac                                                                                                                     <  1:373666‑M2/251‑116 (MQ=255)                                                                                                                                cATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACCGGTTTACGGCCCGACAGTACAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgacaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtacctcctgttcagctactgacggggtggtgcgtaac                                                                                                                     <  1:398772‑M2/251‑116 (MQ=255)                                                                                                                                     tGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggca                                                                                                                 <  1:152807‑M2/250‑120 (MQ=255)                                                                                                                                        tCCTCAGTGCTGACAACGGTTGTGATCCGCCCTGGACCGGTACTGAAAACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGAAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaaattccagatttagtgtatgatggtgtttttgaggtgctccagtggcttttgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggaaca                                                                                                               <  2:555308‑M2/249‑122 (MQ=255)                                                                                                                                        tCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggcggtgcgtaacggcaaa                                                                                                               <  1:375793‑M2/249‑122 (MQ=255)                                                                                                                                           tCACTGCTGAACACGGTTGTGATCCGCCCTGGACCGGTACTGCCCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcac                                                                                                          <  2:592880‑M2/251‑127 (MQ=255)                                                                                                                                             actgctgACCACGGTTGTGATCTGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTCCTGGGTTAAGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGAGACACCTTCGCGggtgatgctgcaaacttactgctttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccg                                                                                                        <  1:405525‑M2/248‑129 (MQ=255)                                                                                                                                             actgctgACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccg                                                                                                        >  1:549050‑M2/4‑123 (MQ=255)                                                                                                                                                gCTGACCACTGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgcc                                                                                                      >  1:267927‑M2/1‑120 (MQ=255)                                                                                                                                                 cTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccgg                                                                                                    <  1:171702‑M2/251‑133 (MQ=255)                                                                                                                                                    aCCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggac                                                                                                  <  1:334784‑M2/250‑135 (MQ=255)                                                                                                                                                         ggTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagc                                                                                            >  1:233983‑M2/1‑111 (MQ=255)                                                                                                                                                          gTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcg                                                                                           >  1:587815‑M2/1‑110 (MQ=255)                                                                                                                                                                  ccGACCTGGACCGGTACTGACCACACGCGTGAACCCATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctct                                                                                   <  2:557364‑M2/251‑150 (MQ=255)                                                                                                                                                                       cTGGACCGGTACTGCCCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgctttagtgtatgatggtgtttttgaggagctccagtggcttctgtttctatcagctgtccctcctgttcagctactgtcggggtggggcgtaacggcacaagcaccgccggacatcagcgctatctctgccct                                                                              <  2:382040‑M2/251‑155 (MQ=255)                                                                                                                                                                         ggACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgcgccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctca                                                                            >  2:361189‑M2/1‑95 (MQ=255)                                                                                                                                                                            ccGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactg                                                                         <  1:424735‑M2/251‑160 (MQ=255)                                                                                                                                                                               gTACTGACCACACCCGTGAACCCATTCCGGTACTGGTTTACGGCACGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggggatgctgcaacttactgatttagtgtatgatggtgtttttgaggtgccccagtggcttcggtttctatcagctgtccctcctgttcagctactgtcggggtggtgcgtaacggtaaaagcaccgcttgacatcagcgctatctctgctctcactgccgt                                                                      >  1:360543‑M2/1‑89 (MQ=255)                                                                                                                                                                                      ccACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaaca                                                               >  1:259215‑M2/1‑82 (MQ=255)                                                                                                                                                                                         cacGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatgg                                                            <  1:251304‑M2/251‑173 (MQ=255)                                                                                                                                                                                         cacGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctggccctcctgttcagctactgacggggtggtgcgtaacggccaaagcaccgccggacataagcgatatcactgctctcactaccgtaaaacatgg                                                            <  2:193515‑M2/251‑173 (MQ=255)                                                                                                                                                                                             cgTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaac                                                        >  1:144981‑M2/1‑75 (MQ=255)                                                                                                                                                                                               tGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatc                                                                                                                                                               <  1:442374‑M2/146‑74 (MQ=255)                                                                                                                                                                                               tGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgagggtgtttttgaggtgctccagtggcttctgtttctatc                                                                                                                                                               >  2:442374‑M2/1‑73 (MQ=255)                                                                                                                                                                                                gAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgc                                                     >  1:4471‑M2/1‑72 (MQ=255)                                                                                                                                                                                                gAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactg                                                      >  1:563932‑M2/1‑72 (MQ=255)                                                                                                                                                                                                 aaCACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTTGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttttatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcacaagcaccgccggacatcagcgctatctctgctctaactgccgtaaaacatggcaactgcc                                                    <  2:215503‑M2/251‑181 (MQ=255)                                                                                                                                                                                                   cacaTTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggagatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaaaggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagt                                                  <  2:316947‑M2/251‑183 (MQ=255)                                                                                                                                                                                                    acaTTCCGGTATTGGTTTACGGCGCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAATCCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgcgccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggggcgtaccggcaaaagcaccgccggacatcagcgctagctctgctctcactgccgtaaaacatggcaacggcagtt                                                 >  2:276073‑M2/1‑68 (MQ=255)                                                                                                                                                                                                        tCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggggtttttgaggtgctccagtggcttctggttctatcagctgtccctcctgttcagctactgacggggtggggcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctatctgccgaaaaaaatggcgaatgcagttcact                                             <  1:332853‑M2/251‑188 (MQ=255)                                                                                                                                                                                                         ccGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcactt                                            >  2:477435‑M2/1‑63 (MQ=255)                                                                                                                                                                                                             tACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaccatggcaactgcagttcacttacac                                        <  1:446380‑M2/251‑193 (MQ=255)                                                                                                                                                                                                                tGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgc                                     <  1:154459‑M2/251‑196 (MQ=255)                                                                                                                                                                                                                   tttacggcccgaAAGTAACACCGGGCTCACTTGGTCACCGTTAAACCTTCGCGggagatgctgccaacttactgatttagggtatgatggggtttttgtggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacgtcacaagcaccgccgtacctcagcgctatctcttctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttc                                  <  2:277444‑M2/246‑199 (MQ=255)                                                                                                                                                                                                                   tttacggcccgaAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggccatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttc                                  <  2:277438‑M2/246‑199 (MQ=255)                                                                                                                                                                                                                      acggccCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgct                                                                                         >  2:490636‑M2/3‑50 (MQ=255)                                                                                                                                                                                                                      acggccCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgct                                                                                         <  1:490636‑M2/191‑144 (MQ=255)                                                                                                                                                                                                                       cggcCCGAAAGTAAAACCGGGTCCACTGGGTCACCGTGAATCCTTCGCGggtgatgctgccaacttacggatttagtgtatgatggtgttttggaggtgctccagtggcttctgtttctatcagctgtccctccggttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgctgctcaa                              >  2:202734‑M2/2‑49 (MQ=255)                                                                                                                                                                                                                        ggCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggcggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaac                             >  2:375793‑M2/1‑48 (MQ=255)                                                                                                                                                                                                                          cccGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgagggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaaccc                           >  1:187183‑M2/1‑46 (MQ=255)                                                                                                                                                                                                                                aGTAAAACCGGGCTCAATGGGTCACCGGGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgccccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggttcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcatcccggtacg                     >  2:569303‑M2/1‑40 (MQ=255)                                                                                                                                                                                                                                         gggCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacctcagcgctatctctgctctcactgccgtaaaccctggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaa             <  2:577918‑M2/250‑220 (MQ=255)                                                                                                                                                                                                                                             tCACTGGGTCACCGTGAAACCTTCGCGgttgatgctgccaacttaatgatttagtgtatgatggtgtttttgaggcgctccagtggcttctgtatctatcagctgtccctcctgttcagctaatgacggggtgatccgtaacggcaacagcacagccggacatcagcgctatctctgctctcactgccgttaaacatcgcaactgcaattcacttacaccgcttctcaacccggtacgacccagaaaatc         <  1:260706‑M2/250‑224 (MQ=255)                                                                                                                                                                                                                                                  gggTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgccctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattg     >  2:540357‑M2/1‑22 (MQ=255)                                                                                                                                                                                                                                                   ggTCACCTTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattgata  >  1:338484‑M2/1‑21 (MQ=255)                                                                                                                                                                                                                                                   ggTCACCGTGAAACCTTCGCGcgtgatgctgccaacttactgattcagtgtatgatggtgcttttgaggtgctccagaggcttctgtttctctcagctgtccctcctgttcagatactgacggggtggtgcgtaagggcaacagcaccgccggacttcagcgctatctctgctctcactgccgtaaaaaatggcaactgcagttcacttacaccgcttctcaacacggtccgcaccagaaaatcattatt   <  2:606877‑M2/250‑230 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      TCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAA  >  NC_000913/4620495‑4620987

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: ‑