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	97.9%	125.2 / ‑6.4	47	N409S (AAC→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (14/32);  minor base C (1/0);  total (15/32)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTC  >  NC_000913/4619212‑4619699                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         gaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAgg                                                                                                                                                                                                                                                  <  1:172383/251‑1 (MQ=255)            aGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGaaa                                                                                                                                                                                                                                       <  1:264142/251‑1 (MQ=255)             gTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCGCTAAAGACGAAAACAGCTGGCAGGAAGCGGCGaa                                                                                                                                                                                                                                        >  1:307633/1‑249 (MQ=255)                      gACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGaa                                                                                                                                                                                                                                        <  1:273906/240‑1 (MQ=255)                            gATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCCCGCGAAAAACGACAACACCTGGGAGGGAACGGCGAAAGCGGGG‑AAAACGG‑Ca                                                                                                                                                                                                                      >  1:303129/1‑250 (MQ=255)                                       aaGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGcc                                                                                                                                                                                                         <  1:325780/251‑1 (MQ=255)                                       aaGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGACAGCGTCCGCTGGCGGTTATACAAGGGAAAGAAGAAAACAGCGGGCAGGAGGCGGAGAAAACGGGG‑AAAGCGGGCAATTAAA‑‑TTGcc                                                                                                                                                                                                         >  1:246717/1‑251 (MQ=255)                                            gTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCTCGCTGGGGATGGCATTGGGTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCATGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATaa                                                                                                                                                                                                    <  2:398955/251‑1 (MQ=255)                                                 ggcggATACCGAAGGTTGGGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGTCGGCAGGCATCGGACACCATTGATTACATCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCAc                                                                                                                                                                                               <  2:363682/247‑1 (MQ=255)                                                 ggCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCAc                                                                                                                                                                                               <  2:162842/250‑1 (MQ=255)                                                     gATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCg                                                                                                                                                                                             <  1:483270/249‑1 (MQ=255)                                                      aTACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTATCCACGCGAAAGACGGAAACAGCTCGCAAGAAGCGGCGACAGCGGTGGAAAG‑GG‑CAATTAAAGCTT‑CCGTTAAAACACCGGa                                                                                                                                                                                           >  1:242313/1‑250 (MQ=255)                                                       ggcagAAGGTTTTGTCTGTGAAATGGATACCCGCGCGCTGGTGATGGCAGGGGGTGCAATGGACGGCGGTCGTCGTCAGGCATCTGACACCATTGATGACAGCGTCGGCGTTACTGATATGGCGCGTGTGGGAGATCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCATGAAGCGGCGAAAGCGGTG‑AAAGTGG‑CAATTAAA‑CTTGCCGATAAAGCGCCGGaca                                                                                                                                                                                         <  2:428508/247‑3 (MQ=255)                                                        aCCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGaaa                                                                                                                                                                                         <  1:93898/250‑1 (MQ=255)                                                          cGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTTTCCCCCCGAAAGACGGAAACCGCTGGCCGGGAGCGGCGAAAGCGGGG‑GAAGCGG‑CAAATAAA‑CTTGCCGATAAAGCCCCgggaaacc                                                                                                                                                                                      >  1:311984/1‑245 (MQ=255)                                                          cGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCGGGCGGGTATCCCCGCGAAAGAAGAAAACAGCCGGGCGGAAGCGGCGGAAGGGGTG‑AAAGGGG‑GAATAAAA‑CTTTACGATAAAACACCGGaaaaa                                                                                                                                                                                       >  2:425870/1‑248 (MQ=255)                                                               gTTTTGTCAGTGAAATGGATACCCGCGCGTTGGGGATGTCAGTGGGTGCAATGGTCGGCGGGCGGTGGCAGTCATCTGACACCAGTGATTACAGCGTCGGCTTACCTGATACGGCGCGTGTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGCTATTCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCtcacta                                                                                                                                                                                 <  2:133399/251‑3 (MQ=255)                                                                 tttGGCAGGGAAAGGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                                                                                                <  2:388590/250‑1 (MQ=255)                                                                       aGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGc                                                                                                                                                                                                                                                     <  2:471459/178‑1 (MQ=255)                                                                        gTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAAGAAAAAAGCTGGCAGGAAGCGGCGAAAGAGGGG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAACCCACCAACTGTccat                                                                                                                                                                        >  1:80098/1‑248 (MQ=255)                                                                          gAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGG‑GAAGCGG‑CAATTAAA‑CTTTCCGATAAAGCAACGGGAAGCACCCCAACCGTCtttcc                                                                                                                                                                      >  1:242568/1‑247 (MQ=255)                                                                           aaaTGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCCACCGTCTATccc                                                                                                                                                                     >  1:6631/1‑249 (MQ=255)                                                                             aTGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCg                                                                                                                                                                   <  2:451083/251‑1 (MQ=255)                                                                               ggATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTa                                                                                                                                                                 <  1:446711/251‑1 (MQ=255)                                                                                    cccGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCTTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGc                                                                                                                                                            <  2:435103/251‑1 (MQ=255)                                                                                        cgggcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGACAGGCATCTGACACCATTGATTACAGCGTCTGCTTGACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGGCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGGCTATCGCCGTATCAGCGAAt                                                                                                                                                        <  2:90088/248‑1 (MQ=255)                                                                                          cgcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGGCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGcc                                                                                                                                                                    <  2:185159/237‑1 (MQ=255)                                                                                              ggggATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCAGTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCTAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                                                                                                    <  2:459529/249‑1 (MQ=255)                                                                                                  aTGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCCCGGAAGAACAAAACAGCTGGGAGGAAGCGGCGGAAACGGGG‑AAAAGGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAAAAAACCCACTGGCTATCCCCGTATCaccggcaaaagggatac                                                                                                                                              >  1:311959/1‑235 (MQ=255)                                                                                                       aGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatc                                                                                                                                         <  2:6631/251‑1 (MQ=255)                                                                                                         tGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCCCGCGAAAGACGAAAACAGCCGGCCGGAAGCGGCGAAAGCGGGG‑AAAACGG‑CCATTAAA‑CTTGCCCATAAAGCACCGGGAAGCACACCAACCGTCCATCGCCGGATCAGCGGATAACGGGATAccgatctg                                                                                                                                       >  1:388420/1‑251 (MQ=255)                                                                                                                 gtggtgggcggTCGTCGTCCGGCATCTGACACCATTGAATACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAAGAACTTTATACTGATCTGATCATTTa                                                                                                                               <  2:244998/245‑1 (MQ=255)                                                                                                                     gcggtgggcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGTGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTAGCAGCGAATAACGGTATACTGATCTGAGCATTTAAAtt                                                                                                                           <  2:311984/240‑1 (MQ=255)                                                                                                                     gcggcgGTCGGCGGCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAtt                                                                                                                           <  2:388420/251‑1 (MQ=255)                                                                                                                           gtcgtctgcAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGc                                                                                                                     <  2:198413/243‑1 (MQ=255)                                                                                                                                    aGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCCCCGGAAGACGAAAAAAGCTTGCCGGGAGCGGCGAAAGCGGGG‑AAAGGGG‑CGAATAAA‑ATTTCCGATAAAGCACCGGAAAGCACACCAAATGTCTTTCGCCGATTCCGCGAATAACGGGTTACTTATCTGATACTTTTAATTTGTAGCaatgtgaac                                                                                                            >  2:57398/1‑243 (MQ=255)                                                                                                                                    aGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTAc                                                                                                            <  1:57398/251‑1 (MQ=255)                                                                                                                                            gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGGGGACGGAGAAc                                                                                                     >  1:153723/1‑250 (MQ=255)                                                                                                                                            gACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                                                                                                        <  2:6987/247‑1 (MQ=255)                                                                                                                                                     gATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACg                                                                                           <  1:125511/251‑1 (MQ=255)                                                                                                                                                        tACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTg                                                                                         <  1:248331/250‑1 (MQ=255)                                                                                                                                                                 ggCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTtattat                                                                               <  2:153723/251‑1 (MQ=255)                                                                                                                                                                       aCTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATACACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACAACAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAAAGTGCATTTATTATGGTGCt                                                                         >  1:21981/1‑251 (MQ=255)                                                                                                                                                                              gggCGCGTTTGGGCGACCAGTGAGACGGTCAGCGTCCGCTGGCGGATAGCCTCGCGAAAGCCGAAAACAGCTGGCAGGAGGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAGTTGAAGCACTGAGTTCGGAGAACATATGAAACGTGCATTTCTTATGGTGCTGGACTCa                                                                  <  2:69636/250‑1 (MQ=255)                                                                                                                                                                                                                                      aaaGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGacac          >  1:200483/1‑251 (MQ=255)                                                                                                                                                                                                                                           cGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTAGTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTgg     <  2:99819/251‑1 (MQ=255)                                                                                                                                                                                                                                              aaaCAGCTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGGCGGGGCTGACACCCTGGGTc  <  1:286136/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         GAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTG‑AAAGCGG‑CAATTAAA‑CTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTC  >  NC_000913/4619212‑4619699

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