Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,585,196	G→A	R431C (CGC→TGC)	ydeP ←	putative oxidoreductase YdeP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,585,196	0	G	A	100.0%	131.0 / NA	40	R431C (CGC→TGC)	ydeP	putative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/17);  total (23/17)

CGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATCCCCATTCCGTAACAGATAATGGTGCGTTCGGCAGCGGCATATGCGTCAGCCAGTTCG  >  NC_000913/1584952‑1585440                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     cGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCAcc                                                                                                                                                                                                                                                <  1:217660/251‑1 (MQ=255)       aGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACg                                                                                                                                                                                                                                          >  2:379356/1‑251 (MQ=255)        ggACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCACCCATGCAGATCAATGCTCGAGCCTGCCCCGTACCTATTGCTTGCATGCTGGCACTTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACg                                                                                                                                                                                                                                          >  1:248576/1‑250 (MQ=255)            cGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGcc                                                                                                                                                                                                                                      >  1:62272/1‑250 (MQ=255)              gCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGcc                                                                                                                                                                                                                                      >  1:54381/1‑248 (MQ=255)                    tataGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAActct                                                                                                                                                                                                                             >  1:327109/1‑251 (MQ=255)                            ggccggggggtccaagaaTGAGAGCTGTTAAGCTTAGTCGCTACGTGTACCGCCAAATCTTAGTGCGTTAACGGTACAGCGCTTGCTTCCCGATCTGGCATTGCCAGCTTAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTTCTTGCATGCTGGCAATTGCAGCATGTCCAGGGGCTTGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGTCAGAAACTCTGCAGCCgg                                                                                                                                                                                                                     <  2:294401/234‑1 (MQ=255)                                gTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtc                                                                                                                                                                                                                 >  1:297311/1‑251 (MQ=255)                                   cGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCCATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCg                                                                                                                                                                                                              >  2:80029/1‑251 (MQ=255)                                          aaGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATAc                                                                                                                                                                                                       <  2:105846/250‑1 (MQ=255)                                            gATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCg                                                                                                                                                                                                     >  1:94282/1‑251 (MQ=255)                                             aTGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCg                                                                                                                                                                                                     >  2:85993/1‑250 (MQ=255)                                             aTGAGAGCGGATAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCACCATGCAGCTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGATTCTCGGTGATACCgc                                                                                                                                                                                                    >  2:231396/1‑250 (MQ=255)                                                                          gtACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgtgt                                                                                                                                                                       <  2:297311/251‑1 (MQ=255)                                                                                              gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCACTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATAcc                                                                                                                                                   <  2:124529/251‑1 (MQ=255)                                                                                                 aaCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGc                                                                                                                                                >  1:317442/1‑251 (MQ=255)                                                                                                     gTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCa                                                                                                                                            <  1:34503/251‑1 (MQ=255)                                                                                                                   tCCCGATCTGGCATTGCCAGCGCAAAATTCCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCCCCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATg                                                                                                                                <  2:97812/249‑1 (MQ=255)                                                                                                                               aTTGCCAGCGAAAAATTGCCCCCCAGGCAGAGCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGGAGGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACACTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCaa                                                                                                                  <  2:85218/251‑1 (MQ=255)                                                                                                                                       cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCACGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAAcagc                                                                                                              >  1:212594/1‑247 (MQ=255)                                                                                                                                           aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGAc                                                                                                      <  1:296854/251‑1 (MQ=255)                                                                                                                                                 ccccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCACTTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTg                                                                                                >  2:376018/1‑251 (MQ=255)                                                                                                                                                         gCAGATCATTGCTCGATCCTGCCGCGTACATATTGCTGTCATGCTGGCAAGTGAAGCGTGGCCAGGCGCCAGGGGGGGGGTGAGGCCGTAGCACCCACCCAGGCGAGCCATAAACTCTGCAGCCGTTTTCTCGGTGATACCGGCGGTTCGGTCGCCCTGTACATTAGAGTGTCCACTTAGTGGACAGATACCCGCACCAGGCTTGTCATTTTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACg                                                                                         <  2:15310/250‑1 (MQ=255)                                                                                                                                                           aGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGttt                                                                                      >  1:162026/1‑250 (MQ=255)                                                                                                                                                           aGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTc                                                                                      >  1:97539/1‑251 (MQ=255)                                                                                                                                                                 aTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGt                                                                                 >  1:19889/1‑250 (MQ=255)                                                                                                                                                                    ctagAGCCTGCCCCGTACATTTTGCTTGCATGATGGCATTTGCCGCATGTCCCAGAGAATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGCTACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATAACCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCCGATTGACCAGTTGCTGTACGTGCTGGGTAccc                                                                             >  2:75252/4‑250 (MQ=255)                                                                                                                                                                       gAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCCCCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGt                                                                          <  2:188118/251‑1 (MQ=255)                                                                                                                                                                       gAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGGTACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGt                                                                          >  2:70301/1‑251 (MQ=255)                                                                                                                                                                       gAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGGGAAGCCATAGCACTCAACCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGGCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCGGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGAGTACCATGt                                                                          >  1:353414/1‑251 (MQ=255)                                                                                                                                                                               cccGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGAGAAGCCATAGCACTCAACCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGGTCGGTCGCCATGTACATTAGAGTATCAACGTAGTGGGCAGATACCCACACCAGGCTTGCCCATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGttattgcgg                                                                  >  2:285651/1‑244 (MQ=255)                                                                                                                                                                                      atatTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATc                                                           >  1:175044/1‑251 (MQ=255)                                                                                                                                                                                               gcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATCCCCATTCCg                                                  <  2:94282/251‑1 (MQ=255)                                                                                                                                                                                                               agcaTGTCCAGGTGCATGAGGTTGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATCCCCATTCCGTAACAGATAATGGTGc                                  <  2:327109/251‑1 (MQ=255)                                                                                                                                                                                                                   tGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTAACCAAGACGAGACAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGAAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGGGGGATCCCCATTCCGTAACAGATAATGGGGCgatc                              >  1:108643/1‑248 (MQ=255)                                                                                                                                                                                                                     tCCAGGTGCATGAGGTGGGGTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGTCGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGAGTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATCCCCATTCCGTAACAGATAATGGTGCGTTCg                             <  2:21475/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggTGAAGCCATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATCCCCATTCCGTAACAGATAATGGTGCGTTCGGCAGCGGCATATGCGTc            >  1:264739/1‑249 (MQ=255)                                                                                                                                                                                                                                           aacncATAGCACTCATCTAGACGAGCCAGAAACGCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTAAATTAGAGTGTCCACGTAGTGGGCAGATACCCGCCCCAGGCTCGCCAATGTTACCTTTCATCACCAGCAGACTGACCAGTTGCTGTCCGTGCTGGGTACCATGTTCGTGCTGAGTGATCACCACTCCGTAACAGATAATGGTGCGTGCGTCAGCGGCATATGCGccagccag      <  2:353414/247‑1 (MQ=255)                                                                                                                                                                                                                                              ccATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATCCCCATTCCGTAACAGATAATGGTGCGTTCGGCAGCGGCATATGCGTCAGCCAGTTc   <  1:351505/251‑1 (MQ=255)                                                                                                                                                                                                                                               cATAGCACTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATCCCCATTCCGTAACAGATAATGGTGCGTTCGGCAGCGGCATATGCGTCAGCCAGTTCg  <  1:379356/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGGCTTGCCAATGTTACCTTTCATCAACAGCAGATTGACCAGTTGCTGTACGTTCTGGGTACCATGTTCGTGCTGAGTGATCCCCATTCCGTAACAGATAATGGTGCGTTCGGCAGCGGCATATGCGTCAGCCAGTTCG  >  NC_000913/1584952‑1585440

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: ‑