Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,285,238	G→A	intergenic (‑25/+175)	yjcH ← / ← acs	DUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,285,238	0	G	A	100.0%	243.3 / NA	75	intergenic (‑25/+175)	yjcH/acs	DUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/36);  total (39/36)

ATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCA  >  NC_000913/4284993‑4285476                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                aTCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGATAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGAtt                                                                                                                                                                                                                                            >  1:446232/1‑251 (MQ=255)    acaATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGtt                                                                                                                                                                                                                                         <  1:186907/251‑1 (MQ=255)         cACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTaa                                                                                                                                                                                                                                     <  1:16118/250‑1 (MQ=255)         cACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAAt                                                                                                                                                                                                                                    <  1:108714/251‑1 (MQ=255)          aCTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGGTCAGTAAAATAAAGCCGCTATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGGCATGTAAGGATTGTTTTAAtt                                                                                                                                                                                                                                   >  1:186190/1‑251 (MQ=255)             ccAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTccc                                                                                                                                                                                                                                <  1:232070/251‑1 (MQ=255)              cAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGGTCAGTAAAATAAAGCCGGTATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGCCTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCg                                                                                                                                                                                                                               >  2:479766/1‑251 (MQ=255)                      atgggAATGCCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCTGGCCCAGCCAGCCGGGCGCGAAGGCTATAAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTATCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCAGTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCACTTa                                                                                                                                                                                                                       <  1:433595/248‑1 (MQ=255)                              aCCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGgaga                                                                                                                                                                                                                <  1:430354/250‑1 (MQ=255)                                cgcgTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTGGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAgg                                                                                                                                                                                                             <  1:127305/251‑1 (MQ=255)                                          gCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCCGTAAAATAAAGCCGATATACACTGCCAGCATAATCATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTTAAATGCGCATTGGCTTCTATCCCCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGCGAAgggtcaagcacg                                                                                                                                                                                                   >  2:449489/1‑241 (MQ=255)                                             ggTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGt                                                                                                                                                                                                >  1:145722/1‑251 (MQ=255)                                                      aTTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGGGAAGGTTAACGCTCGGGTAACCCTTGc                                                                                                                                                                                       >  2:14459/1‑251 (MQ=255)                                                         aaGCGTCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGGAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCg                                                                                                                                                                                     <  2:430440/249‑1 (MQ=255)                                                                     cagccagcCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCgg                                                                                                                                                                         >  1:392747/1‑250 (MQ=255)                                                                       gccagccGGTCGAGAAGACGCTCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGAGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGACATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGGGCAGGTTAACGCTCGGGTAACACTTGCCGAATGTAGGCCGGAt                                                                                                                                                                       >  2:228882/1‑250 (MQ=255)                                                                        ccagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGACGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATa                                                                                                                                                                      >  2:161808/1‑250 (MQ=255)                                                                           gccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGc                                                                                                                                                                  <  1:349348/251‑1 (MQ=255)                                                                              gggCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCCTTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGCGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTAGGGTAACCCTTGCCGCATGTAGGCCGGATACGGCGtt                                                                                                                                                               >  2:352273/1‑251 (MQ=255)                                                                                           aTCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGATAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                                  <  2:446232/251‑1 (MQ=255)                                                                                            tCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACCCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCAAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCATTATGGGAGAAGGTTAACGCTCGGGAAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCCCATCCgg                                                                                                                                                 >  1:284523/1‑251 (MQ=255)                                                                                                aaaaaTAAAGCCGATATAAACTGCCAGCATAAAAATCAACAGGATGGTGGCAAACCGTTGCCGTTTTTCGCCGAACTCCCTGAAATGCGCATTGTCTTCTAGACGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATTGGAGAAGGTGAACGCTCGGGTAACCCTTTCCGAATGTAGGCCGGATAAGGCGTTTACGCTGCATCCGGcaat                                                                                                                                             <  2:41389/250‑1 (MQ=255)                                                                                                  aaaTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatc                                                                                                                                            >  2:471721/1‑250 (MQ=255)                                                                                                   aaTAAAGCCGCTATAAACTGCCAGCATAATAATCGCCAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACACCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGccctca                                                                                                                                           >  2:61206/1‑245 (MQ=255)                                                                                                     aaaaGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATg                                                                                                                                        <  1:292701/250‑1 (MQ=255)                                                                                                          ccGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGa                                                                                                                                   <  1:177202/251‑1 (MQ=255)                                                                                                                taAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAGATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGGAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑Ac                                                                                                                            >  1:399879/1‑251 (MQ=255)                                                                                                                      gCCAGCATAATAAGCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGGTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTc                                                                                                                      <  2:249070/251‑1 (MQ=255)                                                                                                                           cagaacaatCGACAGGATGGGGGCAAACCGTGGCCGTTTTTCGACTAACGCCCTGAAATGCGCAGTGTCTTCTATCCGCTGATAAATAGTGCAAGTCATCACAGATTCTCCCGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTATCCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCTCAGCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTc                                                                                                                 <  2:221669/245‑1 (MQ=255)                                                                                                                              aataatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTAGCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGGAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCtt                                                                                                               <  2:186190/250‑1 (MQ=255)                                                                                                                                taatCGACAGGATGGTGGCACACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGCTAAGGCGTTTACGCCGCCTCCGGCAATCAATGCCTGATTCG‑ACGCTGTCGCGTCTTATc                                                                                                            >  2:10821/1‑251 (MQ=255)                                                                                                                                taatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATc                                                                                                            >  1:384104/1‑251 (MQ=255)                                                                                                                                  atCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGcc                                                                                                                                                                           >  1:91502/1‑187 (MQ=255)                                                                                                                                  atCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGcc                                                                                                                                                                           <  2:91502/187‑1 (MQ=255)                                                                                                                                    cGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGc                                                                                                        >  2:345310/1‑251 (MQ=255)                                                                                                                                     gACAGTATGGTGTCAAACCGTTGCAGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGTGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGcc                                                                                                       <  2:399879/251‑1 (MQ=255)                                                                                                                                     gACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGTTAAGGCGTTTACGCCGCATCCGGCCATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGcc                                                                                                       >  2:290973/1‑251 (MQ=255)                                                                                                                                     gACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGACATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGcc                                                                                                       >  1:305074/1‑251 (MQ=255)                                                                                                                                      aCAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCt                                                                                                      >  1:376669/1‑251 (MQ=255)                                                                                                                                        aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCt                                                                                                                <  1:77025/239‑1 (MQ=255)                                                                                                                                        aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCt                                                                                                                >  2:77025/1‑239 (MQ=255)                                                                                                                                          gATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACaa                                                                                                  <  1:14459/251‑1 (MQ=255)                                                                                                                                             ggtggCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAAcc                                                                                               <  1:345310/251‑1 (MQ=255)                                                                                                                                                 gCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATGGTCTTCTAGCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTa                                                                                           <  2:112722/251‑1 (MQ=255)                                                                                                                                                  cAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCGTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGTCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                          <  1:212053/251‑1 (MQ=255)                                                                                                                                                    aaCCGTTGCCGTTTTTCTACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTGTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCg                                                                                        <  2:172300/251‑1 (MQ=255)                                                                                                                                                       cGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCTTTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCGCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTGAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACt                                                                                     <  1:449489/251‑1 (MQ=255)                                                                                                                                                        gTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGCTAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAAAAGGTTAACGCTCGGGTAACCCTTGCCGCATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGACTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTCCCACCCGTTACCGACTc                                                                                    >  2:170454/1‑251 (MQ=39)                                                                                                                                                          tGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCg                                                                                   <  2:55764/250‑1 (MQ=255)                                                                                                                                                             cGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGGCTTATCAGGCCTACAAACCGTTACCGACTCGCATc                                                                               <  1:136418/251‑1 (MQ=255)                                                                                                                                                               tttttCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGTTAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCGCATCgg                                                                             >  2:26399/1‑251 (MQ=255)                                                                                                                                                                 tttCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCggg                                                                            <  1:45787/250‑1 (MQ=255)                                                                                                                                                                                  aactGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATg                                                          >  2:32357/4‑251 (MQ=255)                                                                                                                                                                                  aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                         >  2:125284/1‑157 (MQ=255)                                                                                                                                                                                  aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                         <  1:125284/157‑1 (MQ=255)                                                                                                                                                                                  aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATg                                                          <  2:376669/251‑1 (MQ=255)                                                                                                                                                                                   aaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATgg                                                         >  2:436014/1‑251 (MQ=255)                                                                                                                                                                                    aTGCGCCTTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAACCCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGc                                                        >  2:179510/1‑251 (MQ=255)                                                                                                                                                                                           ttGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATcgcg                                                  >  2:25629/1‑250 (MQ=255)                                                                                                                                                                                                 tctATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCCCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTGACCGACGCGCATCGGGCAATTGTGGGTTACGATGGCATCACGATAGcc                                            <  2:5813/250‑1 (MQ=255)                                                                                                                                                                                                      ccGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTtc                                      <  1:266951/251‑1 (MQ=255)                                                                                                                                                                                                         cTGATAAATAGTGCCATTCATCACAGCTTCTCCAGAGGAAATGAAAGGGTTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTCCAACCCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGTCATCGCGATAGCCTGCTTctct                                   >  1:50020/1‑251 (MQ=38)                                                                                                                                                                                                            aTAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTACCGCTCCCTTATGGGCGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTACGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCGCACCGGGCAATTGTGGGTTACGATGGCATCGCTATAACCTtcttctcttc                                 >  2:315467/1‑250 (MQ=38)                                                                                                                                                                                                             tAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCa                                >  1:262714/1‑250 (MQ=255)                                                                                                                                                                                                               aaTAGTGCCATTCATCCAAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTCTGGGCGAAGGTTAACGCTCGGGTAACCCTTGCCGACTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGCTGCG‑CCGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                          >  2:220931/1‑190 (MQ=255)                                                                                                                                                                                                               aaTAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                          <  1:220931/190‑1 (MQ=255)                                                                                                                                                                                                                aTAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGCCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagc                             >  1:177646/1‑250 (MQ=255)                                                                                                                                                                                                                      ccATTCATCACAGATTCTCCAGGGGGAAGGGAAGGAGTGTTTTAATTCCCGCTCCCTTAGGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACGCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCATGCAGCTtct                      <  2:201840/251‑1 (MQ=38)                                                                                                                                                                                                                        aTTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATc                                                     >  2:238413/1‑218 (MQ=255)                                                                                                                                                                                                                        aTTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATc                                                     <  1:238413/218‑1 (MQ=255)                                                                                                                                                                                                                             tcaCAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑CCGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGAAAGCCTGATTTTCTTCAACCAGCTTCTCGACTAc               >  2:362830/1‑251 (MQ=255)                                                                                                                                                                                                                             tcaCAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTAc               >  1:436101/1‑251 (MQ=255)                                                                                                                                                                                                                                        ccAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATc      <  1:25629/249‑1 (MQ=255)                                                                                                                                                                                                                                           gagGTAATATAATGATTGGTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAAGGCTCGGTTATCCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGTAATCAATGCCTGATGCGCAC‑CTGTCACGTCTTCTCAAGACTACCAACCGTACCCAACTCGCATCGGTCAATTGTGGTTTACGATGGCATCGCGATAACATGCTACTCTTCACGCAACTTTTCAACTACACCATAATCGGCa  >  2:367975/1‑250 (MQ=255)                                                                                                                                                                                                                                            agGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCCGGGTCGGCa  >  2:189444/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                ATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCG‑ACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCA  >  NC_000913/4284993‑4285476

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