Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,717	T→G	A163A (GCT→GCG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,717	0	T	G	96.1%	145.2 / ‑5.0	51	A163A (GCT→GCG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (32/17);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAC  >  NC_000913/4618472‑4618947                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            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)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            GTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCC‑GGATGACAACCGTTTCCGC‑GAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACG‑CGACCCGTGATATTACCGCAACCGTGGAC‑TCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGC‑TGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTT‑ATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAA‑GCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAC  >  NC_000913/4618472‑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: ‑