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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,702	0	T	C	91.7%	92.9 / 0.2	37	S158S (AGT→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (18/15);  total (20/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.34e-01

GAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAC  >  NC_000913/4618459‑4618947                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         gAATGGCCCGATTGTTTATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACt                                                                                                                                                                                                                                                     <  2:323926/249‑1 (MQ=255) gAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTgg                                                                                                                                                                                                                                                   <  1:169944/251‑1 (MQ=255)     ggCCCGATTGTCGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACATCCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTc                                                                                                                                                                                                                                               <  1:397288/250‑1 (MQ=255)                       cattCCACTGGTTGCTTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTa                                                                                                                                                                                                                            <  1:9/248‑1 (MQ=255)                       cattCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTa                                                                                                                                                                                                                            <  1:354455/248‑1 (MQ=255)                                      gTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTAAGCGACCCGTGatat                                                                                                                                                                                                             <  1:273074/251‑1 (MQ=255)                                        cGGCGATGGGACTTCGCTGATGTTTGGTCCGATGGTCGCATCCTGCGGCGGCTTTATTCCGATGATCTGTGGTCGCGGTCTCGGTCATACTGGCTGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGGCAGTTTCCCGGATGACAACCGTTT‑CCTCGAAATTAGTAAAGACGTCGGGGTGGCGATTATCGGTCAGACCATCTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGCCCCCTGatatc                                                                                                                                                                                                            <  2:323521/250‑2 (MQ=255)                                                 gACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAAcc                                                                                                                                                                                                   <  1:142713/250‑1 (MQ=255)                                                 gACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCg                                                                                                                                                                                                  <  1:252259/251‑1 (MQ=255)                                                             gTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATc                                                                                                                                                                                       <  1:8671/250‑1 (MQ=255)                                                                     cGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATc                                                                                                                                                                              <  1:349293/251‑1 (MQ=255)                                                                                           cTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAAAGTTTCTACGCGACCCGTGTTATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGGCTCTATTTTGGCga                                                                                                                                                           >  1:425267/1‑248 (MQ=255)                                                                                                          ctctGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGAGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGGAGAAACTTGCGGAAGGTc                                                                                                                                         >  1:3722/1‑251 (MQ=255)                                                                                                                     gCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGCAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGTTATTCCCGCAACCGTGGACTCCCTCCCGCGGGTCACCGCCTCTATTCTGGCGATGAAACTGGCGGGAGGGCTGGAcgcggtg                                                                                                                              >  1:97849/1‑248 (MQ=255)                                                                                                                       cTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGAAGACAACCGTTT‑CCGCGAAATTATTAAAAACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTggg                                                                                                                            >  1:367642/1‑250 (MQ=255)                                                                                                                          ggTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGGGAt                                                                                                                         >  1:290559/1‑251 (MQ=255)                                                                                                                             cATACTGGCGGTACGCTCGACAAACTAGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCCGTCAGACCAACTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATAACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGaagg                                                                                                                       >  1:213759/1‑247 (MQ=255)                                                                                                                                aCTGGCGGTACGCTCGACAAACTGGAACCCATCCCGGGCTTCGACATTTCCCCGGATGACAACCGTTC‑CCGCGAAATTATTAAAGACGTCGGTGTGGCGTTTATCGGCCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGCCCCCACCCCTCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGt                                                                                                                   <  1:85794/251‑1 (MQ=255)                                                                                                                                       gTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTg                                                                                                            >  1:63378/1‑251 (MQ=255)                                                                                                                                           gCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCAcc                                                                                                                                                                           >  1:66024/1‑184 (MQ=255)                                                                                                                                           gCTCGACAAACTGGAATCCAACCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGGCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCAGCCCGCTGATCAcc                                                                                                                                                                           <  2:66024/184‑1 (MQ=255)                                                                                                                                                                    ggttccGACATTTTCTCGGATGACACCCGTTT‑CCGCAAATTTAGTAAAGCCGGCGGCGTGGCGATTATCGGTCAGACCGCCCCACTGG‑CCCCGGCGGATAAACGTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGCCTCCTTCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                                                                 <  1:232072/244‑1 (MQ=255)                                                                                                                                                                        tcgtcATTTCCCCGGATGACAACCTTTC‑CCGAAAAATTATTAAAGCCGTCGTCGTGGCGATTATCGGTCAGACCAGCTCACTGC‑CTCCGGCGGATAAACTTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAct                                                                            <  1:78211/246‑1 (MQ=255)                                                                                                                                                                               tttCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTACGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGGGGCGCGTTTATGCCGACCTCCGAAatctccga                                                                      >  1:68503/1‑246 (MQ=255)                                                                                                                                                                                  cccGAATGACA‑CCGTTTTCCGCGAAATTATAAAAGGCGCGGGCGGGGCGTTTATCGGCCAGCCCAGCTCACTGT‑CTCCGGCGGAAAACCTTTTCTACGCGCCCCGTGATATTCCCGCAACCGTGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGccc                                                                <  1:446112/251‑1 (MQ=255)                                                                                                                                                                                                  tt‑CCGCGAAATGATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCAGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGGGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGGCGCGGTGGTGATGGGCGTGCAAGTGGGGAGCGGCGCGGTTATGCCGATCTGCGAACTCTC‑TGAAGGCCATGCCGacacggttgt                                                >  2:468943/1‑242 (MQ=255)                                                                                                                                                                                                      cgcgAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGACGACCTACGAACTCTC‑TGAAGACCCTGCCGA‑AGCGGTTGTTgg                                            >  1:294957/1‑251 (MQ=255)                                                                                                                                                                                                        cgAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTCAGAAGTCTC‑TGAAGCCCTTGCCGA‑GACGCTTGTTGGCg                                          >  1:134121/1‑251 (MQ=255)                                                                                                                                                                                                         aaaaTTAGTAAAGACGGCGGCGTGGCGTTTAGGGGACAGACC‑GGCCCCTGGGCCCCGGCGAATAACCGTTTCTACGCGACCCGTGATATTACCGCAACCGGGGACTCCACCCCGCGGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGt                                         <  1:328498/250‑1 (MQ=255)                                                                                                                                                                                                             ttattaAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGAACCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTCTTCAGGCGAAGAAACTGGCGGAAGGTCTGTACGCGGTGGTGATGGACGTGGAAGGGGGTAGCCGCGCGTTTATGCCGACCTACGCACTCTC‑TGGAGCCCTTGCCGA‑AGCGCTTGTTGGCGTgggt                                     >  1:62730/1‑249 (MQ=255)                                                                                                                                                                                                                              gTGGCGGTTATCGGTCAGACCAGCTCACTGG‑ATCTGGCGGGTAAACGTTTCTATGCGACCCGTGATATTACCGCACCCGTGGACTCCCTCCAGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGCAGGTCTGGACGCGCTGGTGATGGACGAGCAAATGGGTAGCGGCGCGGTTATGCCGACCTATGAAC‑CTCGTGAACCCCTTGCCGAGAGC‑ATTTTTGGCGGGGCTAACGGCTCTggggagcg                    >  2:319973/1‑245 (MQ=255)                                                                                                                                                                                                                                ggCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGGAGGTCTGGACGCGCTGGTGATGGACGTGAAAGGGGGCAGCGGAGCGTTTATGCCGCCCTACGAACTCTC‑TGAAGCCCTTGCCGG‑AGCGATTGTTGGCGGGGCTAACGGCGCTGGCggcggcc                  >  2:427915/1‑245 (MQ=255)                                                                                                                                                                                                                                   gATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGcacc                >  1:147039/1‑250 (MQ=255)                                                                                                                                                                                                                                    aTTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCTCGTTTATGCCGACCTACGAACTCTC‑TGAAGGCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGGGCTGGCGTGCGcaccac              >  1:275808/1‑251 (MQ=255)                                                                                                                                                                                                                                       aTCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAACGTGGGTAGCGGCGCGTTTATGCCGACCTACGGACTCTC‑TGAAGCCCTTGCCGG‑AGCGATTGTTGGCGGGGCGAACGGGGCTGGCGTGCGCACCCCcgc           >  2:187896/1‑251 (MQ=255)                                                                                                                                                                                                                                              aGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATCTTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTATTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTc    >  1:406497/1‑251 (MQ=255)                                                                                                                                                                                                                                                aCCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGGGCCCCACCGCGCTGCTCAc  >  1:66729/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         GAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAC  >  NC_000913/4618459‑4618947

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