Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,560,309	T→A	P124P (CCA→CCT)	yjiH ←	uncharacterized protein YjiH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,560,309	0	T	A	100.0%	189.1 / NA	60	P124P (CCA→CCT)	yjiH	uncharacterized protein YjiH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (28/32);  total (28/32)

GGCTCCGACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAAC  >  NC_000913/4560063‑4560548                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                ggCTCCGACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGa                                                                                                                                                                                                                                              >  1:20500/1‑250 (MQ=255)  gCTCCGACAAATTTAAACACCAGGATGACGGCGAAAGGTACGAGTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGAGGGCGCTACCGCTGGTGTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCGTCCTTCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAAt                                                                                                                                                                                                                                            <  1:409379/251‑1 (MQ=255)     ccGACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGcc                                                                                                                                                                                                                                         <  1:616387/251‑1 (MQ=255)         cAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCcaac                                                                                                                                                                                                                                     <  2:4747/251‑1 (MQ=255)                    acCAGGGTGACGGCTAAAGGCAAGATTACGGAAGAGCCCAGAACGGCAAACACCGAGACGCCGGCAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGGAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCCTACCGGCGGCGGAATCGGTGCTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGACTGCCCAACAACGGtttt                                                                                                                                                                                                                            >  2:383080/1‑249 (MQ=255)                      cAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGGAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTACCGCGTTCGGTCATTTTACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGAGCCCAGGAATGCCCAACAACGGTTTTAGca                                                                                                                                                                                                                        <  2:183621/251‑1 (MQ=255)                                      ggTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAActg                                                                                                                                                                                                        <  2:490776/251‑1 (MQ=255)                                          cGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCCACAACGGTTTTAGCCCAGGTGTCATTAACtgctgc                                                                                                                                                                                                    >  1:276390/1‑251 (MQ=255)                                             ttACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGACGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcg                                                                                                                                                                                                 >  1:562749/1‑251 (MQ=255)                                                  gACGTACCCAGAACGGCAAAAACAGCGCCGCCGGAAGCGAAGTCGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcggcgcg                                                                                                                                                                                            >  2:319690/1‑251 (MQ=255)                                                   aCGTGCCCAGAAAGGAAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCGTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTGTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGATTGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGt                                                                                                                                                                                           <  1:402316/251‑1 (MQ=255)                                                     gTGCCCAGAAAGTCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGGTTGATATGCGGCGAAAATCTCTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTaa                                                                                                                                                                                         <  2:331410/251‑1 (MQ=255)                                                           aGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGGGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAAGTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCAGTAACTGCTGCGCGGCGCGTAATccacc                                                                                                                                                                                   <  2:562749/251‑1 (MQ=255)                                                                ggCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGAc                                                                                                                                                                              <  2:20500/251‑1 (MQ=255)                                                                    aaCACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATc                                                                                                                                                                          <  2:440170/251‑1 (MQ=255)                                                                      cacCGCGACGCCGGAAGAGACGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCgg                                                                                                                                                                        >  2:361437/1‑251 (MQ=255)                                                                       acCGCGGCGCCGGGAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGCGTTTTGCAGGATAGCGCTCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGATTTAGCACAGGTATCATTAACTGCTGCTCGGCGCGTAATTCACCCAGACCATCGGt                                                                                                                                                                       >  2:81136/1‑251 (MQ=255)                                                                                          agtagtTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATcccc                                                                                                                                                       >  1:50818/1‑248 (MQ=255)                                                                                             tgtTGGTGAGTATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGGCATTAACTGCTGCGCTGCGCGTAATCCACCAAGGCCATCGGTAATTGAAATAATCCCCagaga                                                                                                                                                  <  2:38701/249‑1 (MQ=255)                                                                                                       gatgGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGGCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTGAGCACAGGTGTGATGAACTGCTGCGCGGCGCGGAATCCGCCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAAtgact                                                                                                                                        <  1:622122/250‑1 (MQ=255)                                                                                                          ggCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCAACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACg                                                                                                                                    >  1:287125/1‑251 (MQ=255)                                                                                                          ggCGCGACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAAtgactgac                                                                                                                                     <  1:383080/250‑1 (MQ=255)                                                                                                             gcTACCGCTGGTTTGCTATGCGGCGAAAATCACTTTATCGCGTTCGGTAATCTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCCGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTg                                                                                                                                 >  1:412544/1‑251 (MQ=255)                                                                                                              cTACCGTTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGAGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGCGAATGACTGACGGTGc                                                                                                                                >  1:124407/1‑251 (MQ=255)                                                                                                                   gCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCt                                                                                                                           <  2:235542/251‑1 (MQ=255)                                                                                                                           gATATGCGGCGAAAATCACTTTATCGCGTTCGGCAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGCGTTTTGCCGGTTAGCGATCAGCGCCCGCGGGCAGATCCCAGGAATTCCCAACAACGGTTTTAGCACAGGTGTCCTTAACTGCTGCGCGCCGCGTAATCCACCAAGACCATCGGGAACTGAAATACTCCCCAGAGAGAGAATGACTGGCGGTGCCAGCTCCAgcgc                                                                                                                    >  2:4106/1‑250 (MQ=255)                                                                                                                             tatGCGGCGAACATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGCAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCaa                                                                                                                  >  2:55802/1‑250 (MQ=255)                                                                                                                                   gcgAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGCCATCGGTGTTTTGCAGGTTAGCGATCAGCCCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGAGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATACCCAGAGAGAGCATGACTGACGGGGCCAGCTCCAGCGCAAAGAGaa                                                                                                            >  1:607008/1‑250 (MQ=255)                                                                                                                                     gAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGGGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGc                                                                                                         >  2:332587/1‑251 (MQ=255)                                                                                                                                     gAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGc                                                                                                         >  2:455852/1‑251 (MQ=255)                                                                                                                                     gAAAATCACATTATCTCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGGTAGCGATCAGCGCCAGCGAGCAGATCCCAGGCATGCCCAACAACGGGTTTAGCACAGGTGTCATTACCTGCTGCGCGGCGCGTAATCCACCAAGACCATCGTGAATTGAAATCATCCCCAGAGGGAGAATGCCTGACGGTGCCCGCTCCAGCGCACAGAGCAAGc                                                                                                         >  2:115957/1‑251 (MQ=255)                                                                                                                                      aaaaTCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGcc                                                                                                        >  2:119615/1‑251 (MQ=255)                                                                                                                                                   cgcgTTCGGTCATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCTGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATTCCCAACATCGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTTACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGcccc                                                                                            <  1:515350/250‑1 (MQ=255)                                                                                                                                                     cgTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCa                                                                                                       >  2:93631/1‑237 (MQ=255)                                                                                                                                                     cgTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCa                                                                                                       <  1:93631/237‑1 (MQ=255)                                                                                                                                                       ttCGGTAATTTAAACTTCCTGCGCCCGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCgcgc                                                                                       >  2:40275/1‑251 (MQ=255)                                                                                                                                                           gTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCAcc                                                                                   <  1:259548/251‑1 (MQ=255)                                                                                                                                                             aaTTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTGAGCCCCCGCTCCAccc                                                                                  <  2:287125/250‑1 (MQ=255)                                                                                                                                                                ttCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGAGGCGGCCTCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAACGAGAAAGCCATCTTTAGCCCCCGCGCCACccgcgc                                                                              >  2:393606/1‑251 (MQ=255)                                                                                                                                                                  cACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACccgcgccg                                                                            <  1:364580/251‑1 (MQ=255)                                                                                                                                                                                  ctctTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGGGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGtggtgg                                                            >  2:101413/1‑251 (MQ=255)                                                                                                                                                                                    cgttGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGAGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAAGTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGtggtgg                                                            <  2:276390/247‑1 (MQ=255)                                                                                                                                                                                        ggTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCCCAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAAGCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCAAGAGCAAAGAGAAAGCCATCGTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGt                                                      <  2:26577/251‑1 (MQ=255)                                                                                                                                                                                           cATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCAccc                                                                                  >  2:488432/1‑220 (MQ=255)                                                                                                                                                                                           cATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCAccc                                                                                  <  1:488432/220‑1 (MQ=255)                                                                                                                                                                                                    ggcggcATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGCGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGc                                           >  1:416370/1‑250 (MQ=255)                                                                                                                                                                                                       ggcATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGtt                                       <  1:332587/251‑1 (MQ=255)                                                                                                                                                                                                          aTCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGcc                                    <  1:361437/251‑1 (MQ=255)                                                                                                                                                                                                                gtTTCGCACGTTAGCGATCAGCGCCTGCGAGCAGAGCCCAGTAATGCCCAACAACGGTTTGAGCACAGGTGTCATTAACTGCTGCTCGGCGCGTAATCCACCAAGACCATCGGTCATTGAAAGAATCCCCAGAGAGAGGATTACTGACGCTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTAAGCCCCCGCGCCTCCCGCTCCGCGAAATGAGGGGGTGTCTCCGTTCTCACCCGGCAGTTGCCCAAAc                               <  2:610581/250‑1 (MQ=255)                                                                                                                                                                                                                    tGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCATGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAAGTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGTCAGTTTCCCAAACGATc                           <  2:316470/250‑1 (MQ=255)                                                                                                                                                                                                                     gCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCg                         >  1:215410/1‑251 (MQ=255)                                                                                                                                                                                                                             gCGATCAGCGCCATGGAGCAGATCCAAGGAATGCCCAACAACGGTTGTAGCACAGGTGTCATTAACGGATGCGCGGCGCGTAATCCACCAAGACCATCGGAAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCGCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCGCCGTTCGCACCAGGAAGTTGCCCAAACGATCTGTTCAGAAc                 <  2:22011/251‑1 (MQ=255)                                                                                                                                                                                                                              cGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCGCCAGCGCAAAGAGAAAGCCATCTGTAGCCCCCGCGCCGCCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTACAGAACg                <  2:311265/251‑1 (MQ=255)                                                                                                                                                                                                                              cGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAAGCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGTTGGCTCCGTTCGCACCAGGCAGTTGCCCATACGATCCGTTCAGAACg                <  2:573694/251‑1 (MQ=255)                                                                                                                                                                                                                                  cAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGaa             <  2:50818/250‑1 (MQ=255)                                                                                                                                                                                                                                       ccAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGAGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGAGGCTCCGTTCGCACCAGGCAGCTGCCCAAACGATCCGTTCAGCACGGaaaattc        >  1:548053/1‑247 (MQ=255)                                                                                                                                                                                                                                        cAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTATCTGCTGCGCGGCGCGTAATCCACCAGGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGGCTGACGGTGCCAGCTCCAGGTCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGGGGTGGCGCCGGTCGTACCAGGCAGTTGACCAAACGATCCGTTCAGAACGGAAAAATCaa      <  1:31049/251‑1 (MQ=255)                                                                                                                                                                                                                                          gCGAGCAGATCCCAGGAAGTCCCATCAACGGTTTTGGCTCAGTTGTCAATAACTGCTGCGCCGCGAGTACTACACAAAGACCATCGGTAATTGAAATAATCCCCAGAGATATAATGACTGCCGTTGCCAGCCCCAGCGCAAAGAGAAAGCCTGCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGGGGCTACGTTCGCACCAGGCAGTTGCCAAACCAATCCGTCCAGATCAGAAATATCaaaa    <  2:273668/251‑1 (MQ=255)                                                                                                                                                                                                                                           cGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCaaaa    >  2:392254/1‑250 (MQ=255)                                                                                                                                                                                                                                            gAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAAc  >  2:178317/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                GGCTCCGACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAAC  >  NC_000913/4560063‑4560548

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