Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,621,236	C→T	R118R (CGC→CGT)	deoD →	purine nucleoside phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,236	0	C	T	97.6%	132.6 / NA	43	R118R (CGC→CGT)	deoD	purine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/18);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAAT  >  NC_000913/4620993‑4621459                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                ccGTGAAGGGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGGGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTAGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGGTAACCGTATCCGtt                                                                                                                                                                                                                          >  2:489808/1‑251 (MQ=255)              aacGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGa                                                                                                                                                                                                             <  2:191584/251‑1 (MQ=255)                 gttagcgGTATGCTGGGTTTCACCGGTATTTACAAAGTCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGTTCCCCGAGTGCGTCGTGAAGAAAATTATCCGCGTGGGTTCCGGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCTTTATCGGTATGGGTTCCGGCCCCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATgccgt                                                                                                                                                                                                          <  2:264961/247‑5 (MQ=255)                    cgcgGTATGTTAGTCTTCACCTGTACTGTCAAAGTCCTCAAAATTTCCCTAATGGTTCTCGTTATGCGTAACCCGTCCTTCTCCATCTACAACATAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCTGCGTGGTTTCCTGTGGCGCAGCTCTGCCGGACGTAAAACTGCTCGACGTCGTTATCGGTATGGTTGCCTTCCTCGATTCCACAATTAACCTTATCCGTTTTAAAGCCCATGACTTTGc                                                                                                                                                                                                       <  2:46373/251‑1 (MQ=255)                              tGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAACCTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGGATACGTTTTAAAGACCATGACTTTGCAGCTATCGCt                                                                                                                                                                                             >  1:97286/1‑251 (MQ=255)                                ggCTTCACCGGTACTTACAAAGTCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGTCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCTCCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGa                                                                                                                                                                                           <  2:408832/251‑1 (MQ=255)                                                   aaGGCCGCAAAATTTCCGTAATTGGTCACGGTCTGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCCGTTCTGCCTCACGTAACACTGCGCGACGTCGTTCTCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCTCTATCGCTGACTTTGACATGGTCCGTAAc                                                                                                                                                                        >  2:393051/1‑251 (MQ=255)                                                     ggCCGCAAACTTTCCGTAATGGGTCACGGTATGCGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGCTCACCGATTTCGGCGTGCAGAAAATTATCCGCGTGGGTTCCTGTGGCGCTGTTCTGCCGCACGTACAACTGCGCGACGTCGTTATCGGTCTGGGTGCCTGCACCGATTCCACAGTTAACCGTATCCGTTTTAAAGGCAATGACTTTGCCGCTATCACTGGCTTTGACCTGGTGGGTAACgg                                                                                                                                                                      >  2:216769/1‑250 (MQ=255)                                                        cGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCGTCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGACGCTATCGCTGACTTTGACATGGTGCGTAACGcggt                                                                                                                                                                   >  1:148278/1‑248 (MQ=255)                                                            aaaatcccatAATTGGTCACTGTATTGGTATCACGTCCAGCTCCATCTACACCAAAGAACTGATCACCGATTGCGGCGTGAAGAAAATTATCCGCGTGGGTGCCTGTGGCGTATTTCTGCCGCACGTAAAACTTCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCGATCGCTGACTTTGACAGGGTGCGTACCGCGGTAGac                                                                                                                                                               <  2:320015/242‑2 (MQ=255)                                                                   cGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAACTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGATTTCGCTGACTTTGACATGGTGCGTAACGCGGTAGacgaggct                                                                                                                                                         >  2:171812/1‑243 (MQ=255)                                                                              aCGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTTAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATTGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGt                                                                                                                                             <  1:315684/251‑1 (MQ=255)                                                                                         aTCCCGTCCTGCTCCATCTACACCAAAGAACTGATCCCCGAGTGCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTTCCGCACGTAAAACTGCGCGACGTCGTTAGCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTTTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTc                                                                                                                                   <  2:294848/250‑1 (MQ=255)                                                                                                ccTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGt                                                                                                                           >  1:326775/1‑251 (MQ=255)                                                                                                  tGCTCCATCTACAACAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTa                                                                                                                          >  1:367154/1‑250 (MQ=255)                                                                                                   gCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTGGAGGCTCGCGTGGGTAAc                                                                                                                        <  1:142056/251‑1 (MQ=255)                                                                                                     tCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGAGTGGGTAACCt                                                                                                                      >  1:340995/1‑251 (MQ=255)                                                                                                      ccATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGGTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTg                                                                                                                     <  2:469876/251‑1 (MQ=255)                                                                                                           tACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTtctc                                                                                                                >  1:362736/1‑251 (MQ=255)                                                                                                                 aaaGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTg                                                                                                           >  1:106409/1‑250 (MQ=255)                                                                                                                  aaGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGAc                                                                                                         >  1:397611/1‑251 (MQ=255)                                                                                                                   aGAACTGATCCCCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGAcc                                                                                                        <  1:380078/251‑1 (MQ=255)                                                                                                                              ccGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGGTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTActct                                                                                             <  1:487621/251‑1 (MQ=255)                                                                                                                               cGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGCCCTGTTCTACtctc                                                                                            >  1:226388/1‑251 (MQ=255)                                                                                                                                               aaaaTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGtt                                                                            <  1:210113/251‑1 (MQ=255)                                                                                                                                                    tATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGAc                                                                        <  2:106409/250‑1 (MQ=255)                                                                                                                                                        cgcgTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAACGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAATCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGAt                                                                   >  1:182984/1‑251 (MQ=255)                                                                                                                                                        cgcgTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGCCCTGTTCTACTCTCCGGACGGCGCACTGTTCgccgcga                                                                    >  2:100621/1‑244 (MQ=255)                                                                                                                                                                  ccTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATAc                                                         >  1:340948/1‑251 (MQ=255)                                                                                                                                                                  ccTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGGTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTGTTAAAGACCATTACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGTTGCTCGCGTGGGTAACCTGTTCTCCGCGGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATAc                                                         <  1:501456/251‑1 (MQ=255)                                                                                                                                                                     gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATgg                                                                 >  1:308409/1‑240 (MQ=255)                                                                                                                                                                     gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATgg                                                                 <  2:308409/240‑1 (MQ=255)                                                                                                                                                                         cgcAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCAtcc                                                  <  1:277667/251‑3 (MQ=255)                                                                                                                                                                             gTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTCGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGCTGGAAGAATACTGCATCCTCgg                                              >  2:116660/1‑251 (MQ=255)                                                                                                                                                                                            aaaCTGCTCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCAGCCGTTGTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCGAAAGCGCTGGGTATTGATGCGCGCGGGGGTAACCTGTGCTCCGCTGACCTGTTCGACTCGCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGGGGAAATGGAAGc                               <  2:49904/251‑1 (MQ=255)                                                                                                                                                                                               cTGCGCGACGTCGTTATCGGTATGGGAGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCACTGACCTGTTCTACTCTCCGGACGGCGAACTGTTCGCCGTGGTGGGAAAATAGCGCATCCTCGGCGGGGAAATGGATGCGGc                            >  2:198614/1‑251 (MQ=255)                                                                                                                                                                                                        gtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGCTGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCt                    >  1:83385/1‑250 (MQ=255)                                                                                                                                                                                                             tATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTTTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGCACTGGAAGCGGCTggggtttatggcg              >  2:356145/1‑240 (MQ=255)                                                                                                                                                                                                                cGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGCCCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGGATCTAAGGCGTcc           >  2:71460/1‑250 (MQ=255)                                                                                                                                                                                                                     tGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGc       >  1:258474/1‑250 (MQ=255)                                                                                                                                                                                                                         ggCCTTCACCGATTCCAAAGTTAACCGTATCCGGTTTAAAGACCATGACGTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGGTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCGCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCGGCAGagt  <  2:207056/250‑3 (MQ=255)                                                                                                                                                                                                                                    aTTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGaaaa                                                             <  1:497244/181‑1 (MQ=255)                                                                                                                                                                                                                                    aTTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGaaaa                                                             >  2:497244/1‑181 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                CCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAAT  >  NC_000913/4620993‑4621459

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