Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,454	A→G	N409S (AAC→AGC)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,454	0	A	G	98.1%	137.3 / NA	53	N409S (AAC→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (15/37);  minor base T (0/1);  total (15/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

GTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACG‑CCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTC  >  NC_000913/4619207‑4619699                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        gTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTg                                                                                                                                                                                                                                                      <  1:338281/251‑1 (MQ=255)      gacagacgagAAGTATCTGCCGACAGGGATGCTGAGTAAAGCAGATC‑AGTCTGATACCGAAGGGTTTGTCAGTGGAATGGATTTCCGCGGGCTGGGGAGGGCTGTGGTTGCCAGGGGCGGCGGTCG‑TCGTCAGGCAGCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCGGAGTCCGCTGGCGGTTAT‑CCACGCGAAAGCCGAAAACAGCTGGCAgg                                                                                                                                                                                                                                                 <  2:157354/242‑1 (MQ=255)        aCTACGCGAAGTATCTGCCGGCAGCGATGCTGACGAAAGCAG‑TCTATGCTGATTCCGAAGGTTTTGTCAGGGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGaa                                                                                                                                                                                                                                               <  2:622657/251‑1 (MQ=255)                aaGTATCTGCCGATAGCGATGCGGACGAAAGCCG‑TCTATGCAGAGACCGAAGGTTGTGTCGGTGAAATGGATCCCCTCGCGCTTGGGATGGCAGTGGTTGCAATGGGCGGCGGGCG‑TCTTCAGGCATCTGACTGCATTGATTACAGCGTCCGCTTTACTGATAGGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCTCTGGCGGTTAT‑CCACGCGAAAGACGAAAACATCTGGAAGTAAGCGGCGaa                                                                                                                                                                                                                                       <  2:57357/251‑1 (MQ=255)                    aTCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑TCCCGCGAAAAACGGAAAAAGCGGGCAGGAAGCGGGGAAAGCg                                                                                                                                                                                                                                   >  1:678698/1‑251 (MQ=255)                       tGCCGACAGCGATGCTGACGAAAGCAG‑TTTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGt                                                                                                                                                                                                                                 <  2:485386/250‑1 (MQ=255)                        tccGACAGCGATGCTGACTAAAGCAG‑GCTATGCTGATACCGAAGGTGTTGTCAGTGAAAGGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGTCAGGAAGCGGCGAAAGCGGt                                                                                                                                                                                                                                 <  2:189242/248‑1 (MQ=255)                            aCAGCGATGCTGACGAAAGCCG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGAAGCCGGTCAGCGTCCGCTGGCGGGTTT‑CCACGCGAAAGACGGAAACAGCCGGCCGGAAGCAGCGGAAGCGGAGaaga                                                                                                                                                                                                                            >  2:39956/1‑248 (MQ=255)                                          gAAAGCAG‑TCTATGCTTATACCGAAGGGTTTGTCAGTGAAATGGATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTTAGCGTCCGCTGGCGTTTAT‑CCACGCGAAAGACTAAAACAGCTGTCAGGAAgcgg                                                                                                                                                                                                                                           <  2:571573/221‑1 (MQ=255)                                             aGCAG‑TCTATGCTGATCCCGAATGGTTTGTCAGTGAAATGGGTACCCGCGCGCTGTGGATGGCAGGGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCg                                                                                                                                                                                                         <  2:252251/251‑1 (MQ=255)                                                                       tttGTCAGTGAAATGGATACCCGCGCGCGGGGGATGGCAGTGGTTGCAATGGGCGGCTGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                                                                                                 <  2:48413/250‑1 (MQ=255)                                                                         tgtacgTGAAATGGAGACCCGCGCGCTGGGGATGGTAGTGGTTGCAATGGGCGGCTGTTGTATGT‑AGGCATCTGACACCATTGATTACAGCGGCGGCGTTACTGATTTGGCGCGTTTGGGCGACCAGTTAGACGGGCAGCGTCCGATGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑TAATTAAACTTGCCGATAAAGCACCTGAAAGCACATCaaat                                                                                                                                                                               <  2:473201/245‑3 (MQ=255)                                                                            cagggAAATGGATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCTGTCT‑TCGTCAGGCATCTGAAACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTTCCGATAAAGCACCGGAAAGCACCCCAACTGTCt                                                                                                                                                                           <  2:649014/247‑1 (MQ=255)                                                                                                cgcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                                                                                                       <  1:530687/251‑1 (MQ=255)                                                                                                 gcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAAAAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGG‑CCATTAAACTTGCCGATAAAGCACCGGGAAACACACCAACAGTCTATTGCCGTATCAGCGAATaa                                                                                                                                                       >  1:338733/1‑250 (MQ=255)                                                                                                       gaaggCATTGGTTGCAATGGGCGGTGCTCT‑TCGTCAGGCATCAGACCCCATTGATTACAGCGTCGGCTTTACGGATATTGCGCGTTTGGGCGACCAGGTAGACGGTTAGCGTCCGCTGGCGGTTAT‑CCACTCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAATCACACCAATTGTCTATCGCCGTATCAGCGAATAACGGtata                                                                                                                                                <  2:117209/248‑1 (MQ=255)                                                                                                          ggCAGTGGTTGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActg                                                                                                                                             <  1:137083/251‑1 (MQ=255)                                                                                                              gTGGTTGCAATGGTCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCTAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCATCTGTCTATCGCCGTATCAGCGAATAACGGTATActgatct                                                                                                                                         <  2:573779/251‑1 (MQ=255)                                                                                                                ggTTACAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGTGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCTCCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctga                                                                                                                                       <  1:70123/251‑1 (MQ=255)                                                                                                                 ggtGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctgat                                                                                                                                      <  1:33388/249‑1 (MQ=255)                                                                                                                   tGCAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACtgatctgatc                                                                                                                                     <  1:190777/250‑1 (MQ=255)                                                                                                                     cAATGGGCGGCGGTCG‑TCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTAT‑CCCCGCGAAAGGCGAAAACAGCTGGCCGGAAGCGGCGGAAGCAGTGAAAGCGGGCAAT‑AAAAGTGCCGGTCAAGCACCGGGAAGCACACCCACTGTCTATCCCCGGATCAGCGAATAACGGGAAACTGATATGATCAtt                                                                                                                                  >  1:654888/1‑251 (MQ=255)                                                                                                                         gggCGGCGGTTG‑TCGTCAGGCATCTGACACAATTGAGTACAGCGTCGGTATTACTGGTATGGCTCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTTATCTGATCATTTaa                                                                                                                               <  2:665587/250‑1 (MQ=255)                                                                                                                               gcggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGa                                                                                                                         <  1:101410/243‑1 (MQ=255)                                                                                                                                 ggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGGATCAGCGAATAACGGGATACTGATCTGATCATTTAAATTTGa                                                                                                                         >  1:417007/7‑249 (MQ=255)                                                                                                                                             gCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCacac                                                                                                                                                                                    <  2:482602/178‑1 (MQ=255)                                                                                                                                             gCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAACGG‑CAATTAAAATTGCCGATAAAGCACCGGGAAGCacac                                                                                                                                                                                    >  1:482602/1‑178 (MQ=255)                                                                                                                                                   gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACa                                                                                                    <  2:289460/251‑1 (MQ=255)                                                                                                                                                    acacCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAAc                                                                                                                                                      <  1:102069/201‑1 (MQ=255)                                                                                                                                                       ccATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATg                                                                                                <  1:24847/251‑1 (MQ=255)                                                                                                                                                           tgATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCCGGGAGCGGCGAAAGCGGGGAAAGGGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGAACACCAACTGTCTATACCCGAATCAGCGAATAACGGGATACTGATCTGATCATTTTAAATTGAAACACTGAATACGGAGA‑ACATATggaa                                                                                             >  1:42788/2‑247 (MQ=255)                                                                                                                                                            gATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACg                                                                                           <  1:122481/251‑1 (MQ=255)                                                                                                                                                                aCAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGG‑CAATTAAAATTGCCGATAAAGCACCGGAAAGCACACCCACTGTCTATCGCCGTATCCGCGAATAACGGGATACTGATCTGATCATTTAAATTTTAAGCACTGGGGACGGAGA‑AAATATGAAACGTGc                                                                                        >  1:487037/1‑250 (MQ=255)                                                                                                                                                                   gCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCAttt                                                                                    <  1:39956/251‑1 (MQ=255)                                                                                                                                                                                    atGGCGCGCTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTCATACCA‑GCGAAAGACGAAAACAGCACGGAGGCAGCGAAGAAAGCGGTGAAAACGG‑CAATTAAAAATGCCGCTAAAAGACCGGAGAGAAACCCAACTGGCTATACCCGAATCAGCGAATCAAGGGATACCGTTATGATAATTTAACTTTTGAGCACAGAGTACGGAGACACA‑ATGAAACGTGCATttattttagagcttgactc                                                                   >  1:46902/1‑237 (MQ=255)                                                                                                                                                                                               gggCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGAGCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCg                                                        <  1:427501/251‑1 (MQ=255)                                                                                                                                                                                                 gCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGgc                                                      <  2:487037/251‑1 (MQ=255)                                                                                                                                                                                                  cGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGGGA‑ACATATGAAACGTGCATTTATTATGGGGCTGGGCTCATTCGGCATCGgcg                                                     >  1:489074/1‑251 (MQ=255)                                                                                                                                                                                                     ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTa                                                  <  1:259137/251‑1 (MQ=255)                                                                                                                                                                                                        ggTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACag                                               <  1:261079/251‑1 (MQ=255)                                                                                                                                                                                                           aGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACagaag                                            <  1:327934/251‑1 (MQ=255)                                                                                                                                                                                                                      gTCCGCTGGCGGTTAT‑CCACGCGCAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAAATTGCCGGTAAAGCAACGGAAAGCACACCAACTGTTTATCGCCGTAACAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTAAGGAGA‑ACATATGAAACGTGCATTTAATATTGTGCTGGACTCATTCGGCAGCGGCGCTACCGAAGATGAAGAACGc                                 >  2:521141/1‑251 (MQ=255)                                                                                                                                                                                                                            gggCGGTTAT‑CCACGCTAAAGCCGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGATATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCTGCGCTACAGAAGATGCAGAACGCTTTGGt                           <  2:42788/250‑1 (MQ=255)                                                                                                                                                                                                                                 gTTAT‑CCACGCGAAAGACGCAAACAGCTGGCAGGAAGCGGCGAAAGCGGAGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGAATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAAAGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAAAGCTTTGGTGACGt                      >  2:515292/1‑251 (MQ=255)                                                                                                                                                                                                                                       ccACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCggg                  <  1:359049/250‑1 (MQ=255)                                                                                                                                                                                                                                             gAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGaca           <  1:77552/251‑1 (MQ=255)                                                                                                                                                                                                                                              aaaGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGacac          <  1:487495/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aCGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGCCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGGACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTg      <  2:489074/251‑1 (MQ=255)                                                                                                                                                                                                                                                   cGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTgg     <  1:74686/251‑1 (MQ=255)                                                                                                                                                                                                                                                   cGAAAACAGCGGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTAGCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATGCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTgg     <  2:678698/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gcaaaCAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCTCCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTAAAGACGATGCAGACAGCTTTGGTGACGTCGGGGCTGACACACTggg    >  2:149453/3‑251 (MQ=255)                                                                                                                                                                                                                                                     aaaaCAGCTGGCAGAAAGCGGCGAAAGCGGTGAAAGCGG‑CCATTAAACTTGCCGCTACAGCAACGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTCTTATGGTGCTGGACTAATTCGGCATCGGCGCTACAGAAGATGCAGAACGATTTGGTGACGTCTGGGCTGACACCCTgggg   >  2:317203/1‑250 (MQ=255)                                                                                                                                                                                                                                                      aaaCAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACAACAACTGTCTATCGCCGTATCAGCGAATAACGGAATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGAATTTATTATGGTGCTGGCCTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCggggcc  >  2:94963/1‑249 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        GTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAG‑TCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACG‑CCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTAT‑CCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGG‑CAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGA‑ACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTC  >  NC_000913/4619207‑4619699

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