Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,709	T→C	Y369Y (TAT→TAC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,709	0	T	C	98.8%	275.3 / ‑4.1	85	Y369Y (TAT→TAC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (39/45);  minor base A (0/1);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

TGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGC  >  NC_000913/4620468‑4620952                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    tGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGaa                                                                                                                                                                                                                                             >  1:512661/1‑250 (MQ=255) ggAAAAAAGCGGCTGATACCACCATCGTCTCCACCAACGTCGTTGACGTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCTGTTATGCTTCGGGTCTGGAGCTGTGCGTCCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCCCGCTGACGACATCCGGATCCTCACTGCTGACCACGGTTTTGAGCCGACCTGGACCGGTACTGTCCACGCGCGTGACCACATTGCGGTACTGGTTTACGGCCTGaaa                                                                                                                                                                                                                                            <  2:56211/250‑1 (MQ=255)   aaaGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCATGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGCCCGGTACTGACCACACGCGTGAACACCTTCCGGAACTGGTTTACGGCCCGAAAGt                                                                                                                                                                                                                                          >  1:389566/1‑251 (MQ=255)           gggTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCgg                                                                                                                                                                                                                                  <  1:559371/251‑1 (MQ=255)             gTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGc                                                                                                                                                                                                                                <  2:345695/251‑1 (MQ=255)               gATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTc                                                                                                                                                                                                                              <  1:156994/251‑1 (MQ=255)               gATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTACGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTc                                                                                                                                                                                                                              >  2:309750/1‑251 (MQ=255)                aTAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTc                                                                                                                                                                                                                              >  1:393586/1‑250 (MQ=255)                   acacCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATGCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTg                                                                                                                                                                                                                          <  2:96920/251‑1 (MQ=255)                     acCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTggg                                                                                                                                                                                                                        >  1:171225/1‑251 (MQ=255)                     acCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTggg                                                                                                                                                                                                                        <  1:287228/251‑1 (MQ=255)                        aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGCGCCACCGTCGCGACGTCGACGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGCTGTCTCTGCTGCGCGCTGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGCCCTGGACCGGTCCTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCACGAAAGTAAAACCGGGCTCACTGGGTc                                                                                                                                                                                                                      >  2:547830/1‑250 (MQ=255)                          cggcTTCACCAACTTCGTTGACTTCGACTCTCCCTGGGGCAACCGTCGCGACGTCTCCGGGTATGCTGCGGGTCTGGAGCTGTTCGCCCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAACCCGGGCTCACTGGGTCacc                                                                                                                                                                                                                   <  2:140476/248‑3 (MQ=255)                            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)                                             gACTTCGACTCTTCCTGGGGCCCCCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCa                                                                                                                                                                                                                     <  2:179997/230‑1 (MQ=255)                                             gACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCa                                                                                                                                                                                                                     >  1:179997/1‑230 (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                                                                                                                >  2:374873‑M2/1‑160 (MQ=255)                                                                                                                                   gATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcc                                                                                                                <  1:374873‑M2/245‑86 (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)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    TGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGC  >  NC_000913/4620468‑4620952

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