Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,142	T→C	intergenic (+22/‑36)	creC → / → creD	sensory histidine kinase CreC/putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,142	0	T	C	100.0%	195.6 / NA	63	intergenic (+22/‑36)	creC/creD	sensory histidine kinase CreC/putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/38);  total (25/38)

ACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATA  >  NC_000913/4637896‑4638382                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 gcccTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTc                                                                                                                                                                                                                                              <  2:329851/248‑1 (MQ=255)      tAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTAATTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGGCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTa                                                                                                                                                                                                                                           <  2:238888/249‑1 (MQ=255)       aaGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGAATcc                                                                                                                                                                                                                                        >  2:315329/1‑251 (MQ=255)        aggtgcgtgATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACGCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTTAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCct                                                                                                                                                                                                                                       <  1:610181/243‑1 (MQ=255)        aggagCTGGATTCCGGTAGTGGTATTCCTGGCTACGCGCTGTCACGTATTTTTGAACGCGTTTCCTCTTTTCCGCGTGCAAAGGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGGGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTACACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCct                                                                                                                                                                                                                                       <  2:391613/247‑1 (MQ=255)           tGCTGGATACCGGTAGTGGAATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGAGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCctgct                                                                                                                                                                                                                                    >  2:53258/1‑251 (MQ=255)            gCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCctgctg                                                                                                                                                                                                                                   <  1:304411/251‑1 (MQ=255)                 aTACCGTTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCGTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCAtt                                                                                                                                                                                                                              <  1:37643/251‑1 (MQ=255)                    ccGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTTGCGTGCTGGCCTCGCTTCGACTTCGCCGTCACTTCACATAGCTGCAAATTCGTCCCACATAGCCGTCGTATCCTGCTGCCATTGCa                                                                                                                                                                                                                           <  2:602001/251‑1 (MQ=255)                     cGGTAGTGGTTTTCCTGACTACTCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCaa                                                                                                                                                                                                                          <  1:447449/251‑1 (MQ=255)                                      aCTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTg                                                                                                                                                                                                         <  1:511761/251‑1 (MQ=255)                                          cgcgCTGTCACTGATTTCTGAACGCGTGTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTGCGTCAGTGAGTTCGCCCGTTTGTGTAACGACGAAGTCACGCTGCGCAACGTGCAGTAAGTTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                                                                                                                     <  1:79741/251‑1 (MQ=255)                                          cgcgCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTTCCAATGGGCCAAAACGCAGCGGTCTGGTTTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTCACGTCGAAGTCACTCTGCGCCACGGGCAGTAAGGTGGCGTGCTGGCCTCGCTTCGACTTCCCCGTCACTTCACATAGCTTCAAATTCTTCCCCCATAACCCTCGTATCCTGCTGCCATTGCAAAGGCGAAGCCTATTTTCAAAt                                                                                                                                                                                                     >  2:533357/1‑251 (MQ=255)                                             gctgtCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAATGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccc                                                                                                                                                                                                  >  2:152267/5‑251 (MQ=255)                                                gtaaCGTATTTTTGAACGGTTTGACTCTTTGCCGCGGGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACGTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATcccc                                                                                                                                                                                                 <  2:105271/246‑1 (MQ=255)                                                gtCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccccc                                                                                                                                                                                                <  1:91964/249‑1 (MQ=255)                                                   aCGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGt                                                                                                                                                                                            >  1:488695/1‑251 (MQ=255)                                                   aCGTATTTTTGAACGCTTGTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGt                                                                                                                                                                                            <  1:140084/251‑1 (MQ=255)                                                     gTATTTTTGAACGCTTTTACGCTTTGCAGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGTGTTGTCGTTCGTCAGTGAGGTCGCCCGTTTGTGTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATATCGTCAAATTCTTCCCACATAGCCTACGTATCCTGCTGCCATTGCAAAGTAGAAGAGTATGTTGAATTCCCCCCTGTTc                                                                                                                                                                                          <  2:330188/251‑1 (MQ=255)                                                       aTTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCATGTTCt                                                                                                                                                                                         <  1:123140/250‑1 (MQ=255)                                                           ttGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGaa                                                                                                                                                                                     <  1:262077/250‑1 (MQ=255)                                                             ggaCGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTGCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGaaaa                                                                                                                                                                                   <  1:379221/248‑1 (MQ=255)                                                                   ttttACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGAccac                                                                                                                                                                            >  2:220783/1‑248 (MQ=255)                                                                          ctTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGAccaccttgttt                                                                                                                                                                     >  1:161447/1‑241 (MQ=255)                                                                              gccgcGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGt                                                                                                                                                                  <  2:458573/246‑1 (MQ=255)                                                                              gccgcGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTg                                                                                                                                                                 <  1:474262/247‑1 (MQ=255)                                                                              gccgcGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCGTCGTATCCTGCTGGCATTGCAAAGGAGACGACTATGTTGAAATCCCCCCTGTTCTGGAAAATgctccccttgtttggtg                                                                                                                                                                 <  2:280968/247‑17 (MQ=255)                                                                                           gggCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGGCACTTCACATAGCTTCAAATTCTTCCCACATAGCCGTCGTATCCTGCTGCCATTGCAAAGGAGAAGTCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCtgt                                                                                                                                                     <  2:261448/250‑1 (MQ=255)                                                                                              cAAAAAAGCAGCGGTCTGGGGTTGGTGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGTAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTgttgtt                                                                                                                                                 <  2:265195/251‑1 (MQ=255)                                                                                              cAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTgttgtt                                                                                                                                                 <  1:466440/251‑1 (MQ=255)                                                                                                  aaaacagcGGTCTGGGGTTGGAGTTCGTCAGTGAGGTCGCACGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGATGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGCAACTGACCACCTTGTTTGGTGCAGTATTGCTGttgttgttt                                                                                                                                             >  2:380060/5‑248 (MQ=255)                                                                                                  aaaGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGTAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                                                                                                                             >  1:350049/1‑251 (MQ=255)                                                                                                    agcagcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                                                                                                                           >  1:496634/1‑251 (MQ=255)                                                                                                    agcagcGGTCTGGGGGTGGCGTTCGGCAGGGAGGTCGCCCGTTGGTTGGACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGTTTCGACTTCACCGTCACTTCACAGAGCGGCAAATTCTACCCACATAGCCTTCACACCCTGCTGCCATTGCAAAGGACAAGACTACGTTGAAATCCCCCCTTTTCTTGGAAATTACCACCGTGTGTGGGGCAGTATGTCTGTTGTGGATTcc                                                                                                                                           <  2:254122/251‑1 (MQ=255)                                                                                                        gcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTGGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                                                                                                                       <  1:397785/251‑1 (MQ=255)                                                                                                           gTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCACGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGACAATGACCACCTTGTTTGTAGCAGTATAGCTGTTGTTGATTCCGATactg                                                                                                                                    >  1:525228/1‑248 (MQ=255)                                                                                                              tGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                                                                                                                                 >  2:525646/1‑251 (MQ=255)                                                                                                               ggggTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGTCCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                                                                                                                                 <  1:241288/250‑1 (MQ=255)                                                                                                                   atGGCGTTCGTCAGAGAGGCAGACCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGAGGCGTGCTGGCCTCGCTTCGACTTCACCGTAACTTCACATAGATTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCGGATAATGCTGATTCg                                                                                                                            >  2:76957/2‑251 (MQ=255)                                                                                                                      gCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGACGGGGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAATTTCTTCCCACATAGCCTTCGTAGCCTGCTGCCATTGCAAAGGAGAAGCCTATGTTGAAATCCTCCCTGTTCTGGAAAATGACCACCTTGTTGGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                                                                                                                         <  2:607296/251‑1 (MQ=255)                                                                                                                              cAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGAtt                                                                                                                  <  1:593237/250‑1 (MQ=255)                                                                                                                                            gtttgtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGAt                                                                                                   <  1:451284/251‑1 (MQ=255)                                                                                                                                                gtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTAcc                                                                                               <  1:57617/251‑1 (MQ=255)                                                                                                                                                 tttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTAt                                                                                                                                                           <  1:475317/190‑1 (MQ=255)                                                                                                                                                 tttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTAt                                                                                                                                                           >  2:475317/1‑190 (MQ=255)                                                                                                                                                                 cGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGa                                                                                  >  2:173608/1‑247 (MQ=255)                                                                                                                                                                  gCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATg                                                                             >  2:88778/1‑251 (MQ=255)                                                                                                                                                                   cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGc                                                                            >  2:239448/1‑251 (MQ=255)                                                                                                                                                                    tGCGCAACGCGCAGGCAGGTGGCGTGCTGGCATCGCTTCGCCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGACAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGGTTACCGGAGCGATGCGGCAGATGCg                                                                           >  2:626610/1‑251 (MQ=255)                                                                                                                                                                       gcAACGTGCAGGAATGTGGCGTGCTGGCCTCGCTTCGTCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTTTTCTGTAAAATGACCTCCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGAtt                                                                        <  1:226168/251‑1 (MQ=255)                                                                                                                                                                            gTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGcc                                                                    >  2:302400/1‑250 (MQ=255)                                                                                                                                                                                      gTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATTACCACCTTGTTTGGTGCAGGATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGc                                                         <  1:626610/251‑1 (MQ=255)                                                                                                                                                                                                                           cacaTAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTc                    >  1:82772/1‑251 (MQ=255)                                                                                                                                                                                                                               tAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAACCTCGTTGGGCCGCTTATCGCTATTCCtgt                >  2:107097/1‑251 (MQ=255)                                                                                                                                                                                                                                  cTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGAc             >  2:126207/1‑251 (MQ=255)                                                                                                                                                                                                                                    tCAAATTCTTCCCACATAGCCTTCTTATCCTGCGGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGAcc            <  2:200301/250‑1 (MQ=255)                                                                                                                                                                                                                                     cAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGa          <  1:53258/251‑1 (MQ=255)                                                                                                                                                                                                                                      aaaTTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCg           >  1:448762/1‑249 (MQ=255)                                                                                                                                                                                                                                       aaTTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGc        <  2:189622/251‑1 (MQ=255)                                                                                                                                                                                                                                         ttcttcCCACATAGCCTTCGTATCCTGCTGACATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTACGCTAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTCTCGCTATTCCTGTGACCGAGCtt      >  1:538537/1‑251 (MQ=255)                                                                                                                                                                                                                                              cccACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTtata  <  1:220783/250‑1 (MQ=255)                                                                                                                                                                                                                                                   aTAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGTGCCGCTTATCGCTAt                      <  2:519944/225‑1 (MQ=255)                                                                                                                                                                                                                                                   aTAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTAt                      >  1:519944/1‑225 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 ACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATA  >  NC_000913/4637896‑4638382

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