Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,026	A→G	E266E (GAA→GAG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,026	0	A	G	96.8%	184.9 / ‑2.5	62	E266E (GAA→GAG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (29/31);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

CCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATA  >  NC_000913/4618779‑4619267                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  ccGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGCATCAGGTACTGGCCTCCCGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAgtat                                                                                                                                                                                                                                                >  1:363682/1‑251 (MQ=255) ccGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAgtat                                                                                                                                                                                                                                                <  1:262023/251‑1 (MQ=255)  cgccgcTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATc                                                                                                                                                                                                                                               <  2:389216/251‑1 (MQ=255)        tATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCCGGTAACGCGGTTGAAGTTCGTGAAGCGGGGC‑AGTTCCTGACGGGTGGGTATCGTAAcc                                                                                                                                                                                                                                         >  1:112292/1‑251 (MQ=255)         atgcTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTaa                                                                                                                                                                                                                                           <  1:176533/245‑1 (MQ=255)         aTTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCTCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAAGCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAAc                                                                                                                                                                                                                                          <  2:396823/249‑1 (MQ=255)            cTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCCCCGACATGAATCAGGTACTGGCCTCCAGTGCCGGGAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGGGAGTATCGTAACCcgcg                                                                                                                                                                                                                                     >  2:45486/1‑251 (MQ=255)               gCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCt                                                                                                                                                                                                                                  <  1:393980/251‑1 (MQ=255)               gCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCACTT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCt                                                                                                                                                                                                                                  <  1:284241/251‑1 (MQ=255)                  aagaaACTGGCGGAAGGTCTGGACGCGCTGGTGTTGGACGTGAAAGTGGGTAGCCGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGTATCAGGGACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGtt                                                                                                                                                                                                                               >  1:398955/1‑251 (MQ=255)                                   tCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGACGAAT‑TTATGCCGACCTACGAACTCTCTGAAACCCTTGCCGACGCGATTGTTGGCGTGGCTAACGGCGCTGGCGGGCGCACCACCGCGCTGCTCACCGAAATGAATCAAGTACTGGCCTTCAGGGCAGGTAACGCGGTTGGAGTTCGTGAAGCGGTGA‑CGTGCCTGACGGGTGAGTATCGTAACCCGCGTCGGTTTGCAGTCACGCTggagg                                                                                                                                                                                                              >  1:133399/1‑248 (MQ=255)                                    cTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTAGCCTCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑CATTCCTGACGGATGAGAATCGTAAACCGCGACTGATTGATGttccgtgggagct                                                                                                                                                                                                             >  1:240151/1‑239 (MQ=255)                                     tGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGGTGGAGTTCGGGAAGCGGTGC‑AGTTCCTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGGTGTAACGATGGGGCtt                                                                                                                                                                                                            >  2:234984/1‑250 (MQ=255)                                       gACGCGCTGGTGATGGACGTGAAAGGGGGTGGGGGCG‑GTGTTTTGCCGACCTACTAACGCTGTGAAGCCCTTGCCTAAGCTAGTGTTGGCGTGTCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTTCCCGACTTGAATCAGGTACTGGCCGCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCGGACGGGTGAGTATCGTAACCCGCGTCTGTTGGATGTCACGATTGCGCtgtt                                                                                                                                                                                                          <  2:274852/251‑2 (MQ=255)                                           cgcTGGTGATGGACGGGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTg                                                                                                                                                                                                      <  1:103966/251‑1 (MQ=255)                                              gggTGATGGACTTGAAAGTGGGTAGCGGCGCTT‑GTATGCCGACCTCCGAACTCGCTGAAGCCCCTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTTCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCTGTGGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGgag                                                                                                                                                                                                   <  2:406221/250‑1 (MQ=255)                                                     ggACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCtt                                                                                                                                                                                             <  1:162307/250‑2 (MQ=255)                                                     ggACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCtt                                                                                                                                                                                             >  1:24726/1‑249 (MQ=255)                                                         gTGAAAGGGGGCAGCGGCGCGT‑TTATGCCGACCTCCGAACTCCCTGAAGCCCGTGCCGAATCGAGTGATGGCGTGGCTATCGGCGCTGGCGTGCGCACCGTCGCGCGGCTCACCGACATGAATCAGTTACTGGCCTCCAGTGCAGGTAACGCGGTGGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTATCCCTCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTAtctc                                                                                                                                                                                        <  2:379735/251‑1 (MQ=255)                                                                  ggaaGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGAACTGGCCTCCAGTGCAGGACACGCGGTTGAAGTTCGTGAAGCGGTGCCAGT‑CCTGGCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTAACGCTGGCGCCGTGCGTGGAGATGCTTATCTCCGGCaaaat                                                                                                                                                                               >  2:125511/4‑249 (MQ=255)                                                                                         ccTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACGGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑CGTTCCTGACGGGAGAGTATCGTAACCCGGGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAcgcg                                                                                                                                                         >  2:172383/1‑251 (MQ=255)                                                                                          cTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGAAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgc                                                                                                                                                        >  1:451083/1‑251 (MQ=255)                                                                                           tACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcg                                                                                                                                                       >  1:10154/1‑251 (MQ=255)                                                                                            aCGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcgc                                                                                                                                                      >  1:237042/1‑251 (MQ=255)                                                                                            aCGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcg                                                                                                                                                       <  1:277250/250‑1 (MQ=255)                                                                                               aaCTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGcgcgcgcg                                                                                                                                                   <  1:275876/251‑1 (MQ=255)                                                                                                                gCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCt                                                                                                                                  >  1:198413/1‑251 (MQ=255)                                                                                                                    aaGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGAc                                                                                                                              >  1:93133/1‑251 (MQ=255)                                                                                                                        gATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACg                                                                                                                          <  1:185994/251‑1 (MQ=255)                                                                                                                         gttgttgGCGTGGCTAACGGCGCTGGCGTGCGCCCCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACgg                                                                                                                         <  2:59938/250‑1 (MQ=255)                                                                                                                         gttgttgGCGTGGCTAACGGCGCTGGCGGGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACgg                                                                                                                         <  2:150030/250‑1 (MQ=255)                                                                                                                              tgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTaaa                                                                                                                     <  1:481296/250‑1 (MQ=255)                                                                                                                                            gggttGGCGTGCCCACCGCCGCGCGGCTCATCGCCATGAACTAGGTACTGGCCTCCAGTGCAGGGAACGCGGTTGAAGTTCGTGAAGGGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCttt                                                                                                      <  2:254799/247‑1 (MQ=255)                                                                                                                                               cTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGt                                                                                                   >  1:129078/1‑251 (MQ=255)                                                                                                                                                 ggCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGt                                                                                                   <  1:45486/249‑1 (MQ=255)                                                                                                                                                          accaccGCGCTGCTCACCGACATGAATCGGGTACTAGCCTCCAATGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTCACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGCAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGCAGTCTTTGGTCGTATGGTAGc                                                                                        >  1:428508/1‑251 (MQ=255)                                                                                                                                                              ccGCGCTGCTCACCGACATGAAACAGGTACGGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGGGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTATCGGca                                                                                    <  2:200148/251‑1 (MQ=255)                                                                                                                                                                   ctgctCACCGATATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCATGAAGCGGTGC‑AGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGGCTTTGGGCGTATGGTAGCGGCACaaa                                                                                >  1:98242/1‑250 (MQ=255)                                                                                                                                                                         cccgacaggagttaGGTACTGGCCACCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGTccc                                                                         <  2:319490/238‑1 (MQ=255)                                                                                                                                                                             aCATGAATCGGGTACTGGCCTCCATGGCAGGTAACGCGGTGGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTATCCCGCTGCTGTTTGATGTCACGATGGCGCCGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCccaacc                                                                     <  2:469309/251‑1 (MQ=255)                                                                                                                                                                             aCATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCccaacc                                                                     <  1:355513/251‑1 (MQ=255)                                                                                                                                                                               aTGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGACGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGGCTTTGGTCGTATGGTAGCGGCACACAAAGGCccaaccg                                                                    >  2:407745/1‑250 (MQ=255)                                                                                                                                                                                   aTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGGTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTc                                                               <  2:10154/251‑1 (MQ=255)                                                                                                                                                                                     cAGGTTCTGGCCGCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCACCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCGGTGGTCGTATGGTAGCGGCACAAAAAGGACCCACCGACTTCg                                                              <  2:59595/250‑1 (MQ=255)                                                                                                                                                                                                cTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGGCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGAGGAGATGCTTCTCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAACGGCCCAACCGCCTTCGTTGAGAACTAcg                                                  >  2:81056/1‑251 (MQ=255)                                                                                                                                                                                                  ccAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTAcgcg                                                >  1:1083/1‑251 (MQ=255)                                                                                                                                                                                                      tGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGt                                            >  1:6987/1‑251 (MQ=255)                                                                                                                                                                                                      ggtagggaacgcggtgtaattgcgggAAGCGGTGC‑TTTTCCCGACGGGGGAGTATCGTAACCCGCGTCTTTTTGATGTCACGATGGCGCTGGGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGt                                            <  2:399139/226‑1 (MQ=255)                                                                                                                                                                                                        cAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTAt                                          >  1:478238/1‑251 (MQ=255)                                                                                                                                                                                                         aGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATc                                         >  1:435103/1‑251 (MQ=255)                                                                                                                                                                                                          ggTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGAAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAACAGGCCCCACCGACTTCGTTGCGAACTACGCGAGGTATCt                                        >  2:267944/1‑251 (MQ=255)                                                                                                                                                                                                           gTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTg                                       <  2:24726/251‑1 (MQ=255)                                                                                                                                                                                                             aaCGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACAGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                     >  1:90088/1‑251 (MQ=255)                                                                                                                                                                                                                 gggTTGAAGTTCTGGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACa                                 <  2:359745/250‑1 (MQ=255)                                                                                                                                                                                                                         gTTCGTGAAGCGGTGC‑AGTTCCTGACGTGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCt                         <  2:331050/251‑1 (MQ=255)                                                                                                                                                                                                                         gTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCt                         <  2:155417/251‑1 (MQ=255)                                                                                                                                                                                                                         gTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCt                         >  1:331050/1‑251 (MQ=255)                                                                                                                                                                                                                            cGTGAAGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGAc                      >  1:108423/1‑251 (MQ=255)                                                                                                                                                                                                                             atGAAGCGGTGC‑AGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACg                     <  2:98242/250‑1 (MQ=255)                                                                                                                                                                                                                                aaGCGGTGC‑AGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCCCAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGCAGTATCTGCCGACAGCGATGCTGACGaa                   >  2:248331/1‑250 (MQ=255)                                                                                                                                                                                                                                            ttCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGCTGTCACGATGGCGCTGGGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGCCGAAAGAAGTCTATGc       >  2:483270/1‑251 (MQ=255)                                                                                                                                                                                                                                                tGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGAt   >  1:280967/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGa    <  1:234984/249‑1 (MQ=255)                                                                                                                                                                                                                                                 gACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATa  <  2:196691/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  CCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGT‑TTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGC‑AGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATA  >  NC_000913/4618779‑4619267

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