Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,428,087	C→T	E214K (GAG→AAG)	insH5 ←	IS5 family transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,428,087	0	C	T	100.0%	174.5 / NA	56	E214K (GAG→AAG)	insH5	IS5 family transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/34);  total (22/34)

GGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGG  >  NC_000913/1427797‑1428374                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                  ggCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCAt                                                                                                                                                                                                                                                                                           >  1:490206/1‑299 (MQ=255)     aGGGTGAATAACATCGCCAGTTGGTTATCGTTTTGCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGTGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCTGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGc                                                                                                                                                                                                                                                                                      <  2:229514/300‑1 (MQ=255)        gTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACNTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGa                                                                                                                                                                                                                                                                                   >  1:131156/1‑300 (MQ=255)        gTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACNTAGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGa                                                                                                                                                                                                                                                                                   >  1:132130/1‑300 (MQ=255)        gTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTACATGAAGCAGa                                                                                                                                                                                                                                                                                   >  2:506945/1‑300 (MQ=255)        gTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTACATGAAGCAGa                                                                                                                                                                                                                                                                                   >  1:213528/1‑300 (MQ=255)        gTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTAGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGa                                                                                                                                                                                                                                                                                   >  1:213048/1‑300 (MQ=255)        gTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCATTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGCTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTACATGAAGCAGa                                                                                                                                                                                                                                                                                   >  2:414537/1‑300 (MQ=255)        gTGAATAACATCGCCAGTTGGTTATCGTTTTTAAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTACATGAAGCAGa                                                                                                                                                                                                                                                                                   >  2:90976/1‑300 (MQ=255)              aaCATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTAccc                                                                                                                                                                                                                                                                             <  1:404972/300‑1 (MQ=255)                    gCCATTTGGTTATCGTTTTTCAGCAGCCCCGTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGTTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGGTGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTg                                                                                                                                                                                                                                                                        <  2:574350/299‑1 (MQ=255)                       aGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTTAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTg                                                                                                                                                                                                                                                                    <  1:28501/300‑1 (MQ=255)                           gtttATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTTAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGt                                                                                                                                                                                                                                                                <  2:214636/298‑1 (MQ=255)                           ggTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGt                                                                                                                                                                                                                                                                <  1:84014/300‑1 (MQ=255)                                     tttCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGtt                                                                                                                                                                                                                                                      <  1:11876/300‑1 (MQ=255)                                         agcagcCTCTTGTATCTGGCGTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑TGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGc                                                                                                                                                                                                                                                   <  2:300542/299‑1 (MQ=255)                                         agcagcCCCGTGTATCTGGCTTTCACGAAGCCAAACTGCCACTTGATGATGCGAGACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGGTGGCCGTTTTGTGCTTACGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGc                                                                                                                                                                                                                                                   <  2:300534/299‑1 (MQ=14)                                                ccTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCggt                                                                                                                                                                                                                                             >  1:529391/1‑298 (MQ=255)                                                    gTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTACATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGAc                                                                                                                                                                                                                                       >  2:374634/1‑300 (MQ=255)                                                     tATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGAcc                                                                                                                                                                                                                                      <  1:308237/300‑1 (MQ=255)                                                     tATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGAcc                                                                                                                                                                                                                                      >  1:574507/1‑300 (MQ=255)                                                         tGGCTTTCACTAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCTCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGc                                                                                                                                                                                                                                  <  2:14731/300‑1 (MQ=255)                                                         tGGCTTTCACGAAGCCGAGCTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTTAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGc                                                                                                                                                                                                                                  <  2:267964/300‑1 (MQ=255)                                                         tGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGc                                                                                                                                                                                                                                  <  1:31313/300‑1 (MQ=255)                                                                                   ttGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATc                                                                                                                                                                                                       <  2:155647/299‑1 (MQ=255)                                                                                          gaaacgggtncggCTCCACCCTGGCACGNCGATGCTCGGGCTTTCATGTATTCGATGTTGATTGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGNGTGGGTCAGGCCActc                                                                                                                                                                                                               <  1:505285/280‑1 (MQ=14)                                                                                             gAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAAtgt                                                                                                                                                                                            <  2:352543/300‑1 (MQ=255)                                                                                             gAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAAtg                                                                                                                                                                                             <  1:424957/299‑1 (MQ=255)                                                                                              aaaggggggtccaaCCCTGGCACG‑‑AATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCActc                                                                                                                                                                                                               <  1:202269/267‑1 (MQ=255)                                                                                                  gggTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGc                                                                                                                                                                                        <  1:170266/299‑1 (MQ=255)                                                                                                  gggTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGc                                                                                                                                                                                        <  1:460376/299‑1 (MQ=255)                                                                                                       cTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTTACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCAt                                                                                                                                                                                  >  1:569453/1‑300 (MQ=255)                                                                                                            cccTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCNGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCaa                                                                                                                                                                             >  2:472773/1‑300 (MQ=255)                                                                                                                ggCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGt                                                                                                                                                                          >  1:415272/1‑299 (MQ=255)                                                                                                                        gATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGAC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTATATACATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACt                                                                                                                                                                    >  2:11308/1‑299 (MQ=255)                                                                                                                          gatgctagctcccatatATTCGATGTTGATGGCAGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTCACCTTGAC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTAGGCCAGCTCATCGCACTGTGGCGCTCATTGGTAGCCGGCATCGGCCGAGACAACTTGCTTCTATCCATGAAGCAGATTACCCAGCTGATTGCGATCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAACGTGCCACt                                                                                                                                                                    >  2:11287/17‑299 (MQ=14)                                                                                                                                         gTATTCGATGTTGATGGCCGTTTTATTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCATCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTACATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCGAATGATTGCCTTTCTTGGt                                                                                                                                                    >  2:68968/1‑300 (MQ=255)                                                                                                                                                  gTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCAt                                                                                                                                           >  2:536519/1‑300 (MQ=255)                                                                                                                                                            gTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTTACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGAt                                                                                                                                   <  2:569453/298‑1 (MQ=255)                                                                                                                                                                       gcgcGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑TGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGctc                                                                                                                      <  1:288790/300‑1 (MQ=255)                                                                                                                                                                                             tttACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCTGCCAGCTCCTCGCGCTGTGGCGCTCCTTTGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGTTAGAGCTGGGt                                                                                                <  1:90372/300‑1 (MQ=255)                                                                                                                                                                                                                   ggCGATCAGCCAGTCCACATCCACCTCTGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTTAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCAt                                                                             <  2:516752/299‑1 (MQ=21)                                                                                                                                                                                                                                   acaTCCACCTAGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTACGGATCGCGTTGTTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCtt                                                              >  2:17835/1‑298 (MQ=255)                                                                                                                                                                                                                                     aTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTTAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGt                                                          <  1:357797/300‑1 (MQ=255)                                                                                                                                                                                                                                               gCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTTAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATAATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACg                                                <  1:11287/300‑1 (MQ=255)                                                                                                                                                                                                                                               gCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACg                                                <  1:11308/300‑1 (MQ=39)                                                                                                                                                                                                                                                  aGCTCCTCGCGCTTTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGc                                               <  2:389441/298‑1 (MQ=255)                                                                                                                                                                                                                                                               gtgGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGgccagcc                                <  2:562672/300‑1 (MQ=39)                                                                                                                                                                                                                                                                     cTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGcc                                                                                                                                                                                       >  2:324244/1‑143 (MQ=21)                                                                                                                                                                                                                                                                     cTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCAAGCCAGCGAt                          >  1:45841/1‑300 (MQ=255)                                                                                                                                                                                                                                                                         ttGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGAt                      >  1:152197/1‑300 (MQ=37)                                                                                                                                                                                                                                                                           gttAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGTTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGAttg                    <  2:524532/298‑3 (MQ=255)                                                                                                                                                                                                                                                                             tAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATTAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGTTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTc                  <  1:414537/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             tAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTc                  <  1:90976/300‑1 (MQ=37)                                                                                                                                                                                                                                                                                     aTCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTTAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAAc            <  1:76925/298‑1 (MQ=37)                                                                                                                                                                                                                                                                                             agaCAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGTCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACggg  <  2:490206/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                  GGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACG‑‑GATGCT‑‑GGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCC‑‑GGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCT‑‑GTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGG  >  NC_000913/1427797‑1428374

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