Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,532,063	C→T	intergenic (‑412/+374)	sgcX ← / ← yjhP	putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,532,063	0	C	T	100.0%	221.0 / NA	69	intergenic (‑412/+374)	sgcX/yjhP	putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/31);  total (38/31)

GGCTGGAGAGACCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCACTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTG  >  NC_000913/4531816‑4532306                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    ggCTGGAGAGACCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCAttct                                                                                                                                                                                                                                                  <  1:291751/251‑1 (MQ=255)  gCTGGAGAGACCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTTGTGCCCAGCATAATGTAAGAAAGCATCTTATCTCCGGAATAGATGCATAGGGAGTATAGATAATCAGCGCGAATCATtctc                                                                                                                                                                                                                                                 <  2:589034/251‑1 (MQ=255)    tGGAGAGACACGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGTTCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCg                                                                                                                                                                                                                                               >  2:598223/1‑251 (MQ=255)     ggAGAGACCAGGGAGTATTGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCGCTTTGTTGCAGGACGCTCCCTCTTCAGTTCGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGAGCTCCTGAATACATTCAGAGGGATTATAGATAGTCAGCTCAAATCATTCTcc                                                                                                                                                                                                                                                <  2:403471/249‑1 (MQ=255)        gagaCCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCttt                                                                                                                                                                                                                                           <  2:476687/251‑1 (MQ=255)            ccAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGa                                                                                                                                                                                                                                       <  1:542584/251‑1 (MQ=255)              aGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTGTGTTGCAGGACGCACCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCATCTCAAATCATGCTCCGCTTCATGAGc                                                                                                                                                                                                                                     <  2:418341/251‑1 (MQ=255)                 gATTATGGTGGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCa                                                                                                                                                                                                                                  <  1:254161/251‑1 (MQ=255)                          tGAAATAATGAGGGTGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGAAGGACGCTACCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTgg                                                                                                                                                                                                                            >  2:565822/1‑248 (MQ=255)                          tGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGAGCGCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTgg                                                                                                                                                                                                                            <  1:565822/248‑1 (MQ=255)                           gAAATAGTGAGAATGAATCCAACGCGGCAACTGTTCTAGCACCAATTAGATGCAGCACGCTCCAACTTCAATACGACCTGCACGTGGATACCGCAGCATCCCCATCAGCGGCTCTTCGCATGTAATGCCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACAACGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCaat                                                                                                                                                                                                                        >  2:241198/1‑249 (MQ=255)                              ataatGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGcc                                                                                                                                                                                                                     >  1:487635/1‑251 (MQ=255)                                aatGAGGATGAATAAAACGCGGGAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCGTCATGACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGTGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACAGCGGATCGCCTGAATAGATGCAGAGGGAGTATAGAGAATCAGCTCAAATCATTCTCCGCTGGATGGGCGCAGGCTGGCAGTGCCGt                                                                                                                                                                                                                   <  2:586568/251‑1 (MQ=255)                                     ggATGCATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTACCTCTTCATTACGTCATGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATACTGTAAGAAAACATCGGATCTCCTGAATAGATTCCGAGGGCGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTc                                                                                                                                                                                                                  >  2:258424/1‑247 (MQ=255)                                         gAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACg                                                                                                                                                                                                          >  1:505245/1‑251 (MQ=255)                                                 cgcgGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCTTGCACGTGGATCCCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGTCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATTGATAATCAGCTCAAATCCTTCTCCGCTTTATGAGCCCAGGCTTTCATTGCCGTCAGATAACGTTTTGCAc                                                                                                                                                                                                  <  2:355829/251‑1 (MQ=255)                                                        aCTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTAtg                                                                                                                                                                                           <  1:459126/251‑1 (MQ=255)                                                        aCTGTTATAGCACCACTTTGATGCAGGACGCTCCCTCTTCATTACGTCCTCCACGTGGAGACCACAGCATTTCCATCAGTGGATCTAAGAATGTACTGGCGCGAGGAAGGTCGTTGATAACCAGCATAATGTAAGAAAAAAACGGATCTCATGATTAGATTCAGAGGGAGTATAGATAATCAGCTCAAATGATTCTCCGCTTTATGAGCCCAGGTTGGGATTGCCGACAGCTAAAGTTTTTCACGCTTAtg                                                                                                                                                                                           >  1:393619/1‑251 (MQ=255)                                                                        tttGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGCAAACATCGGATCTCCTGCATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATACCGATTTGCACGCTTATGTGCTTCCTCCGGGCCg                                                                                                                                                                           >  1:323287/1‑251 (MQ=255)                                                                             tgCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACttt                                                                                                                                                                      >  1:346495/1‑251 (MQ=255)                                                                              gCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACttt                                                                                                                                                                      <  1:299090/250‑1 (MQ=255)                                                                                  gACGCTCCCTCTTCATTACGTCCTGCACGTGGATACACCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCAAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGcc                                                                                                                                                                 >  2:458248/1‑251 (MQ=255)                                                                                   aCGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCCTCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGcc                                                                                                                                                                 >  1:416155/1‑250 (MQ=255)                                                                                    cGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGcc                                                                                                                                                                            <  1:136703/238‑1 (MQ=255)                                                                                    cGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGcc                                                                                                                                                                            >  2:136703/1‑238 (MQ=255)                                                                                         ccTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATc                                                                                                                                                          <  1:510721/251‑1 (MQ=255)                                                                                               cATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCATGTCGTTGGTTCCCAGCATAATGTAATAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTTCCGTCAGATAACGTTTTTCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGaa                                                                                                                                                    <  1:624298/251‑1 (MQ=255)                                                                                                    cGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCGTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTgg                                                                                                                                               <  2:87372/251‑1 (MQ=255)                                                                                                           cACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAAcc                                                                                                                                        >  1:442259/1‑251 (MQ=255)                                                                                                            aCGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCt                                                                                                                                       <  1:471311/251‑1 (MQ=255)                                                                                                            aCGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCt                                                                                                                                       <  1:67433/251‑1 (MQ=255)                                                                                                               tGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCAAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTAAGAGGGAGTATAGATAATAAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGAt                                                                                                                                    >  1:557666/1‑251 (MQ=255)                                                                                                                           cATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTCAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAAACTGATCCCGCTTCTTcc                                                                                                                        >  1:275648/1‑251 (MQ=255)                                                                                                                           cATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGACTAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTcc                                                                                                                        >  1:328155/1‑251 (MQ=255)                                                                                                                           cATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTcc                                                                                                                        >  2:468705/1‑251 (MQ=255)                                                                                                                               tcCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAgg                                                                                                                    >  1:75448/1‑251 (MQ=255)                                                                                                                                ccATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAggg                                                                                                                   >  1:182924/1‑251 (MQ=255)                                                                                                                                ccATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAgg                                                                                                                    <  2:530163/250‑1 (MQ=255)                                                                                                                                    cAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGttt                                                                                                                <  1:5280/250‑1 (MQ=255)                                                                                                                                    cAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTGCCCAGGGttt                                                                                                                <  2:416155/250‑1 (MQ=255)                                                                                                                                     aGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGACGCTGGCTACACCTGATACCGCTTCTTCCCAGGGTTTg                                                                                                               >  2:612437/1‑250 (MQ=255)                                                                                                                                       cGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTgaga                                                                                                            <  1:467371/251‑1 (MQ=255)                                                                                                                                          aTCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATCATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTg                                                                                                         >  1:563644/1‑251 (MQ=255)                                                                                                                                              tAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGcc                                                                                                       <  2:473897/249‑1 (MQ=255)                                                                                                                                               aGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGtt                                                                                                    <  2:346495/251‑1 (MQ=255)                                                                                                                                                gAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTg                                                                                                   >  1:391043/1‑251 (MQ=255)                                                                                                                                                     gaagggCGCGTGGCAGGTCGTTGGTTCCCAGCATAAAGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCGCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTAGGTGCTTCCTCCGGGCCGACTTTTGGCCGGGAAGCAAGGAAGCTGGCTAGACCTGATACCGCTTCTTCCCAGGGTTTGAGAGTGCCGTTGCTGTc                                                                                              <  2:300919/247‑1 (MQ=255)                                                                                                                                                      aaTGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTc                                                                                              <  2:300928/250‑1 (MQ=255)                                                                                                                                                       aTGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGc                                                                                            >  1:170500/1‑251 (MQ=255)                                                                                                                                                            gcgTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCAt                                                                                       <  2:328251/251‑1 (MQ=255)                                                                                                                                                              gTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACCTCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATAACGCTTCTTTCCCGGGTTTGAGATTGCCGTTGCTGTAGCTGCAtt                                                                                      >  2:352323/1‑250 (MQ=255)                                                                                                                                                                 gCAGGTCGTTGGTTCCCAGAATAATGTAAGAACCCATCGGATCTCCTGAATAGATTCAGAGGGCGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGc                                                                                  >  1:496493/1‑251 (MQ=255)                                                                                                                                                                  cAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGc                                                                                  >  1:33161/1‑250 (MQ=255)                                                                                                                                                                    ggTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGcac                                                                                >  1:69841/1‑250 (MQ=255)                                                                                                                                                                    ggTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTACGCTTTATGAGCCCAGGCTGGCATTGCCGTCTGATAACGTTTTGCTCGCTTATGTGCTTCCTCCGGGACGACTTTTTGCCAGGAATAATGGAAGCTGGCTAAACCTGATACCGCTTCTTCCACGGGTTTGAGATTGCCGTTGCTGTCGATGCCTTACGcac                                                                                >  2:505467/1‑250 (MQ=255)                                                                                                                                                                               cccAGCATAATGTAAGAAAACATCGGATCTCCTGACTAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCAc                                                                    >  2:140287/1‑251 (MQ=255)                                                                                                                                                                                 cAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTg                                                                  >  1:234359/1‑251 (MQ=255)                                                                                                                                                                                  aGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTgg                                                                 >  1:530203/1‑251 (MQ=255)                                                                                                                                                                                             aGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGc                                                       >  1:391447/1‑250 (MQ=255)                                                                                                                                                                                                aaaCATCGGATCTCCTGACTAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACATTTTGAACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGGTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCt                                                   >  2:1962/1‑251 (MQ=255)                                                                                                                                                                                                      cGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGcc                                             >  1:394291/1‑251 (MQ=255)                                                                                                                                                                                                            cccTGAATAGATTCAGAGGGAGGATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTGTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATc                                       <  1:86599/250‑1 (MQ=255)                                                                                                                                                                                                             ccTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATc                                       <  2:451083/250‑1 (MQ=255)                                                                                                                                                                                                               tGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGc                                    >  1:346115/1‑251 (MQ=255)                                                                                                                                                                                                                 aaTAGATTCAGAGGGAGAATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGCCATCGGCAt                                  >  1:206820/1‑251 (MQ=255)                                                                                                                                                                                                                           gagGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGcc                        <  2:187787/251‑1 (MQ=255)                                                                                                                                                                                                                                aGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGATGGCATTGCCGTTAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGAc                    >  1:62086/1‑250 (MQ=255)                                                                                                                                                                                                                                      gATAATCAGCTCAAATCATGCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTg             <  1:90894/251‑1 (MQ=255)                                                                                                                                                                                                                                                 cAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCAtgttgt   <  1:598223/250‑1 (MQ=255)                                                                                                                                                                                                                                                  aaaTCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTtgtg  >  2:538910/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GGCTGGAGAGACCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCACTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTG  >  NC_000913/4531816‑4532306

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