Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,139,511	Δ1 bp	coding (198/663 nt)	fsaB ←	fructose‑6‑phosphate aldolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,139,511	0	C	.	100.0%	225.8 / NA	53	coding (198/663 nt)	fsaB	fructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (25/28);  total (25/28)

GAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGG  >  NC_000913/4139270‑4139751                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  gAACAGTACGAAAGCCGGTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATGTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGGATTTTAAGTGCTGCCAGACCGTCGGAAGTCACCGGGATGTGCACCCCAATACCCGGAATAGCGTCGCGCAGGCGCGTCGCTTCTTCCACCATCCC‑TGCGCGTcgc                                                                                                                                                                                                                                         <  2:505208/251‑1 (MQ=255)     cGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCt                                                                                                                                                                                                                                     <  1:422192/251‑1 (MQ=255)        aCGAATGCCGGTTCCTCCCTGGGCATCTACGCGTTGAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCAt                                                                                                                                                                                                                                  <  2:504131/251‑1 (MQ=255)          gAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATgg                                                                                                                                                                                                                                >  1:120191/1‑251 (MQ=255)           aaTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGt                                                                                                                                                                                                                               >  1:180380/1‑251 (MQ=255)                 gtttcCGCCCTGCGCATCTACGCGGTTAACATCCGGAGCAAAGTATTTTGCCCCTGCCAGTGCGGATAATACCCATTGTGCATCGCTATATACAGCAGGGCCAAGAGTAGAAATACCCTCTTTTTTACGTATTTTAATTGCTGACAGACCTTCGGAAGTCACCGGGATTTTCACCACAATCACCGGAAAAGCGTCGTGCAGACGCTTCGCTTCTTCCACCATGCC‑TGCGCGTCGCGGCTCATAGTCTg                                                                                                                                                                                                                            >  2:509624/4‑248 (MQ=255)                     ccGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCNC‑TGCGCGTCGCGGCACATGGTCTGAGCAAAc                                                                                                                                                                                                                     >  2:248992/1‑251 (MQ=255)                       gCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAg                                                                                                                                                                                                                   <  1:414438/251‑1 (MQ=255)                                    aCGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGGCTGAGCAAACAGAAGGCCCtcatcc                                                                                                                                                                                                      >  2:439996/1‑250 (MQ=255)                                                  cGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAAGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTg                                                                                                                                                                                        <  1:111892/251‑1 (MQ=255)                                                    gAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGc                                                                                                                                                                                       >  2:151351/1‑250 (MQ=255)                                                           gTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCCAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCgg                                                                                                                                                                               >  1:378634/1‑251 (MQ=255)                                                              ttttGCCACTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGTGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑GGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGAcgcgtca                                                                                                                                                                             <  2:169938/250‑4 (MQ=255)                                                               tttGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGgcag                                                                                                                                                                            <  2:85251/250‑1 (MQ=255)                                                                     ccTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACATCAGTGCCAAGTGTAGGAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGTATTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACCGAATGCCCTCATCCCCAATCGCTTTTTGTAGCCGCGGCAGCACTTcc                                                                                                                                                                     >  2:49888/1‑251 (MQ=255)                                                                                 gCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTcc                                                                                                                                                         <  1:71006/251‑1 (MQ=255)                                                                                     ataaCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGTCCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGc                                                                                                                                                     >  1:201037/1‑251 (MQ=255)                                                                                          ccTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGNCAGCACTTCCCATATGGCCTCCTTGCTGGc                                                                                                                                                 >  2:257761/1‑250 (MQ=255)                                                                                             tgtgCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATTCCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGc                                                                                                                                              <  1:307966/250‑1 (MQ=255)                                                                                                gcggcgCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGATTCCTTGCTGGCAGCGAt                                                                                                                                           <  1:149236/250‑1 (MQ=255)                                                                                                    cgcTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCCTGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGGCGCGGCAGCACTTCCCATATGGACTCCTTGCTGCCAGCGATAATGc                                                                                                                                      >  2:519236/1‑251 (MQ=255)                                                                                                            aCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGtt                                                                                                                               <  1:242293/250‑1 (MQ=255)                                                                                                              agcagTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCCCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGttcgt                                                                                                                            >  1:283374/1‑248 (MQ=255)                                                                                                                     ccAAGTGTAGTAATGCCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATTCCCGGAATAGCGTCGCGCAGGCTCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCAcc                                                                                                                     <  2:417400/251‑1 (MQ=255)                                                                                                                       aaGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCAcccc                                                                                                                   >  2:223546/1‑251 (MQ=255)                                                                                                                             agtactaccaactTTTTTCAGAATTTTAATTGCTGCCAGACATTCGGAAATCACCGGGATTTTCACCCCAATACCAGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCCGCCGCGGCAGCACTTCCCATATGGACACCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACACCGGTAAt                                                                                                             >  2:494566/12‑251 (MQ=255)                                                                                                                                  tACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGCAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATgggg                                                                                                         >  2:147943/1‑250 (MQ=255)                                                                                                                                                    ttttAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGtt                                                                                      >  1:542535/1‑251 (MQ=255)                                                                                                                                                        aaTTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GAc                                                                                 >  1:251300/1‑251 (MQ=255)                                                                                                                                                               ggcaggccTTCGGAAGGCACCGGGATTTACCCCACAATACCCGGAAGAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGTGTCGCGTCTCTTGGTCTGAGCAAACAGAATGCCCTCTTCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCTGCGATAATGCTCGGGTTAGTTGTCTCCCCGGCAATGGGGAATAAGCGTGCCAGACGTTT‑GACTTCTgt                                                                           <  2:310238/244‑2 (MQ=255)                                                                                                                                                                 cAGACCTTCGGAAGTAACAGGGATTTTCACCNCAATACCCGGAATAGCGTCGAGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCGCATGGTCTGAGCAAACAGAATGACCGCAGTACCAATCGCGGTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCGGCGAc                                                                        <  2:253792/251‑1 (MQ=255)                                                                                                                                                                   gACCTTGGGAATTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGAGCGTCGCGGCTCATGGTCGGAGCAAACAAAATGCCCTCACCCCCATTCGCTTTGTGCAGACGCGGCAGCACTTCCCATATGGACTCCGTGCTGGCAGCGATAACGCTCGGGTTAGTTGTCACCACGGCAATGGGGCATATTCGTGCCAGACGTTC‑GACTTcggagccgt                                                                      >  2:358597/1‑243 (MQ=255)                                                                                                                                                                    aCCTTTGGCAGTCCCCGGGTTTTTCCCCCCAATGCCCGGAATAGCGTCCCGCAGGCGCTGCGCTGCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGTTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTGTTTGCAGACGCGGCCGCGCTTCCCATTTGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGTGGAATATGCGTGCCAGACGTTC‑TACTTCTGCGACGtt                                                                     <  2:293378/251‑1 (MQ=255)                                                                                                                                                                      cTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGtt                                                                     >  2:88767/1‑249 (MQ=255)                                                                                                                                                                         cGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGCGCTAACAGAATGCCCTCATCACCAATCTCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGATGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATACGTGTGCCCGTCGTTCTGA‑TTCTGCGACGTTGGCgg                                                                >  2:361904/1‑251 (MQ=255)                                                                                                                                                                           gAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGgtg                                                              <  1:423776/251‑1 (MQ=255)                                                                                                                                                                              gTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTcc                                                           <  1:220657/251‑1 (MQ=255)                                                                                                                                                                              gTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTcc                                                           >  1:374703/1‑251 (MQ=255)                                                                                                                                                                                   cGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCTATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGGAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATa                                                      <  1:240533/251‑1 (MQ=255)                                                                                                                                                                                   cGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATa                                                      <  1:448500/251‑1 (MQ=255)                                                                                                                                                                                    gggATTTTCACCCCAATACCCGGAATAGCGTCGCGCAGGCTCTTCGCTTCTCCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATAc                                                     <  2:374703/251‑1 (MQ=255)                                                                                                                                                                                                  aaTACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtt                                       >  2:91176/1‑251 (MQ=255)                                                                                                                                                                                                   aTACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGttt                                      <  1:151351/251‑1 (MQ=255)                                                                                                                                                                                                     ccccGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTcc                                    <  2:120191/250‑1 (MQ=255)                                                                                                                                                                                                      cccGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCt                                   <  1:12563/251‑1 (MQ=255)                                                                                                                                                                                                         ggAATAGCGTCGCGCAGGCGCTTCGCTTCTTNCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtttt                                <  2:251282/251‑1 (MQ=255)                                                                                                                                                                                                              agcttCGCGCAGGCGCTTCGCTTTATCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCAACACCAATCGCTTTTTGCAGACGCGGCAGCTCTTCCCATATGGACTCCTTGCTGGCAGCGTTAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGTCGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTGCCATGATGTTTCCTTTTAcagt                           <  2:210838/247‑1 (MQ=255)                                                                                                                                                                                                              aGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGGGCCAGACGTTC‑GACTTCTGCGACGTTAGAGGTGTCCAGATACAGTTCCATGCTGTTTCCTTTTAcagt                           >  2:59083/1‑251 (MQ=255)                                                                                                                                                                                                                   cgcgCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCAGCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagttccag                      <  2:291039/251‑1 (MQ=255)                                                                                                                                                                                                                        aGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCAt                 <  1:433861/251‑1 (MQ=255)                                                                                                                                                                                                                                   ttcttcCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAgg      >  1:241259/1‑251 (MQ=255)                                                                                                                                                                                                                                   ttcttcCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAgg      >  1:159093/1‑251 (MQ=255)                                                                                                                                                                                                                                        ccaccaTCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTgg  <  1:257761/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  GAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTC‑GACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGG  >  NC_000913/4139270‑4139751

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