Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,804	T→G	L192L (CTT→CTG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,804	0	T	G	96.1%	161.1 / ‑3.3	51	L192L (CTT→CTG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (33/16);  total (34/17)
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.24e-01

CTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGAT  >  NC_000913/4618564‑4619050                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          ctctGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGAGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGGAGAAACTTGCGGAAGGTc                                                                                                                                                                                                                                                 >  1:3722/1‑251 (MQ=255)            gCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGCAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGTTATTCCCGCAACCGTGGACTCCCTCCC‑GCGGGTCACCGCCTCTATTCTGGCGATGAAACTGGCGGGAGGGCTGGAcgcggtg                                                                                                                                                                                                                                      >  1:97849/1‑248 (MQ=255)              cTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGAAGACAACCGTTT‑CCGCGAAATTATTAAAAACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTggg                                                                                                                                                                                                                                    >  1:367642/1‑250 (MQ=255)                 ggTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGGGAt                                                                                                                                                                                                                                 >  1:290559/1‑251 (MQ=255)                    cATACTGGCGGTACGCTCGACAAACTAGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCCGTCAGACCAACTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATAACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGaagg                                                                                                                                                                                                                               >  1:213759/1‑247 (MQ=255)                       aCTGGCGGTACGCTCGACAAACTGGAACCCATCCCGGGCTTCGACATTTCCCCGGATGACAACCGTTC‑CCGCGAAATTATTAAAGACGTCGGTGTGGCGTTTATCGGCCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGCCCCCACCCC‑TCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGt                                                                                                                                                                                                                           <  1:85794/251‑1 (MQ=255)                              gTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTg                                                                                                                                                                                                                    >  1:63378/1‑251 (MQ=255)                                                           ggttccGACATTTTCTCGGATGACACCCGTTT‑CCGCAAATTTAGTAAAGCCGGCGGCGTGGCGATTATCGGTCAGACCGCCCCACTGG‑CCCCGGCGGATAAACGTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGCCTCCTTCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                                                                                                                                                                         <  1:232072/244‑1 (MQ=255)                                                               tcgtcATTTCCCCGGATGACAACCTTTC‑CCGAAAAATTATTAAAGCCGTCGTCGTGGCGATTATCGGTCAGACCAGCTCACTGC‑CTCCGGCGGATAAACTTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAct                                                                                                                                                                                    <  1:78211/246‑1 (MQ=255)                                                                      tttCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTACGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGGGGCGCGTTTATGCCGACCTCCGAAatctccga                                                                                                                                                                              >  1:68503/1‑246 (MQ=255)                                                                         cccGAATGACA‑CCGTTTTCCGCGAAATTATAAAAGGCGCGGGCGGGGCGTTTATCGGCCAGCCCAGCTCACTGT‑CTCCGGCGGAAAACCTTTTCTACGCGCCCCGTGATATTCCCGCAACCGTGGACTCCACCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGccc                                                                                                                                                                        <  1:446112/251‑1 (MQ=255)                                                                                         tt‑CCGCGAAATGATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCAGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGGGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGGCGCGGTGGTGATGGGCGTGCAAGTGGGGAGCGGCGCGGTTATGCCGATCTGCGAACTCTC‑TGAAGGCCATGCCGacacggttgt                                                                                                                                                        >  2:468943/1‑242 (MQ=255)                                                                                             cgcgAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGACGACCTACGAACTCTC‑TGAAGACCCTGCCGA‑AGCGGTTGTTgg                                                                                                                                                    >  1:294957/1‑251 (MQ=255)                                                                                               cgAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTCAGAAGTCTC‑TGAAGCCCTTGCCGA‑GACGCTTGTTGGCg                                                                                                                                                  >  1:134121/1‑251 (MQ=255)                                                                                                aaaaTTAGTAAAGACGGCGGCGTGGCGTTTAGGGGACAGACC‑GGCCCCTGGGCCCCGGCGAATAACCGTTTCTACGCGACCCGTGATATTACCGCAACCGGGGACTCCACCCC‑GCGGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGt                                                                                                                                                 <  1:328498/250‑1 (MQ=255)                                                                                                    ttattaAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGAACCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTCTTCAGGCGAAGAAACTGGCGGAAGGTCTGTACGCGGTGGTGATGGACGTGGAAGGGGGTAGCCGCGCGTTTATGCCGACCTACGCACTCTC‑TGGAGCCCTTGCCGA‑AGCGCTTGTTGGCGTgggt                                                                                                                                             >  1:62730/1‑249 (MQ=255)                                                                                                                     gTGGCGGTTATCGGTCAGACCAGCTCACTGG‑ATCTGGCGGGTAAACGTTTCTATGCGACCCGTGATATTACCGCACCCGTGGACTCCCTCCA‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGCAGGTCTGGACGCGCTGGTGATGGACGAGCAAATGGGTAGCGGCGCGGTTATGCCGACCTATGAAC‑CTCGTGAACCCCTTGCCGAGAGC‑ATTTTTGGCGGGGCTAACGGCTCTggggagcg                                                                                                                            >  2:319973/1‑245 (MQ=255)                                                                                                                       ggCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGGAGGTCTGGACGCGCTGGTGATGGACGTGAAAGGGGGCAGCGGAGCGTTTATGCCGCCCTACGAACTCTC‑TGAAGCCCTTGCCGG‑AGCGATTGTTGGCGGGGCTAACGGCGCTGGCggcggcc                                                                                                                          >  2:427915/1‑245 (MQ=255)                                                                                                                          gATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGc‑acc                                                                                                                       >  1:147039/1‑250 (MQ=255)                                                                                                                           aTTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCTCGTTTATGCCGACCTACGAACTCTC‑TGAAGGCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGGGCTGGCGTGCGc‑accac                                                                                                                     >  1:275808/1‑251 (MQ=255)                                                                                                                              aTCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAACGTGGGTAGCGGCGCGTTTATGCCGACCTACGGACTCTC‑TGAAGCCCTTGCCGG‑AGCGATTGTTGGCGGGGCGAACGGGGCTGGCGTGCGC‑ACCCCcgc                                                                                                                  >  2:187896/1‑251 (MQ=255)                                                                                                                                     aGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATCTTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTATTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTc                                                                                                           >  1:406497/1‑251 (MQ=255)                                                                                                                                       aCCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGGGC‑CCCACCGCGCTGCTCAc                                                                                                         >  1:66729/1‑251 (MQ=255)                                                                                                                                           gctCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTA‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCACGGCTCACCGa                                                                                                      >  1:406831/3‑250 (MQ=255)                                                                                                                                               catggctgccggcgaataatcttttCCACGCGACCCGTGATTTTACCGCCCCCGTGGACCCCATCCC‑GCTTACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATg                                                                                                  <  1:144991/229‑1 (MQ=255)                                                                                                                                                         ggcggATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTgg                                                                                                                                                                                                                   >  1:338676/5‑131 (MQ=255)                                                                                                                                                         ggcgaaaaaaCTTTTCTACGCGCCCCGTGATTTTACCGCACCCGTGGACTCCACCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTAc                                                                                        <  1:208677/244‑1 (MQ=255)                                                                                                                                                           cggATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTg                                                                                      >  1:342225/3‑251 (MQ=255)                                                                                                                                                            ggATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGATGGCGTGCGC‑AACACCGCGATGCTCACCGACCTGAATCAgcatctgg                                                                                     >  1:31186/2‑244 (MQ=255)                                                                                                                                                              aaaacctttctgtccggccccGTGATTTTCCCGCAACCG‑GGAACACATCCCCGCTGACCTCCTCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGcc                                                                                   <  2:407109/233‑1 (MQ=255)                                                                                                                                                                      ttCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGc                                                                           >  1:185996/1‑251 (MQ=255)                                                                                                                                                                         agcgcagcgcGTGATAGTACCGCAACGGGGGACACCCTCCG‑GCTGAGAGCCGCCGCTATGCTGGCGAAGAAACTGGCGGAAGGTCGGGACGCGCGGGTGAGGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGGGGCTAACGGCGCTGGCGTGCGC‑ACCTCCGCGTTGCTTACCGACATGAATCAGGTACTGGCCTCCAGTGCAgg                                                                        <  1:139979/242‑1 (MQ=255)                                                                                                                                                                             cgACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGAGCTGGTGCTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGCAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTCACGGCGCTGGCGTGCGC‑ACCACCGGGCTGCTCAGCGCCCGGAATCTGGTACTGGCCTCCAAGGCAGGTAAc                                                                    >  2:54370/1‑251 (MQ=255)                                                                                                                                                                                   gTGATATTCCCGCACCCGTGGCCTCCACCCC‑GCTGAACCCCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGACGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCGCCAGTGCAGGTAACGCGGtt                                                              <  2:3722/251‑1 (MQ=255)                                                                                                                                                                                                aaacGTGGAATCCATCCC‑GCTGATCACCGCCTCTCTTCTGGCGACGAAACTGGCGGAAGGGCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGATGGCGTGCGC‑ACCACCGCGCTGGTCACCGAAATGAATCAGGTACTGGCCTCCAGTGCCGGTAACGCGGTTGAAGTTCGTGa                                                   >  1:434412/4‑249 (MQ=255)                                                                                                                                                                                                      ggACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAACGGTGc                                           >  1:340312/1‑251 (MQ=255)                                                                                                                                                                                                          cactcccc‑GCTGAGCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGtt                                       <  2:19956/246‑1 (MQ=255)                                                                                                                                                                                                           cctgcccgcctgacccccGCCTCTATTCTGGCGAAGAAACGGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGtt                                       <  2:41495/235‑1 (MQ=255)                                                                                                                                                                                                            cATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGtt                                       >  1:19956/1‑249 (MQ=255)                                                                                                                                                                                                               ccc‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑CCCCCCGCGCTGCTCACCGGCATTAATCAGGGACTGGCCTCCCGTGCAGGTAACGCGGTTGAAGtttgggaagggggggagtttctgg                                  >  2:113511/1‑229 (MQ=255)                                                                                                                                                                                                                   gCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACgg                               >  1:234539/1‑251 (MQ=255)                                                                                                                                                                                                                    cTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGACTGGACGCGCTGGTGATGGACGGGTAAGTGGGTAGCGGCGCGGTTATGCCGACCTACGAACTCCC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑CCCACCGCGCGGCTCACTGACATGCATCGGGTACTGGCATCCCGTGCAGAGAACGCGGTGGCAGTTCGTGagcaggagcggtgcctgaccgg                              >  2:158341/1‑230 (MQ=255)                                                                                                                                                                                                                      gATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTGGCGGAGCGATTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGAGC‑ACCACCGCGCTGCGCAACGACCTGAATCAGGGACTGGCCTACAGCGCAGGGACCGCGGTGGAAGTTCGTGGAGCGGTGCAGGTCCTGGCGGGTg                            >  2:451813/1‑251 (MQ=255)                                                                                                                                                                                                                       aTCACCGCCTCTATTCTGGCGATGAAACTGGCGGAAGGTCTGGACGCGCAGGTGATGGTCGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCTCTGGCGGGCGCCACC‑CCGCGCTGATCCCCGACATGGATCATGTACTGGCCTCCAAGGCAGGGAACGCGGTTGAAGGTCGGGCAGCGGGGCGATTCCTGActggagg                           >  2:502749/1‑245 (MQ=255)                                                                                                                                                                                                                        tagccGCCTCTATTCTGGCGAAGAAAAGGGCGGAAGGTCTGGACGCGCTGGTGATGGACTTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTTCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTGCCTGACGGGTGa                           <  2:341842/247‑1 (MQ=255)                                                                                                                                                                                                                          cccGCCTCTATTCTGGCGAAGAAACTGTCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGagta                        <  2:342225/250‑4 (MQ=255)                                                                                                                                                                                                                          cccGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTACCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGagta                        <  1:362183/250‑4 (MQ=255)                                                                                                                                                                                                                                 ctATTCTGGGGACGAAAGTGGCGGACGGTCTGGACGCGCTGGCGATGGACGTGACAGAGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGACGCCCTTGCCGA‑CGCGATTGTTgg                                                                                                                                                    >  2:209453/1‑121 (MQ=255)                                                                                                                                                                                                                                 ctATTCTGGCGAAGAATCTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCTTCGCGTTTAGGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTgg                                                                                                                                                    <  1:209453/121‑1 (MQ=255)                                                                                                                                                                                                                                     tCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGCTGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGG‑CGCGCTTGTTGGCGTGGCTAACGGAGCTGGCGTGCCC‑ACCACCCCGCCGCACCCCGACAAGAATCAGGTACTGGCCTCCACGGCAGGTTACGCAGGTGAAGATTGTGGAGCGGTGCCGTTTCTgcagggggaggttccgaaccggc             >  2:228064/1‑229 (MQ=255)                                                                                                                                                                                                                                                aaCTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGAt  <  1:178787/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          CTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCC‑GCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGC‑ACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGAT  >  NC_000913/4618564‑4619050

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