Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,702	T→C	S158S (AGT→AGC)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,702	0	T	C	96.2%	133.0 / ‑4.4	53	S158S (AGT→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (32/18);  total (35/18)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAC  >  NC_000913/4618458‑4618947                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           ggtaTGGCCCGAGTGGTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACAGTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTg                                                                                                                                                                                                                                                      <  2:346321/248‑1 (MQ=255)       gCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACCTTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGGCAGACCAGCCCACTGGCTccc                                                                                                                                                                                                                                                >  1:434747/1‑250 (MQ=255)               gttgATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACGGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGTCATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATaa                                                                                                                                                                                                                                        <  1:465840/251‑1 (MQ=255)                  gATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACg                                                                                                                                                                                                                                     <  1:336345/251‑1 (MQ=255)                            ccACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGGTTCCGC‑GGAAATATTAAAAACGTCGGCGTGGCGGTTTTCGGTCAGAACAACTCCCTGGCTCCGGCGGATAAAAGttttttagcg                                                                                                                                                                                                                           >  1:256210/1‑244 (MQ=255)                             cactggTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGGTGAAAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGGGGCGATTATCGGTCAAAACAGCCCACTGGCTCCGGCGGGGAAACGTTTCTAcg‑cg                                                                                                                                                                                                                          >  1:295400/5‑250 (MQ=255)                                 ggTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGAcc                                                                                                                                                                                                                       <  1:439993/249‑1 (MQ=255)                                            cGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAATGGAATCCATCCCTGGCTTCGACATTTTTCC‑GGGTGACAACCGTTTCCGC‑GAAATTATTTAAGACGTCGGCGTGGCGATTATCGGTCAGAACAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCg                                                                                                                                                                                                          >  1:361638/1‑251 (MQ=255)                                                gtgACTGCTCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAAc                                                                                                                                                                                                      <  2:434529/251‑1 (MQ=255)                                                      tCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTgg                                                                                                                                                                                                 <  1:680775/250‑1 (MQ=255)                                                       cgcggtggTTGGGGCCGATGGTCGCCGCCTTCGGGGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTGCCGC‑GAAATTATTAAAGACGTCGGCGTTGCGATTATCGGGCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGTCCCGTGATATTACCGCAACCGTGGAc                                                                                                                                                                                               <  2:434747/244‑1 (MQ=255)                                                           gATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑Tcc                                                                                                                                                                                           <  1:443023/250‑1 (MQ=255)                                                            aTGTGGGGGCCGCTGGTCGCAGCCTGCGGTGGCTATATTCCGATGATCTCTTGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGAAAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCaa                                                                                                                                                                                         <  2:255932/251‑2 (MQ=255)                                                             tGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGCCATTTTCCC‑GGCTGACAACCGTTTCCGC‑GGAAATATTAAAGACGTCGGCGTGGCGGTTATCGGTCAGACCAGCTCACCGGCTCCCGCGGATAAACGTTTCTACG‑CGGCCCCCGCTATTACCGCAACAGTGGAC‑TCCaac                                                                                                                                                                                        >  1:278250/1‑249 (MQ=255)                                                                     cctctggtcgtgtccTGCTGCGGCTATAGTCCGATGATCTCTGGTCCCGGCCGCGGTCAGACTGGCGGTACTCTCGACAAACTGGAAGCCATCCCGGGCTCCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGATCCGTAATATTACCGCAACCGTGGAC‑TCCATCCCGCTGAt                                                                                                                                                                                <  2:202291/238‑1 (MQ=255)                                                                             cGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCct                                                                                                                                                                        <  1:411470/251‑1 (MQ=255)                                                                                       ggcggcTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGAACAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGGTATTAACGCAACCGGGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGc                                                                                                                                                              >  1:35302/1‑251 (MQ=255)                                                                                                  tCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCC‑ACTGATCACCGCCTCTATTCTGGCGAAGAAACTTGc                                                                                                                                                  >  1:334433/1‑251 (MQ=255)                                                                                                      atgctctcTGGTCGCGGCCCCGGTCATACTGGCGGTACGCGCGACAAACTGGAATCCACCCCGGGTTTCGAAATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATAAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACGGGCTCCGGGGGATAAACGTTTCTCCG‑CGCCCCGTGATATTCCCGAAACCGTGGAC‑TCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGa                                                                                                                                               <  1:265777/247‑1 (MQ=255)                                                                                                         atCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGCTGCCAACCGCTTCCGA‑GAAATTATTAAAGACGTCGGCGTGGCGCTTATCGGTCAGACCAGCTCCCTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑ATCATCCCGCTGATCAACGCCTCTTTTCTGGCGAAGAAAATGGCGGAAgg                                                                                                                                            >  1:287387/1‑251 (MQ=255)                                                                                                              gggTCGCGCCCTCGTTCTTCCTGGCTGTACGCCCGACAAACGGGAACCCATCCCGGGCTTCAACATTTCCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGGGGCGATTATCGGTCAGACCAGCTCACTGCCTCCGGCGGATAAACGTTTCTCCG‑CGCCCCGTGATATTCCCGCAACCGTGGAC‑CCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTgg                                                                                                                                       <  1:129507/250‑1 (MQ=255)                                                                                                                    cgGCCTCGGTCATACTGGCGGTACGCTCGACAAAATGGAATCCATCCCAGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGGTAAACGTTTTTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGGAGGTCTGGgcgcgc                                                                                                                                 >  1:336807/1‑251 (MQ=255)                                                                                                                       ccTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCCGGCTTCGACATTTTCCC‑GGATGACAACCGGTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGGTTATCGGTTCGACCAGCTCCCTGGCTCCGGCGGGGAAAAGTTTCTAAG‑CGACCCGGGATTTTAACCCAACCGTGGAC‑TCCCTCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGGCGGTCTGGgcgcgcggg                                                                                                                              >  2:289626/1‑248 (MQ=255)                                                                                                                        cTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGAAATTACCGCAACCATGGAC‑TCCATCCCGCTGATCGCCGGCTCAATTCTGGCGAAGAAACTGGCGGAAGGGCTTGACGCGC‑Tggg                                                                                                                            >  1:378869/1‑250 (MQ=255)                                                                                                                                 cctggcggaccgtccaaCAATCTGAAATCCCTCCCTGTCTACGAC‑TTTTCCCTGTATGACAACCGTTCCCGC‑AAAATTATTAAAGACGCTGGCGGGGCTATTATCGGTCAAACCGGCTCATGGGCTCCGGCGGAAAACCTTTTCTACG‑CGCCCCGTTATATTACCCCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGt                                                                                                                   <  2:390100/235‑1 (MQ=255)                                                                                                                                            gCTCGACAAACTGGAATCCATCCCTGGCTTAGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGAAATTACCGCAACCGGGGAC‑TCCATCCCGCTGATCACCGCATCTATTCTGGCGAAAAAACTGGCGGAAGGTCTGGACGCGC‑TGGGGATGGGCGTGGAAGAGGGt                                                                                                         >  1:370487/1‑250 (MQ=255)                                                                                                                                                    aaCTGAAACCCATCCCTGGCTTCGACATTTCCCG‑GGATGACACCCGTTCCCGC‑AAAATTATAAAAGCCGGCGGCGTGGCTATTATCGGCCAGCCCAGCTCACTGGCTCCGGCGGATAAACTTTTCTCCG‑CGCCCCGTGATATTCCCGCACCCGTGGCC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGcgcg                                                                                                <  1:289626/251‑1 (MQ=255)                                                                                                                                                      cTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCCGCTCACTGGCTCCGGCGGGTAAACGTTTCTACG‑CGACCCGTGGTATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGACGAAAATGGAGGAAGGACTGGACGGGCGTGG‑GGTGGACGTGGAAGAGGGGAGCGGGGCGtt                                                                                              >  2:398098/1‑251 (MQ=255)                                                                                                                                                           atccatccCTGACTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAACTTATTAAAGACGTCGGCGTGGCGATTATCGGTAAGACCCGCTCACTGGCTCCGGCGGATAAACGGTTCTACG‑CGAACCGTGATATTACCGCAACCGGGGAA‑GTCATCCCGGTGATACCCGCCTCTATTCTGGCGAAGAACCTGGCGGAGGGTCTGGCCGCGC‑TGGTGATGGGCGTGAAAGgggggcggggagcgtgtaggg                                                                                         >  2:262608/1‑234 (MQ=255)                                                                                                                                                                tccCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCTGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGGGGAC‑TCCATCCCGCTGATAACCGCCTCTTTTCTTGCGAAGAAACTGGCGGAAGGGCTGGACGCGC‑GGGCGGTGGAAGGGGAAGTGGGTAacggggggttttttcaggcc                                                                                    >  2:660192/1‑235 (MQ=255)                                                                                                                                                                    gggttccGACATTTCCCC‑GAATGACAACCGTTTCCGC‑AAAATTATTAAAGCCGCCGGCGTGGCGATTATCGGCCAGACCAGCCCACTGGCCCCGGCGGATAACCGTTTCTACG‑CGCCCCGTGATATTCCCGCACCCTGGGCC‑CCCACCCCGCTGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACg                                                                               <  1:124021/245‑1 (MQ=255)                                                                                                                                                                       cTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAAc                                                                            >  1:611348/1‑251 (MQ=255)                                                                                                                                                                            aCATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCCCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGCAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGCGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAActatct                                                                       >  1:157354/1‑251 (MQ=255)                                                                                                                                                                             cATTTTCCCCGGATAACAACC‑TTTCCCCAGAAATAATTAAAGCGGACGGGGTGGCTATTATGGGTCAGCCCTGCTCACTGGCGCGGGCGGA‑AAACCGTTCTTCGGCGCCCCGTGATATTCCCGCACTCG‑GGACCTCCATCCTGCTGTCCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAActctct                                                                       <  2:256210/251‑1 (MQ=255)                                                                                                                                                                              aTTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTg                                                                      >  1:103083/1‑250 (MQ=255)                                                                                                                                                                                   ccc‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑Tg                                                                                                                              >  1:594691/1‑190 (MQ=255)                                                                                                                                                                                                ccGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGGCGCGC‑TGGCGATGGACGTGAAAGTGGGTAGCGGCGCGGTG‑TTGCCGGCCTACGAGATCTCTGAAACCCCTGCCGGA‑GCgg                                                   >  2:209354/1‑250 (MQ=255)                                                                                                                                                                                                 cGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCCACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGACGGTCTGGTCGCGC‑TGGGGATGGACGTGAAAGGTGGTAGCGGCGGGATT‑AATCCGACCTACGAACTCCCTCAAGCCCTTGCCggagggat                                                   >  2:204290/1‑244 (MQ=255)                                                                                                                                                                                                               ttattaAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGATCACTTGCTCCGGCGGATAAACGGTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCACCTCTATTCTGGCGAAGAAACTGGCGGCAGGGCTGGCCGCGC‑TGGTGATAGACGTGACAGTGGGTAGCCGCGCGGTT‑ATTCCGACCTACGAACTCCCTGTAGCCCTGGCCGGA‑ACGGTTGTTGGCGTCGCt                                     >  2:549798/1‑251 (MQ=255)                                                                                                                                                                                                                    aaaGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGGCCTACGAACTCCCTGAAGCCCTTGCCGCA‑GCGGTTGTTGGCGTGGCCAACgg                                >  2:577704/1‑251 (MQ=255)                                                                                                                                                                                                                     aaGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGAGGAAGGTCTGGACGCGC‑TGGAGATGGACGTGAAAGTGGGTAGCGGCGGGTTT‑ATGCCGACCTACGGACTCTCTGAATCCCTTTCCGGA‑ACGGTTGTTGGCGTGGCTATCGgc                               >  2:150452/1‑251 (MQ=255)                                                                                                                                                                                                                       gACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGGGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTTAT‑CCGGCCCACGAACTCTCTGAAGCCCTTTCCGAA‑GCGATTGTTGGCGTGGCTAACggtcc                             >  2:225109/1‑248 (MQ=255)                                                                                                                                                                                                                        aCGTCGGCGGGGCGGTTATCGGTCAGACCAGCTCACTGGCTCAGGCGGATAAACGTTTCTACG‑AGACCCGTGATATTACCGCAACCGTGGCA‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGAGTGGTCGCGC‑TGGTGATGGTCGTGAAAGTGGGGAGCGGCGCGTTT‑CTGGAGACGCACGAAATCTCTGGAGCCCTTGCCGAAAGC‑ATTGTTGTCGTggtcacatgcgcg                            >  2:6286/1‑240 (MQ=255)                                                                                                                                                                                                                         cggcggcggGGCGTTTATCGGTCAGCCCACCTACCTGGCTCCGGGGGAAAACCTTTTCTACG‑CGACCCGTAATATTACCGCAACCGGGGAC‑TCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTg                           <  1:385968/248‑1 (MQ=255)                                                                                                                                                                                                                               cGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGAGGGC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAGGAAACTAGCGGAGAGTCTGGACGCGC‑TGATGAGGGACGGGAAAGGGAGTAGCGGCGAGTTT‑AGGCCGGCCTAGGAACTCTCTGAAGCCCTTTCCGAA‑GCGATTGTTGGCGTGGCTAACGGAGCTGGCGTgg                     >  1:1241/1‑250 (MQ=255)                                                                                                                                                                                                                                      aTTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGctccac              >  1:676389/1‑251 (MQ=255)                                                                                                                                                                                                                                        tATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGGTTGTTGGCGTGGCTAAAGGCGCTGGGGTGCGCACCCCcg            >  2:214461/1‑251 (MQ=255)                                                                                                                                                                                                                                        tATCGGACAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGAACCGTGATATTACCGCAACCGGGGCC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGAGGAAGGTCTGGAAGCGC‑TGGTGATGGACGTGAAAGTGGGACGCGGCGAGTTT‑CTGCCGACCTCCGAACGATGTGACGCCCTGACTGAG‑GCGATTGTTAGCGTGGCTaatggagtgggagggcgaccacccg            >  2:26380/1‑228 (MQ=255)                                                                                                                                                                                                                                         aTCGGTCAGACCAGCTCCCTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCAGACCTACGAACTCTATGAAGCCCTTGCCGGA‑GCGGTTGTTGGCGTGGCTAACGGCGCTGGCggtcggcccaccgc           >  1:99586/1‑238 (MQ=255)                                                                                                                                                                                                                                            ggTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCGTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGGGGGCACCCCCgcgcc        >  2:37382/1‑250 (MQ=255)                                                                                                                                                                                                                                               cAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCgcggct     >  1:189242/1‑251 (MQ=255)                                                                                                                                                                                                                                                aGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTc    >  1:346346/1‑251 (MQ=255)                                                                                                                                                                                                                                                  aCCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAc  >  1:204937/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           TGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAC  >  NC_000913/4618458‑4618947

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