Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,054,126	G→A	L685L (CTG→TTG)	torS ←	sensor histidine kinase TorS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,054,126	0	G	A	100.0%	169.2 / NA	49	L685L (CTG→TTG)	torS	sensor histidine kinase TorS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/17);  total (32/17)

AGCCAATCAAAACCAGCGACGGATATTGCTGTGCCAGTTGTCGGGCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATA  >  NC_000913/1053901‑1054370                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     aGCCAATCAAAACCAGCGACGGATATTGCTGTGCCAGTTGTCGGGCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCAcc                                                                                                                                                                                                                             <  1:179338/251‑1 (MQ=255)   ccAATCAAAACCAGCGACGGATATTGCTGTGCCAGTTGTCGGGCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGc                                                                                                                                                                                                                           >  1:67607/1‑251 (MQ=255)      aTCAAAACCAGCGACGGATATTGCTGTGCCAGTTGTCGGGCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGGGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTg                                                                                                                                                                                                                        >  1:142498/1‑251 (MQ=255)             ccAGCGACTGATATTGCTGTGCCAGTTGTCGGGCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAAAGTAAACCGTAAAGACGCACCGCATGATTGACt                                                                                                                                                                                                                 >  2:311387/1‑251 (MQ=255)                      gATATTGCTGTGCCAGTTGTCGGGCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGc                                                                                                                                                                                                        >  1:182969/1‑251 (MQ=255)                                  ccAGTTGTCGGGCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTg                                                                                                                                                                                            >  1:182501/1‑251 (MQ=255)                                             gCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTaa                                                                                                                                                                                 <  2:219254/251‑1 (MQ=255)                                             gCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTa                                                                                                                                                                                  >  2:299011/1‑250 (MQ=255)                                                   gTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCa                                                                                                                                                                            >  1:298369/1‑250 (MQ=255)                                                     aaTGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGTCAAGc                                                                                                                                                                         >  2:314388/1‑251 (MQ=255)                                                     aaTGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAg                                                                                                                                                                          <  1:114822/250‑1 (MQ=255)                                                        gCCGACGATACCCGGAAGATCAATAGCCTCCAGTGCGGAAGAAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGGAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGt                                                                                                                                                                       <  2:6412/250‑1 (MQ=255)                                                         ccGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTaa                                                                                                                                                                     <  2:225100/251‑1 (MQ=255)                                                           gTCGATATCCGGCAGATCACAATCTACCAGTGCGGCAGCAAACGGTTCGCTATTTTACAGTGTCTCTAAGGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGGGCCGTGGCAACACGTAACGGCAAGCGTAAAc                                                                                                                                                                   >  2:154132/1‑251 (MQ=255)                                                                            cAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAAcc                                                                                                                                                  <  1:55002/251‑1 (MQ=255)                                                                                 tCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGTAAGCGTAAACAAAAACAGCTTCCAACCTCCgg                                                                                                                                             >  2:223365/1‑251 (MQ=255)                                                                                 tCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGAAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCCGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGCCAAGCGTAAACAAAAACAGCTTCCAACCTCCgg                                                                                                                                             >  1:210879/1‑251 (MQ=255)                                                                                    cccaGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGt                                                                                                                                          <  2:393956/250‑1 (MQ=255)                                                                                              catcaAACAGTTCGCTATTTTGCAGTATCTCTAAAGCCTGCGGGGGATTTCCAACAGCAACAATCTGCGCACCACTGGTTTTAAGAATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTACACCGTCAAGACGAACAGCATGATTGACTGATTTTGGCACGGGTACAGTGGCAACAGGTAACGGCAAACGGAGACAAAAACAGCTTCCGACCTCCGACGTGAtggggacg                                                                                                                                >  2:435923/4‑246 (MQ=255)                                                                                              cagcaAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGgcg                                                                                                                                >  1:4110/1‑251 (MQ=255)                                                                                                  aaaCGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCGGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTc                                                                                                                             >  2:264118/1‑250 (MQ=255)                                                                                                     cGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGtt                                                                                                                         >  2:301825/1‑251 (MQ=255)                                                                                                       gTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTcg                                                                                                                       >  2:352089/1‑251 (MQ=255)                                                                                                          cGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCAGTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTcgcc                                                                                                                     <  2:302311/250‑1 (MQ=255)                                                                                                           gCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTcgncgc                                                                                                                   >  2:420470/1‑251 (MQ=255)                                                                                                           gCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTcgccgc                                                                                                                   <  1:283358/251‑1 (MQ=255)                                                                                                                     aGTGTCTGTAAAGACTGCGCGGCATTGCCAACAGCAACAATCTGCGCAAAAATGGTTTTCAGCATCTCAATGGTAACTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGACAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAAACAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCACAGTTCGCCGCCCATCGCCTg                                                                                                         >  2:456558/1‑251 (MQ=255)                                                                                                                     aGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTg                                                                                                         >  1:447203/1‑251 (MQ=255)                                                                                                                            ctAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCCCAGTTCGCCGCCCATCGCCTGTGcc                                                                                                     >  2:108225/1‑248 (MQ=255)                                                                                                                                 gCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGAGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCt                                                                                             >  2:380155/1‑251 (MQ=255)                                                                                                                                        cgGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGAt                                                                                       >  2:248387/1‑250 (MQ=255)                                                                                                                                                     aGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACAGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGAGAGTCAGTCCCAGcc                                                                         >  2:370245/1‑251 (MQ=255)                                                                                                                                                                     ccacTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGttt                                                         <  2:298369/251‑1 (MQ=255)                                                                                                                                                                                 agcacgtcagtGGTAATTCGGTAGGTGAGCGAGTTATCTTCACTAAACAATAGACATAAACCGTCAAGACGCACCGCCTGGTTGACTGCTTTACGCAAGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCACCGCCGTGCTGATGGCGCTCAGTTCGCCGCCCATAGCCTGGCCCAGACGGCGACTGATAGTCAGTCCCAGGCCGGTGACGCCGCGTTTGCCGCTAACCTg                                             >  2:422492/11‑251 (MQ=255)                                                                                                                                                                                        cAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTGCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACaa                                        <  1:435760/249‑1 (MQ=255)                                                                                                                                                                                          aTGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATgg                                    <  2:182501/251‑1 (MQ=255)                                                                                                                                                                                              tAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACATGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTgg                                >  1:235125/1‑251 (MQ=255)                                                                                                                                                                                                aTTCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTAAGTTCGACGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCCGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGa                               >  2:11633/1‑250 (MQ=255)                                                                                                                                                                                                  tCGCTGGGTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAg                             <  2:296781/250‑1 (MQ=255)                                                                                                                                                                                                         gTTAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCATGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGc                     >  1:297804/1‑251 (MQ=255)                                                                                                                                                                                                          ttAGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGcc                    >  2:3607/1‑251 (MQ=255)                                                                                                                                                                                                            aGCGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAg                  <  1:301825/251‑1 (MQ=255)                                                                                                                                                                                                              cGGGTTATCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGGTTTCTGCAAGtt                >  1:284293/1‑251 (MQ=255)                                                                                                                                                                                                               ggggtATCTTCAATTAACAACAAACGTAAACCGGCAAGACGCACCGCCTGATTGACTGGTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGttt               <  1:299850/247‑1 (MQ=255)                                                                                                                                                                                                                     tCTTCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCCAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCGGCCCATCGCCTGGGCCGGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCAGCGTTTGCCGTTTACCTGCACAACTGGCTGGAAGACTTCTGCCAATATCGCggg         >  2:342160/1‑251 (MQ=255)                                                                                                                                                                                                                        tCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGAt       <  1:272727/250‑1 (MQ=255)                                                                                                                                                                                                                        tCAATTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTACCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATc      >  2:58422/1‑251 (MQ=255)                                                                                                                                                                                                                          aaTTAACAACAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCa     >  1:374434/1‑250 (MQ=255)                                                                                                                                                                                                                             tAACAACAAACGTAAACCGTCAAGACGCACCGCCTGAGTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCATAAATGGCTGGAAGATTTCTGCCAATGTCTCGGGATCAATa  <  2:268964/250‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     AGCCAATCAAAACCAGCGACGGATATTGCTGTGCCAGTTGTCGGGCAAGCGTAATGCCGTCGATATCCGGCAGATCAAAATCCACCAGTGCGGCAGCAAACGGTTCGCTATTTTGCAGTGTCTCTAAAGCCTGCGCGGCATTGCCAACAGCAACAATCTGCGCACCACTGGTTTTCAGCATCTCAATGGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATA  >  NC_000913/1053901‑1054370

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