Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,985,726	A→G	intergenic (‑134/‑121)	yqeF ← / → yqeG	putative acyltransferase/putative transporter YqeG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,985,726	0	A	G	100.0%	250.5 / NA	74	intergenic (‑134/‑121)	yqeF/yqeG	putative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/36);  total (38/36)

ATTAACGCTTTCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTT  >  NC_000913/2985480‑2985968                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aTTAACGCTTTCACGACCAGACTGCCAAGTTCCACGGCGGAATGACCCGCTAACGCCCCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCGCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGGGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGtt                                                                                                                                                                                                                                                <  2:116474/251‑1 (MQ=255)               gACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAAcc                                                                                                                                                                                                                                  >  2:27135/1‑251 (MQ=255)                   aGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACccacc                                                                                                                                                                                                                              <  2:278608/251‑1 (MQ=255)                                 cACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCCCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGa                                                                                                                                                                                                                <  1:304767/251‑1 (MQ=255)                                   cggcggAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACAGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGAt                                                                                                                                                                                                               >  2:5891/1‑250 (MQ=255)                                          aTGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGtt                                                                                                                                                                                                       <  1:421045/251‑1 (MQ=255)                                               ccGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAAtct                                                                                                                                                                                                  <  1:438665/251‑1 (MQ=255)                                                 gCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTcc                                                                                                                                                                                                <  1:469609/251‑1 (MQ=255)                                                     aCGCACCACGAACGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCAt                                                                                                                                                                                             >  2:461492/1‑250 (MQ=255)                                                        caccacGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCTCCAATTAACAGTGCAATAGTAAATGATTGCTAATTACAGTTATTTTAATTGTCTAAAAAAAGTGATTTCTATCACAAAGGCAATATGCATGAGCAGCAGTCAGAGTAATAACTCGAACTTAGGGTTTGCAACCACCAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTACCATACtt                                                                                                                                                                                         <  2:286812/251‑1 (MQ=255)                                                               aaaGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATGGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACTTAGGGTTTGCAACCTCTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTgg                                                                                                                                                                                  <  1:220748/251‑1 (MQ=255)                                                                            ggTGTCCGTAACGCCCCGACAATCACAACGTNTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAAtt                                                                                                                                                                     <  2:246790/251‑1 (MQ=255)                                                                                     aaCGCCCCGACAAGCACAACGTCTTTCATAACCACCTCGCACCAATTAACAGTGCAATAGGAAATGATTGTGAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTGTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCAGGTCGGGTTTGCAACCTCTATCACTCCAATAGAGGGGTAGATAGGGCTTTAATCTACCATACGTATCTTGGTTTGTGGTAAATGGCCCTCCAt                                                                                                                                                            <  2:405376/251‑1 (MQ=255)                                                                                            cGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGGGCAATAGTAAATGATTGTTAATTACAGTTAGCTTAAGTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGATCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATATAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAAGTGCCCTCCATTCTGttt                                                                                                                                                     <  2:78184/251‑1 (MQ=255)                                                                                             gACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTa                                                                                                                                                    <  1:84993/251‑1 (MQ=255)                                                                                                  tCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTg                                                                                                                                               >  2:95758/1‑251 (MQ=255)                                                                                                         gTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAgatga                                                                                                                                        <  1:27135/251‑1 (MQ=255)                                                                                                             ttCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAgatgat                                                                                                                                       <  1:36787/248‑1 (MQ=255)                                                                                                             ttCATCAAAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAACAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATATAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACg                                                                                                                                    >  2:130925/1‑251 (MQ=255)                                                                                                              tcatcaCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACg                                                                                                                                    <  1:223629/250‑1 (MQ=255)                                                                                                               catcaCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGtt                                                                                                                                  <  2:57985/251‑1 (MQ=255)                                                                                                                ataacaACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATATCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTc                                                                                                                                 >  2:371729/4‑251 (MQ=255)                                                                                                                  cacaACCTCGCACCAATTAACAGTGAAATAGTAAATGATTGTTAAATACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTAGCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCGCCAATAGAGGGGTAGATATGGCGTTAATCTTCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGAGACGTTCt                                                                                                                                <  2:126942‑M1/250‑2 (MQ=255)                                                                                                                   acaACCTCGAACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTCATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGCCACGTAGGGTTTGCAACCAATAACCCACCAATAGAGGGGAAGATAGGGCGTTAATCTCCCATACTTAACCTAGTTTATGGTAAATTGCCCTCAATTCTGTTTAATTTATAGATGATACGTTCta                                                                                                                               >  1:477069‑M1/1‑248 (MQ=255)                                                                                                                     aaCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaa                                                                                                                             >  2:99457‑M1/1‑246 (MQ=255)                                                                                                                            cACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGAAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTCATTTGTAGATGATACGTTCtaaaaaaacccc                                                                                                                     >  2:12344‑M1/1‑239 (MQ=255)                                                                                                                               cAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATCTGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGTTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCAGTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                  <  1:156443‑M1/251‑16 (MQ=255)                                                                                                                               cAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATATAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTTGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                  <  1:96294‑M1/251‑16 (MQ=255)                                                                                                                                    aaCAGTGCGATAGTAAATGATTGGTCATTCCAGTTATCTTAATTGTTTAAAAAAAGGGATTTTTATCACAAAGCAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATCGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTGTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttt                                                                                                             >  2:34796‑M1/1‑231 (MQ=255)                                                                                                                                     aCAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGAGATTTTTATCACAAAAGAAATATGCCTGAGCAGCAGTTAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAACTTGCCCTACATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttac                                                                                                            >  2:34463‑M1/1‑230 (MQ=255)                                                                                                                                       aGTGCAATAGTAAATGATTGTTAATTACAGTAATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAACGAAATATTCCTGAGCAGCAGTCAGAGACATAACTTGCACGTAGGATTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTACCATACTTAACCTGGGTTATGGTAAATTGCCATCCCTTCTGTTTAATTTGTAGATGATACGTTAtaaaaaaaccccgacatttgccg                                                                                                          >  2:166427‑M1/1‑228 (MQ=255)                                                                                                                                        gTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATGCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgg                                                                                                         <  2:364487‑M1/251‑25 (MQ=255)                                                                                                                                                    aTGATTGTTAATTACAGTTCTCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcat                                                                                             >  1:364389‑M1/1‑215 (MQ=255)                                                                                                                                                    aTGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcat                                                                                             >  1:486294‑M1/1‑215 (MQ=255)                                                                                                                                                      gATTGTTAATTACAGTTATCTTAAGTGGTTAAAAAAAATGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaatccccgacatttgccggggttttgagcatat                                                                                           <  2:132715‑M1/251‑39 (MQ=255)                                                                                                                                                        ttGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacg                                                                                         >  2:40564‑M1/1‑211 (MQ=255)                                                                                                                                                          gTTAATTACATTTATCTTAATTGATTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGAAGGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGAGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaacccgacatttgccggggttgtgagcataacgt                                                                                        >  2:469474‑M1/1‑209 (MQ=255)                                                                                                                                                          gTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgta                                                                                       <  2:289982‑M1/251‑43 (MQ=255)                                                                                                                                                             aaTTACAGTTATCTGAATTGGTTGAAAAAATTGATGTTTATCACAAAGGAAAGATTCCTGAGCAGCAGTCAGTGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGGTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCAGTCTGTTTAATTTGTTGATGATACTGTCgaaaaaaaccccgacatttgccggggttgtgatcataacgtaaa                                                                                     <  2:27701‑M1/250‑45 (MQ=255)                                                                                                                                                              aTTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgc                                                                                   <  1:61745‑M1/251‑47 (MQ=255)                                                                                                                                                              aTTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgc                                                                                   >  2:39735‑M1/1‑205 (MQ=255)                                                                                                                                                                           tAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttac                                                                          <  1:371729‑M1/247‑56 (MQ=255)                                                                                                                                                                             aTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacg                                                                    >  1:9541‑M1/1‑190 (MQ=255)                                                                                                                                                                              ttGTTTAAAAAATGTGATTTTTATCACAAAGGCAATATGCCTGAGCAGAAGTCAGAGCCATAACTGGCACGTAGGGTTTGCAACCACAAACCCACCAATAGCGGGGTAGATAGGGCGTTAATCTCTCATACTTAACCTGGTTTATGGTACATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgcaggggttgtgatcatcacgtaatgcttattttaccggacgc                                                                   >  2:162452‑M1/1‑189 (MQ=255)                                                                                                                                                                               tGTTTAAAAAAAGTTATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACTTAGGGTTTTCAACCACTACCCCACCAATAGAGGGGTAGATAGGGCGGTAATCTCCCATACTTAACCTGGTTTATGTTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggtttggagcataacgtaatgcttattttaccggacgc                                                                   <  2:412365‑M1/250‑63 (MQ=255)                                                                                                                                                                                     aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccg                                                            >  1:216703‑M1/1‑182 (MQ=255)                                                                                                                                                                                     aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccg                                                            >  1:443122‑M1/1‑182 (MQ=255)                                                                                                                                                                                          aGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaac                                                       >  1:133227‑M1/1‑177 (MQ=255)                                                                                                                                                                                            tGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGAAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagaataacgtaatgcttattttacctgacgcatcgccgggaac                                                       >  2:152686‑M1/1‑175 (MQ=255)                                                                                                                                                                                             gATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacaga                                                    <  2:175264‑M1/251‑78 (MQ=255)                                                                                                                                                                                              aTTTTTATCCCAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgggggtgtgagcctaacagaatgcttatttgaccggacgcatcgcctggaacaga                                                    >  2:293‑M1/1‑173 (MQ=255)                                                                                                                                                                                              aTTTTTATCACAAAGGAACTATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaaccga                                                    >  2:414934‑M1/1‑173 (MQ=255)                                                                                                                                                                                                ttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACACACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGATTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaacaaaccccgacatttgccggggttgtgagcataaagtaatgcttattttaccggacgcatcgccgggaacagaat                                                  >  1:196554‑M1/1‑171 (MQ=255)                                                                                                                                                                                                   tATCACAAAGTAAATATGCCTGAGCAGCAGTAAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTGTATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggcgttgtgagcataacgtaatgcttattttaccggacgcatcgccggtaacagaataac                                               <  2:471725‑M1/250‑83 (MQ=255)                                                                                                                                                                                                    aTCACAAAGGAAATAGGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacg                                              <  1:218368‑M1/250‑84 (MQ=255)                                                                                                                                                                                                     tCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtc                                            >  1:410968‑M1/1‑166 (MQ=255)                                                                                                                                                                                                      cacaAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtc                                            <  1:118319‑M1/250‑86 (MQ=255)                                                                                                                                                                                                        caAAGTAAATATGCCTGAGCAGCAGGCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGTGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTTTAGATGATACGTTCtaaaaaaaccccgacatgtgccggggttgtgagcataacggaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcg                                         <  2:369410‑M1/251‑89 (MQ=255)                                                                                                                                                                                                          aaGGAAATATGACTGAGCAGCAGTCAGAGCCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGAAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGACCTCCATTCTGTTTAATTTGTAGAGGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgagg                                        >  2:118254‑M1/1‑161 (MQ=255)                                                                                                                                                                                                              aaaTATGCCTGAGCAGCAGTCAGAGACATAANTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggt                                   <  2:242286‑M1/251‑95 (MQ=255)                                                                                                                                                                                                              aaaTATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggt                                   >  1:31707‑M1/1‑157 (MQ=255)                                                                                                                                                                                                              aaaTATGCCTGAGCAGCAGTCAAAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacctttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgg                                    >  1:331104‑M1/1‑157 (MQ=255)                                                                                                                                                                                                              aaaTATGACTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggt                                   >  2:131873‑M1/1‑157 (MQ=255)                                                                                                                                                                                                                atatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtat                                 <  1:88996‑M1/251‑97 (MQ=255)                                                                                                                                                                                                                        gAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtgtggctgttg                         <  1:34796‑M1/251‑105 (MQ=255)                                                                                                                                                                                                                         agcagcAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggaatggctgttgg                        >  1:74543‑M1/1‑146 (MQ=255)                                                                                                                                                                                                                                gTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaaca                  <  1:400319‑M1/250‑112 (MQ=255)                                                                                                                                                                                                                                gTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggggaaca                  >  2:400319‑M1/1‑139 (MQ=255)                                                                                                                                                                                                                                 tCAGAGACATAACTGGCACGTAGGGATTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagc                >  1:253609‑M1/1‑138 (MQ=255)                                                                                                                                                                                                                                   agagACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgtgatggctgatggtgaaccgca               >  2:458817‑M1/1‑136 (MQ=255)                                                                                                                                                                                                                                     agaCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcatt             >  2:19828‑M1/1‑134 (MQ=255)                                                                                                                                                                                                                                         aTAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccat        >  1:181142‑M1/1‑130 (MQ=255)                                                                                                                                                                                                                                         aTAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccat        >  1:161931‑M1/1‑130 (MQ=255)                                                                                                                                                                                                                                               ggCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtc  <  1:175762‑M1/251‑128 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ATTAACGCTTTCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTT  >  NC_000913/2985480‑2985968

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