Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012971	550,559	C→A	pseudogene (10/364 nt)	ECD_RS02680 →	prophage tail fiber N‑terminal domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012971	550,559	0	C	A	100.0%	193.9 / NA	53	pseudogene (10/364 nt)	ECD_RS02680	prophage tail fiber N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/23);  total (30/23)

CATGGCTGGTAAATGGCTGGTATCCCACAGCAAGTATCAGCGATTTGCTGGTTGTGGTGATGAAGAAAGCCACTGCAGGCATCATGATTAGCTGAATTTTATAACCCAGATACGGGCGCCAGAAATGGCGCCTTTTTTATTTGTGGAGTGAATATGGCAGTACAGATTTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTGGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTTGAGCATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCAGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGT  >  NC_012971/550397‑550842                                                                                                                                                                   |                                                                                                                                                                                                                                                                                            tggccattgaccttgcttatgaaggttccggcagtggtgactggcgcactgacggtttcatcgtgggggtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTAc                                                                                                                                                    >  2:22090‑M1/157‑300 (MQ=255)   gccatttaccttgcgtatgaaggttccggcagtggtgactggcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAANCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAg                                                                                                                                                  >  2:575062‑M1/155‑300 (MQ=255)      attgaccttgcttatgaaggttccggcagtggtgactggcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCAt                                                                                                                                               >  1:79686‑M1/152‑300 (MQ=255)                 ttatgaaggttccggcagtggtgactggcgcacggacggtttcatcgtgggtgtcggttataaattctgattagccangtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACTACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGgatgtcgag                                                                                                                                     >  1:154344‑M1/141‑295 (MQ=255)                          ttccggcagtggtgactggcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCAAGGTACTGGTGAACACGGTAGCATATGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTatggcaagt                                                                                                                           <  2:755412‑M1/169‑5 (MQ=255)                              ggcagtggtgacttgcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccnccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTAc                                                                                                                         >  2:627309‑M1/128‑298 (MQ=255)                              ggcagtggtgacttgcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTAc                                                                                                                         >  2:805674‑M1/128‑298 (MQ=255)                              ggcagtggtgactggcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTAc                                                                                                                         >  2:805424‑M1/128‑298 (MQ=255)                              ggcagtggtgactggcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTAc                                                                                                                         >  1:72106‑M1/128‑298 (MQ=255)                              ggcagtggtgactggcgcactgacgggttcagcgtgggtgtcggtgataaattctgattagccaggtaacacagtgttatgacagcccnccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTAc                                                                                                                         >  2:627299‑M1/128‑298 (MQ=255)                                 agtggtgactggcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGt                                                                                                                    >  2:702314‑M1/125‑300 (MQ=255)                                  gtggtgactggcgcactgacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCAGGGCGTTACAGCATGGATGTCGAGTATGGCAAGTACAGCGTc                                                                                                                   <  2:95022‑M1/177‑1 (MQ=255)                                      tgactggcgcactgactgtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCANTACAGAACTGCACTATTCAACTGAAGGCAAAGNGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCAccc                                                                                                               >  2:616661‑M1/120‑300 (MQ=255)                                               cactgacggtttcatcgtgggtgtcggttataaattctgattagcaaggtaacacagtgttatgacggcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACCGAACTGCACTCTTCAACTGAAGGCAAAGAGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCAACATGCTGGtt                                                                                                       <  2:175014‑M1/189‑1 (MQ=255)                                                   gacggtttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGt                                                                                                                    >  1:542656‑M1/107‑282 (MQ=255)                                                         ttcatcgtgggtgtcggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAANTCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTcc                                                                                            <  1:170990‑M1/200‑1 (MQ=255)                                                                       cggttataaattctgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCAACCTGCTGGTTGAAGGTTTTCCGCCTTCACATgca                                                                               <  2:646520‑M1/213‑2 (MQ=255)                                                                                 ttctgatgatccgggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTTAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAACGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTAc                                                                    >  2:375969‑M1/77‑300 (MQ=255)                                                                                    tgattagccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCg                                                                  >  1:242031‑M1/74‑299 (MQ=255)                                                                                     gattagccaggtaacacagtgttatgacagcccgccgattcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTAACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTc                                                                >  1:665674‑M1/73‑300 (MQ=255)                                                                                          gccaggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGa                                                           <  1:695551‑M1/233‑1 (MQ=255)                                                                                             aggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggggtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGaa                                                          >  1:751072‑M1/65‑298 (MQ=255)                                                                                             aggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTTGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGaa                                                          <  1:319081‑M1/234‑1 (MQ=255)                                                                                             aggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGaa                                                          >  1:649036‑M1/65‑298 (MQ=255)                                                                                             aggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGaa                                                          >  2:44941‑M1/65‑298 (MQ=255)                                                                                             aggtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAgg                                                        >  1:549778‑M1/65‑300 (MQ=255)                                                                                               gtaacacagtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTAAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGACGGCAAAGCGTAACAGCACCAGGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTAACGTATATGAAGGtt                                                      <  2:206416‑M1/238‑1 (MQ=255)                                                                                                      agtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGAcc                                               >  1:89806‑M1/56‑300 (MQ=255)                                                                                                      agtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGAc                                                >  1:101413‑M1/56‑299 (MQ=255)                                                                                                       gtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTAACGTCTATGAAGGTTCCAGacct                                              <  1:99502‑M1/246‑2 (MQ=255)                                                                                                       gtgttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCAGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGACTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGacct                                              <  2:10682‑M1/246‑2 (MQ=255)                                                                                                          ttatgacagcccgccggttcaggcgggcttttttgtggagtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTG‑GGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTa                                          >  1:201264‑M1/52‑300 (MQ=255)                                                                                                                      gaagatatcgtaccggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATg                                 >  1:535885‑M1/40‑298 (MQ=255)                                                                                                                         gatatcgtaccggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGAtttt                            <  1:617557‑M1/264‑1 (MQ=255)                                                                                                                          atatcgtaccggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTtc                           <  1:11122‑M1/265‑1 (MQ=255)                                                                                                                           tatcgtaccggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAAc                                                                          >  1:72330‑M1/35‑252 (MQ=255)                                                                                                                           tatcgtaccggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCATCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAAc                                                                          <  2:72330‑M1/218‑1 (MQ=255)                                                                                                                             tcgtaccggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCGCTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCg                        >  1:417302‑M1/33‑300 (MQ=255)                                                                                                                             tcgtaccggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCg                        >  1:417190‑M1/33‑300 (MQ=255)                                                                                                                               gtaccggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCCGAATGATTTTCTCGGt                      <  1:227137‑M1/270‑1 (MQ=255)                                                                                                                                 acctgtgccgtcagctccgaatatggcaGCATTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCTTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATTTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTCCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGc                    >  2:167036‑M1/29‑300 (MQ=255)                                                                                                                                    ggggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCAt                 <  2:73696‑M1/275‑1 (MQ=255)                                                                                                                                      ggccgtcagctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGa               >  1:483654‑M1/24‑300 (MQ=255)                                                                                                                                              gctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAgat       >  1:793356‑M1/16‑300 (MQ=255)                                                                                                                                              gctccgaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATC‑GGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGAAGGAAgat       <  2:783454‑M1/284‑1 (MQ=255)                                                                                                                                                 gtggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAgatgat    <  2:796282‑M1/288‑1 (MQ=255)                                                                                                                                                   ggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGt  <  2:506929‑M1/290‑1 (MQ=255)                                                                                                                                                   ggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGt  <  1:133670‑M1/290‑1 (MQ=255)                                                                                                                                                   ggatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGAAGGAAGATGATGt  <  1:134478‑M1/290‑1 (MQ=255)                                                                                                                                                   ggatatggcaGCAGTAAAAATCTAAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGCATGATTTTCTCGGTGCCATGACGGAAGATAATGt  <  2:344570‑M1/290‑1 (MQ=255)                                                                                                                                                   gaatatggcaGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTATGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGt  <  2:189752‑M1/290‑1 (MQ=255)                                                                                                                                                   gaatatggcaGCAGTAAAAATCTAAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTATGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGAAGGAAGATGATGt  <  2:186811‑M1/290‑1 (MQ=255)                                                                                                                                                   gaatagggcaGCATTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACGGGATTGTACTATTTAATGGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCCGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGAGGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGt  >  2:699105‑M1/11‑300 (MQ=255)                                                                                                                                                                   |                                                                                                                                                                                                                                                                                            CATGGCTGGTAAATGGCTGGTATCCCACAGCAAGTATCAGCGATTTGCTGGTTGTGGTGATGAAGAAAGCCACTGCAGGCATCATGATTAGCTGAATTTTATAACCCAGATACGGGCGCCAGAAATGGCGCCTTTTTTATTTGTGGAGTGAATATGGCAGTACAGATTTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTGGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTTGAGCATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCAGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGT  >  NC_012971/550397‑550842

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