Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,904,907	C→T	E231K (GAA→AAA)	bglF ←	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,904,907	0	C	T	100.0%	146.1 / NA	46	E231K (GAA→AAA)	bglF	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/25);  total (21/25)

CACAGCGGCACCAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAG  >  NC_000913/3904660‑3905148                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  cacaGCGGCACCAGGCCCCAGTGCAGTCCGACAATGACGAAGATTTGCCAGAAGCCGCCCATCACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAg                                                                                                                                                                                                                                                <  2:481908/251‑1 (MQ=255)            cAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGcaca                                                                                                                                                                                                                                     <  2:461126/251‑1 (MQ=255)             aGGCCCCCGTGCAGTCCGAACATGCCGAAAATCTGCCAAAAGCCTCCCATTACCGCGCCAGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTTCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACaa                                                                                                                                                                                                                                    <  2:106379/251‑1 (MQ=255)                ccccAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATCACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATAAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAAcc                                                                                                                                                                                                                                  >  1:22975/1‑250 (MQ=255)                ccccAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCa                                                                                                                                                                                                                                 <  1:446580/251‑1 (MQ=255)                 cccAGTGAAGTCCGAAAATGACGAAGAGTTGCCAGAAGCCGCACATTACCGCGCCCGCAAAGGCAGTAACAGCTGGATAAAGCCAGAGATAACCGGCGGAAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGATGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAg                                                                                                                                                                                                                                <  2:202348/251‑1 (MQ=255)                  ccAGTGCAGGCCGAACATGACGAAGAGTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCa                                                                                                                                                                                                                                 <  2:430864/249‑1 (MQ=255)                                     aCGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAAt                                                                                                                                                                                                            >  2:524486/1‑251 (MQ=255)                                       gAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATa                                                                                                                                                                                                           >  2:555184/1‑250 (MQ=255)                                          gATTTGCCAGAAGCCGCCCATTCCCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCg                                                                                                                                                                                                       <  1:393748/251‑1 (MQ=255)                                                            cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTcaac                                                                                                                                                                                      >  2:239242/1‑250 (MQ=255)                                                              ttACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGAAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAAt                                                                                                                                                                                   >  1:228901/1‑251 (MQ=255)                                                               tACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAAt                                                                                                                                                                                   <  1:169800/250‑1 (MQ=255)                                                                           aaaTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATAc                                                                                                                                                                      <  1:50261/251‑1 (MQ=255)                                                                             aTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATAccc                                                                                                                                                                    >  2:99487/1‑251 (MQ=255)                                                                               gCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCa                                                                                                                                                                   <  1:242269/250‑1 (MQ=255)                                                                                cAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAgg                                                                                                                                                                 >  1:140889/1‑251 (MQ=255)                                                                                  ggAACCTCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATCCCCAGGaa                                                                                                                                                               <  2:209666/251‑1 (MQ=255)                                                                                     aCCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATc                                                                                                                                                            <  1:211847/251‑1 (MQ=255)                                                                                           tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGccc                                                                                                                                                      >  1:252564/1‑251 (MQ=255)                                                                                            gagaaaGACAGAGACAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTGTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATTCCCAGGAAATCCAGcc                                                                                                                                                       <  2:272029/246‑1 (MQ=255)                                                                                                     agagATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCagcc                                                                                                                                            >  1:504019/1‑248 (MQ=255)                                                                                                          tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCttt                                                                                                                                       >  1:522961/1‑251 (MQ=255)                                                                                                                   gcAATCAGGTCGTTTATCCAGGGTGATAGCGGCCCCACCATCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTGGAGTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCAAAACCATGCAGAAAAAATAATGGGAATACCCTATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCATCCCCAGCGCATCCGCTTTTTGCCCGtt                                                                                                                              <  2:929/251‑1 (MQ=255)                                                                                                                     aaTCAGTTCGATTATCCAGGTTAATAGCGGCCCCACAAGCCGAAAGGTGACGGGTGTGATAACCATCAGACCTAGCCCTGGTGTGAAGAAATTTTTGATTGCCGACGGTCACCACGCATTAAGCGGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATCACCGATGACGAGTACTTCAACAATGTGACCGGACTAACCAGGAAATCCAGCCCCAGAGCATCCGCTTTTGGCCCCTtct                                                                                                                            >  2:209831/1‑251 (MQ=255)                                                                                                                      aacAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAGTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAAAGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTtctc                                                                                                                           <  2:2435/249‑1 (MQ=255)                                                                                                                               cTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATAACCAGGAAATCCAGCCCCAGCGCATCAGCTTTTTGCCCGTTCTCGAAAGCAGt                                                                                                                  >  1:9686/1‑251 (MQ=255)                                                                                                                                  aTCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTc                                                                                                                 >  1:225373/1‑249 (MQ=255)                                                                                                                                     cAGGTTGATAGCGGCCCCACAAGCAGAAAGGTGACGGGTGCGATAACAATCAGACATAGCAATGGAGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGCTGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAAt                                                                                                            >  2:555973/1‑251 (MQ=255)                                                                                                                                          tGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATg                                                                                                       >  1:308175/1‑251 (MQ=255)                                                                                                                                              aGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTATTGGATg                                                                                                   <  2:267756/251‑1 (MQ=255)                                                                                                                                              aGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATg                                                                                                   >  1:197617/1‑251 (MQ=255)                                                                                                                                                cGGCACCACCAGCAGAAAGGTGACGGTTTTGACAACCATCAGACATAGCAAGGGTGTGAAGAAATTTTTGATGGCCGACGGTAACCACGCAGTAAGTCGGCGGTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGAGGACGAGGAATGCAACAATGTGACCGGAATACCCAGGAAATCCAGCACCAGCGCAGCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAGTTAAtgcatgca                                                                                                 <  2:22975/251‑1 (MQ=255)                                                                                                                                                 ggCCCCACCCGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGAGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACACCCAGGCAGAACAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCAGTTCTCGAAAGCAGTCAGATTTAATGGATGCAc                                                                                                >  2:265648/1‑251 (MQ=255)                                                                                                                                                    cccACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTaa                                                                                             <  1:22546/251‑1 (MQ=255)                                                                                                                                                                    tGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCAcc                                                                             >  1:429388/1‑251 (MQ=255)                                                                                                                                                                        gggTGTGATAACCATCAGACATAGCAATGGGGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAAGAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTAAGAATTAATGGATGCACTAACGCTCCACCAATCACCATgg                                                                         <  2:429388/251‑1 (MQ=255)                                                                                                                                                                          gtgtGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACg                                                                                                                                                                                                 <  1:537457/129‑1 (MQ=255)                                                                                                                                                                              gatcaCCATCAGACATAGAAATGGTGTGAAGAAATGTTTGATTGCCTACGGTAACCACTCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAAGCCAGCCCCAGCTCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTaa                                                                   <  2:51802/247‑1 (MQ=255)                                                                                                                                                                                           cATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccgcc                                                      >  1:76308/1‑251 (MQ=255)                                                                                                                                                                                                cAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAAGGGATTACCGCCGAAGc                                                 <  1:477325/251‑1 (MQ=255)                                                                                                                                                                                                              aaTTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGGGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCCCCAATCGCCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGgtg                                   <  1:507602/251‑1 (MQ=255)                                                                                                                                                                                                                     gatggCCGACGGTAACCACGCATTATGGCGGCGTTTCAGAATGAGGCACAACCAGGCAGAAAAAAGAATGGGAATATACGATGACTAGTAATGCAACAATGTGACCGGAATCCCCAGGAAATCCATCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGAGGCACTAACGCTCCACCAATCCCCATGTCAGTAAATGGATTCCCGCCGAAGCGTTTCCCCGCGGTGTATccc                             <  2:168401/246‑1 (MQ=255)                                                                                                                                                                                                                     gATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGACAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGAAATACACAGGAAATCCGGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGCATTAATGGATGCACTAACGCTCCACCAGTCACCATGACAGTAAATGGAATACCGCGGAAGCGTCTCCCCGCGGTGTATccc                             >  1:416459/1‑250 (MQ=255)                                                                                                                                                                                                                                    ccACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGga                >  1:60349/1‑248 (MQ=255)                                                                                                                                                                                                                                               gTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAAGCCAGCCCCAGCGCATCCGCTTTTTTCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAg  <  1:555184/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  CACAGCGGCACCAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAG  >  NC_000913/3904660‑3905148

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