Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,985,816	T→C	intergenic (‑224/‑31)	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,816	0	T	C	100.0%	294.3 / NA	77	intergenic (‑224/‑31)	yqeF/yqeG	putative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/36);  total (41/36)

GACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAG  >  NC_000913/2985572‑2986060                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     gaccaTCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCCCTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGGAAATTGCCCTCCATTCTGTTTa                                                                                                                                                                                                                                                <  2:33426/247‑1 (MQ=255)    aaTCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAAtt                                                                                                                                                                                                                                             >  1:101309/1‑251 (MQ=255)     aTCACAACGTCTTTCATCACAACATCGCACCNATTAACAGTGTAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGTATCAGTCAGAGACATAACTGGCACGTTGGGTTTGCAACCACTTACCCACCGATAGAGGGGTAGATAGGTCGTTAATCTCCCATACTTAACCTGGTTTAGGGTAAATTGCCCTCCGTTCTGTTTAAttt                                                                                                                                                                                                                                            <  2:356095/251‑1 (MQ=255)                      cacaACCTCGCACCAATTAACAGTGCAATAGNAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCta                                                                                                                                                                                                                           <  2:343351‑M1/251‑3 (MQ=255)                          aCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTATTTTGTAgat                                                                                                                                                                                                                                      >  2:474316/1‑236 (MQ=255)                          aCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAgat                                                                                                                                                                                                                                      <  1:474316/236‑1 (MQ=255)                          aCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATGCTGTTTAATTTGTCGATGATACGTTCtaaaaa                                                                                                                                                                                                                       <  2:628593‑M1/251‑7 (MQ=255)                              cGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaacc                                                                                                                                                                                                                   <  1:203418‑M1/251‑11 (MQ=255)                                cACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACAAATAGAGGGGTAGAGAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATACTGTTGAATTTGTAGATGATAAGTTCtaacaaaacccc                                                                                                                                                                                                                 >  1:654253‑M1/1‑239 (MQ=255)                                 aCCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccg                                                                                                                                                                                                                >  2:421312‑M1/1‑238 (MQ=255)                                   cAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                                                                                                              >  2:633552‑M1/1‑236 (MQ=255)                                    aaTTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                                                                                                              >  1:211650‑M1/1‑235 (MQ=255)                                     aTTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGAGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGTTACGTTCtaaaaaaaccccgacat                                                                                                                                                                                                            >  1:454560‑M1/1‑234 (MQ=255)                                       tAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacattt                                                                                                                                                                                                          >  1:707678‑M1/1‑232 (MQ=255)                                        aaCAGTGCAATAGGAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGAAATAACTGGCACGTAGGGTTTGCAACCTCTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAGCCTGGTTTATTTTAAATTGCCCTCCATTCTGTTTAAGTTGTAGATGATACGTTCtaaaaaaaccccggcattt                                                                                                                                                                                                          <  2:238528‑M1/250‑20 (MQ=255)                                                aaTAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggtt                                                                                                                                                                                                 >  1:163551‑M1/1‑223 (MQ=255)                                                aaTAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGAGATTTTTATCACAAAGGAAATATGACTGAGCAGCAGACAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaataaaccccgacatttgccggggtt                                                                                                                                                                                                 >  2:130887‑M1/1‑223 (MQ=255)                                                           aTTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataac                                                                                                                                                                                      >  2:406984‑M1/1‑212 (MQ=255)                                                                  aTTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatg                                                                                                                                                                                >  1:627357‑M1/1‑205 (MQ=255)                                                                  aTTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCATCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgc                                                                                                                                                                               >  1:167923‑M1/1‑205 (MQ=255)                                                                        gTTAACTTAATTGTTTAAAAAAAGTGATTCTTATCACGAAGGAAATATGCCTGAGCAGCATTCAGAGCCATAACTGGAACGTAGGTTTTGCACCCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCGCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaataaaccccgacattggccggggttgtgaggataacgttattcttattt                                                                                                                                                                         <  2:602912‑M1/251‑53 (MQ=255)                                                                         ttATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgcnggggttgtgagcataacgtaatgcttatttt                                                                                                                                                                        >  2:377669‑M1/1‑198 (MQ=255)                                                                              ttatttgtttAAAAAAAGTGATTTTTATCACAAAGGAAATATGACTGAGCAGCAGTCAGAGACATAACTGGCACGTAGAGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGTGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCACTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgaccttttccggggttgtgagcataacgtaatgcttattttaccgt                                                                                                                                                                   >  2:557281‑M1/5‑193 (MQ=255)                                                                               tAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttac                                                                                                                                                                      <  2:126532‑M1/247‑56 (MQ=255)                                                                                aaTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGTGACATAACTGGTACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttagtttaccggac                                                                                                                                                                 <  1:129302‑M1/251‑61 (MQ=255)                                                                                 aTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATAc                                                                                                                                                                                                                                 >  2:392353/1‑186 (MQ=255)                                                                                 aTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATAc                                                                                                                                                                                                                                 <  1:392353/186‑1 (MQ=255)                                                                                       tAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgc                                                                                                                                                          >  2:164644‑M1/1‑184 (MQ=255)                                                                                                  aTTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacaga                                                                                                                                                >  1:139519‑M1/1‑173 (MQ=255)                                                                                                    atttaTCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaat                                                                                                                                              <  1:583028‑M1/249‑80 (MQ=255)                                                                                                        aTCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgt                                                                                                                                         >  2:45409‑M1/1‑167 (MQ=255)                                                                                                          cacaAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcg                                                                                                                                       >  1:179949‑M1/1‑165 (MQ=255)                                                                                                          cacaAAGGAAAGATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtc                                                                                                                                        <  2:555019‑M1/250‑86 (MQ=255)                                                                                                               aGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggat                                                                                                                                  <  1:164644‑M1/251‑92 (MQ=255)                                                                                                               aGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTGTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGTTTTATGGTAAATGGCCCTCCATTCTGTTTAATTTGTAGATGATACGGTCtaaaaaaaccccgacagttgccggggttgggagcataacgtaatgcttatgttaccggacggaacgccggtaacagaataacgtcgcggat                                                                                                                                  <  2:329247‑M1/251‑92 (MQ=255)                                                                                                                 gAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgg                                                                                                                                >  2:96782‑M1/1‑158 (MQ=255)                                                                                                                  aaaTATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggt                                                                                                                               <  1:406984‑M1/251‑95 (MQ=255)                                                                                                                   aaTATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaacccagacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacggcgcggatggtt                                                                                                                              >  2:128541‑M1/1‑156 (MQ=255)                                                                                                                      atGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatgg                                                                                                                           >  1:166297‑M1/1‑153 (MQ=255)                                                                                                                           tGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgctggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgtt                                                                                                                      >  1:368466‑M1/1‑148 (MQ=255)                                                                                                                               cagcagTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtg                                                                                                                  >  1:441439‑M1/1‑144 (MQ=255)                                                                                                                                  cagTCAGAGACAGAACTGGCTCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATCGGGCGTTAATCACCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGGTTAAGTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggtggtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggatgttggtgatc                                                                                                               <  2:616062‑M1/251‑111 (MQ=255)                                                                                                                                      cAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaaccgc                                                                                                            >  1:749782‑M1/1‑137 (MQ=255)                                                                                                                                       agagACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcat                                                                                                          <  1:647699‑M1/251‑116 (MQ=255)                                                                                                                                       agagACATAACTGGCACGTAGGGATTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggggcaacgcct                                                                                                          >  2:583326‑M1/1‑136 (MQ=255)                                                                                                                                        gagaCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcatt                                                                                                         <  1:562058‑M1/251‑117 (MQ=255)                                                                                                                                          gaCATAACTGGCACGTAGGGTTTGAAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATATCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagcataacgtcgcggatggtatggctgttggtgaacagcattac                                                                                                       >  1:619643‑M1/1‑133 (MQ=255)                                                                                                                                                   ggCACGTAGGGTTTGCAACCACTAACCCANCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATGCTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacgg                                                                                                <  2:334700‑M1/249‑126 (MQ=255)                                                                                                                                                   ggCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtc                                                                                              <  1:421312‑M1/251‑128 (MQ=255)                                                                                                                                                     cACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaatccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggac                                                                                                                                                                 <  1:362246‑M1/182‑61 (MQ=255)                                                                                                                                                     cACGTAGGGTTTGCAACAACTAACCCACCAATAGAGGGGTAGATAGGGCGTTACTCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggac                                                                                                                                                                 >  2:362246‑M1/1‑122 (MQ=255)                                                                                                                                                      aCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgat                                                                                           >  1:75368‑M1/1‑121 (MQ=255)                                                                                                                                                               ttGCAACCTCTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTTCCATACTTAACCTGGTTTACGGTAAATTGCCCTCCATTCTGTTGAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaaggcgtattttaccggacgcatcgccgggaacagaataacggcgcggagggtatggctgttggtgaacagcattaccatacggtcgataccaagtc                                                                                   <  2:196705‑M1/250‑139 (MQ=255)                                                                                                                                                                 gCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATTCTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggtgttgtgagcataacgtaatgcttatttttccggacgcatcgccgggatcagaataacgtcgcggatggtatgtctgttgttgaacagcattaccatacggtcgataccaattccc                                                                                 <  1:569‑M1/250‑141 (MQ=255)                                                                                                                                                                 gCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattccca                                                                                <  1:624987‑M1/251‑142 (MQ=255)                                                                                                                                                                 gCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattccc                                                                                 <  1:598483‑M1/250‑141 (MQ=255)                                                                                                                                                                   aaCCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattccc                                                                                 >  1:233963‑M1/1‑108 (MQ=255)                                                                                                                                                                              caccaATAGAGGGGTAGATAGGGCGTTAATCNCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcgg                                                                   <  2:368466‑M1/251‑155 (MQ=255)                                                                                                                                                                                     agagGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTCTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggggaacagcattaccatacggacgataccaattcccagacctgctgtcggcggga                                                              >  2:83086‑M1/1‑90 (MQ=255)                                                                                                                                                                                      gagGGGTAGATAGGGCGTTAATCTCCCATGCTGAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATGTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttatttaaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagc                                                           <  2:131773‑M1/251‑163 (MQ=255)                                                                                                                                                                                       agGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggcacagaataacgtcgcggatggtatggctgttggtgaacagcattaacatacggtcgataccaattaccagacctgctgtcggcgggaagcc                                                          >  2:513455‑M1/1‑88 (MQ=255)                                                                                                                                                                                              gATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttcc                                                   <  2:684994‑M1/251‑171 (MQ=255)                                                                                                                                                                                               aTAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttcca                                                  <  2:490049‑M1/251‑172 (MQ=255)                                                                                                                                                                                                 aGGGCGTTAATCTCCCATTCTGAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccat                                                        <  2:25079‑M1/243‑166 (MQ=255)                                                                                                                                                                                                 aGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccat                                                        >  1:25079‑M1/1‑78 (MQ=255)                                                                                                                                                                                                   ggCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgc                                              >  1:288949‑M1/1‑76 (MQ=255)                                                                                                                                                                                                         aaTCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagcgcggtgac                                        >  2:307382‑M1/1‑70 (MQ=255)                                                                                                                                                                                                                     ttAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGCTGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggggaacagcattaccatacggtcgataccaattcccagacctcctgaccgcggtaaaccaagttccaatgcgagcacgaaaattttaac                            >  2:461757‑M1/1‑58 (MQ=255)                                                                                                                                                                                                                         ccTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggcacagaataacgtcgcggatggtatggctgttggtgaacagcattacaatccggtcgataccaattccccgacctgctgtcggcggtaagccaanttccagtgcggtgaagtaatcttcatcca                          >  2:369230‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                                 aTGGTAAATTGCCCTCCATTCTGGTTAAGTTGTAGATTATCCGTTCtaaaaaaaccccgacaggtgccggggttgtgagcataacgtaatgcttattttaccgtacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgatgccacttcccagacctgctgtcggcggtaagccatggttcagtgcggtgacgtaatcttcatcgtagagcatcgc                <  2:194213‑M1/251‑206 (MQ=255)                                                                                                                                                                                                                                   ggTAAATTGCCCTCCATTCTGTTTAATTTTGAGATGATACGTTCtaaaaaaaccccgacatttgcctgggttgtgagcataacgtaatgcttattttacctgacgcatcgccgggaacagaataacgtcgcggatggtatgtctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcggataacatcgctt              <  1:571267‑M1/251‑208 (MQ=255)                                                                                                                                                                                                                                   ggTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgct               <  1:716008‑M1/250‑207 (MQ=255)                                                                                                                                                                                                                                    gTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcggagaacatcgcttc             <  2:101309‑M1/251‑209 (MQ=255)                                                                                                                                                                                                                                     tAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttcg            >  1:453487‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                                       aaTTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttcgtc          >  1:251731‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                                               tCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgttgggggtaagcccttgtccagtgcggtgacgtaatcttcatcgtagaacatcgcctcgtcgtcacctg  >  2:287777‑M1/1‑32 (MQ=255)                                                                                                                                                                                                                                                ccATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttcgtcgtcacct   <  2:75368‑M1/249‑219 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAG  >  NC_000913/2985572‑2986060

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