Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,041,467	A→G	K31K (AAA→AAG)	zinT →	metal‑binding protein ZinT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,041,467	0	A	G	100.0%	156.6 / NA	44	K31K (AAA→AAG)	zinT	metal‑binding protein ZinT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/18);  total (26/18)

CTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTG  >  NC_000913/2041223‑2041705                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               gcggtcttgctggggacgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcccgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcctcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCtt                                                                                                                                                                                                                                           >  2:327219‑M2/157‑250 (MQ=255)        tgctgggtacgggactgggtggagaaacactgtcgaagattgcatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGATGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAgg                                                                                                                                                                                                                                   >  2:339170‑M2/150‑251 (MQ=255)         gctgggtacgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGt                                                                                                                                                                                                                                  <  2:403758‑M2/103‑1 (MQ=255)               tacgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAAc                                                                                                                                                                                                                             >  2:380313‑M2/143‑250 (MQ=255)                  gggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACaaaa                                                                                                                                                                                                                         <  2:419691‑M2/112‑1 (MQ=255)                            ggagaaacacggtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAAt                                                                                                                                                                                                               <  2:16143‑M2/122‑1 (MQ=255)                                              attgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtaacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGGGTCTTTATTGTTAGCGCTACTGCCTTTTCGAATGGGCATCACTCACACGGCAAGCCCTTAACAGAGGTCGCACAAAAAGATGCTAATGGTGTTTTtgatgatgac                                                                                                                                                                                             >  2:228746‑M2/112‑249 (MQ=255)                                                tgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTttgatgatgtca                                                                                                                                                                                            >  2:247259‑M2/110‑247 (MQ=255)                                                                     aattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGATGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAg                                                                                                                                                                      <  1:339170‑M2/163‑1 (MQ=255)                                                                      attgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGt                                                                                                                                                                     <  2:452226‑M2/164‑1 (MQ=255)                                                                           cttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGCGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGCAGAACAAAAAGCTGCTAATGGCGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTg                                                                                                                                                                >  1:158070‑M2/83‑251 (MQ=255)                                                                                 ggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATgg                                                                                                                                                          >  2:15007‑M2/77‑251 (MQ=255)                                                                                              acgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATcc                                                                                                                                             <  1:380313‑M2/188‑1 (MQ=255)                                                                                              acgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaAGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGACATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTGTTTGATGATGCCAATGGACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATcc                                                                                                                                             <  2:379773‑M2/188‑1 (MQ=255)                                                                                               cgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCg                                                                                                                                            <  2:14846‑M2/189‑1 (MQ=255)                                                                                                gttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGt                                                                                                                                           >  1:284160‑M2/62‑251 (MQ=255)                                                                                                     gaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATc                                                                                                                                      >  1:265462‑M2/57‑251 (MQ=255)                                                                                                              tcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGc                                                                                                                             >  2:31887‑M2/48‑251 (MQ=255)                                                                                                                 tcccggcaaaacccgtagttgcgcctgttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGGTTATCCTTTACTGCaa                                                                                                                           <  2:361918‑M2/206‑1 (MQ=255)                                                                                                                 ccccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTATTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAAAAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGAGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCaa                                                                                                                           >  2:401471‑M2/45‑250 (MQ=255)                                                                                                                  cccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCaa                                                                                                                           >  1:274929‑M2/44‑249 (MQ=255)                                                                                                                              cgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTg                                                                                                              >  1:415340‑M2/32‑250 (MQ=255)                                                                                                                                  gttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGAcccc                                                                                                         >  2:38967‑M2/28‑251 (MQ=255)                                                                                                                                     gcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGt                                                                                                       <  1:22205‑M2/226‑1 (MQ=255)                                                                                                                                       gccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCtt                                                                                                    >  1:280236‑M2/23‑251 (MQ=255)                                                                                                                                             tattcaggactcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCaga                                                                                               <  2:266300‑M2/234‑1 (MQ=255)                                                                                                                                               ttcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCaga                                                                                               >  2:191097‑M2/15‑248 (MQ=255)                                                                                                                                                    gacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCg                                                                                       >  1:43787‑M2/10‑251 (MQ=255)                                                                                                                                                            aCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGa                                                                               <  1:16297‑M2/250‑1 (MQ=255)                                                                                                                                                                                   tGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGaa                                                         >  1:214659/1‑250 (MQ=255)                                                                                                                                                                                           tAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCaa                                                                                                                                                                                           >  1:342343/1‑112 (MQ=255)                                                                                                                                                                                           tAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCGTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCaa                                                                                                                                                                                           <  2:342343/112‑1 (MQ=255)                                                                                                                                                                                              gtgtCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGattat                                             >  2:55727/1‑251 (MQ=255)                                                                                                                                                                                                    ttATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATc                                         <  2:222963/249‑1 (MQ=255)                                                                                                                                                                                                        tGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACaa                                     >  2:415461/1‑249 (MQ=255)                                                                                                                                                                                                         gTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAgg                                  >  1:410452/1‑251 (MQ=255)                                                                                                                                                                                                             gcgcTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTAt                              <  2:280236/251‑1 (MQ=255)                                                                                                                                                                                                                  ccTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCa                           <  1:123985/249‑1 (MQ=255)                                                                                                                                                                                                                  ccTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAAc                         >  2:294510/1‑251 (MQ=255)                                                                                                                                                                                                                       cTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGat                     >  2:158007/1‑250 (MQ=255)                                                                                                                                                                                                                                 ggTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatgat          >  1:24220/1‑251 (MQ=255)                                                                                                                                                                                                                                        actcacACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCAtt   <  2:214659/251‑1 (MQ=255)                                                                                                                                                                                                                                          tcacACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTg  >  2:274989/1‑250 (MQ=255)                                                                                                                                                                                                                                           cacaCGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTg  <  1:38967/249‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               CTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTG  >  NC_000913/2041223‑2041705

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