Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,279,875	Δ1 bp	intergenic (+47/‑105)	ghxP → / → yjcE	guanine/hypoxanthine transporter GhxP/putative transporter YjcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,279,875	0	A	.	100.0%	325.8 / NA	71	intergenic (+47/‑105)	ghxP/yjcE	guanine/hypoxanthine transporter GhxP/putative transporter YjcE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (32/39);  total (32/39)

ATTTTGTTGATGCGATGGCGGGTCTGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGT  >  NC_000913/4279629‑4280120                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       aTTTTGTTGATGCGATTGCGGGTCTGGGTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGGCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTTGCGCTGGTGACCTTCTATGCGGGTGGCTGTGCTATCTAATCTTTCATCGGATTCTGAAATCGGGTGGCAATTGCTGCCTGTTTTT‑TTTtc                                                                                                                                                                                                                                                   <  2:449303/251‑1 (MQ=255)   tttGTTGATGCGATGGCGGGTCTGGTTTCGGTTGTATTCATCGTGCTGACCTGTAACATCTTAACAGGCATCATGATCGGCTTCGCGACGCTGGTGATTGGGCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCAGCGGATTCTGAAAACGGGGGGCAATGGCTGCCCGTTTTT‑TTTtct                                                                                                                                                                                                                                                  <  2:125316/250‑1 (MQ=255)         gATGCGATGGCGGGTCTGGTTACGGCTGTATTCATCGTGCTGACCTGTTACATCGTAACATGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGGTGCAATGGCTGCCCGTTTTT‑TTTTCTCCGcac                                                                                                                                                                                                                                            <  1:151136/250‑1 (MQ=255)         gATGCGATGGCGGGTCTGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACCGGCATCCTGCTCGGCTTCGCGCCTCTGGTGCTTGGTCGTCTGGTTTCCGGCGAATGGCGCACGTTGAACCTCGGTCCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCCTCGGATTCTGAAAACGGGTGGCCCTGGCTGCCCGTTTTT‑TTTTTTCCGcacc                                                                                                                                                                                                                                           >  1:115740/1‑250 (MQ=255)         gATGCGATGGCGGGTCTGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGc                                                                                                                                                                                                                                              <  2:324757/248‑1 (MQ=255)         gATGCGATGGCGGGTCTGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGc                                                                                                                                                                                                                                              >  1:324757/1‑248 (MQ=255)         gATGCGATGGCGGGGCAGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTTGCGAATCTGGGGATTGGTCGTCTGGGTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGGGACCTTCTATGCGGGTGGCTGGGCTATCTAATTTTTCATCGGATTCTGAAAACGGGTGGCCATGGCTGCCCGTTTTT‑TTTTTTCCGcacc                                                                                                                                                                                                                                           >  1:344080/1‑250 (MQ=255)                tGGCGGGTCTGGGTACGGCTGTATTCATCGTGCTGGCCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGGGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGGACGGTCGTTTTCGCCGTGGCGCGGGTGGCCTTCTATGCGGGGGGGCGGGCTATCTAATCTTTCATCGGGTTCTGAAAACGGGGGGCAATGGCTGCCCGtttttttttttccccccattttgtt                                                                                                                                                                                                                                     >  2:14728/1‑236 (MQ=255)                   gggggcTGGTTCCGGCGGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCAGCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGAGTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGc                                                                                                                                                                                                                                 <  2:154452/246‑1 (MQ=255)                    gggTCTGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGTTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTTTGAGCt                                                                                                                                                                                                                                <  2:313614/251‑1 (MQ=255)                     gggcTGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCTCACATTGTGTGAGCtt                                                                                                                                                                                                                               <  1:570412/248‑1 (MQ=255)                       tCTGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCtttt                                                                                                                                                                                                                             >  1:496851/1‑251 (MQ=255)                          ggTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                                                                                                                                                                          <  2:86957/251‑1 (MQ=255)                              aCGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGTCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑GGTTCTCCGCACATTGTGTGAGCTTTTGCGata                                                                                                                                                                                                                       <  1:612030/250‑1 (MQ=255)                                ggCTGCATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGCTCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGGATGGCGCAAGATGAACATCGGTACGGACGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCAACGGATTCTGAAAACGGGGGGAAATGGCTGCCCGTTTTT‑TTTTTTCCCCACGTTGTGTGAGCTTTTGCGATATTc                                                                                                                                                                                                                    >  2:462740/1‑251 (MQ=255)                                 gCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCTGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGGGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCt                                                                                                                                                                                                                   >  1:384973/1‑251 (MQ=255)                                    gTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCTGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGTTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaa                                                                                                                                                                                                                <  1:636957/251‑1 (MQ=255)                                     tATTCCTCGTGATGACCTGGAACATCGTAACAGGCATCATGCTCGGCTTCGAGCCTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGATGAACATCGGTACGGTCGTTATCGCAGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGACTGCCCGTTTTT‑TTTTCTCCGCACATTGTGGGAGCTTTTGCGTTATTTTGaa                                                                                                                                                                                                                >  2:700807/1‑250 (MQ=255)                                          aTCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGCTTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAAt                                                                                                                                                                                                          >  1:17405/1‑251 (MQ=255)                                            cGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATga                                                                                                                                                                                                        <  1:693179/251‑1 (MQ=255)                                            cGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATg                                                                                                                                                                                                         >  2:294292/1‑250 (MQ=255)                                            cGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATg                                                                                                                                                                                                         >  1:549027/1‑250 (MQ=255)                                                 tGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAAtt                                                                                                                                                                                                   >  2:300028/1‑251 (MQ=255)                                                       gTAACATCGTAACAGGCATCTGGATCGTCTTCTCGACTCTTGTGATTGGTCGTCTTGTTTCCGGCGAATGGCGCAAGTTGAACAACTGTACGGTCGTTATCGCCTTGGCGCTTGTGACCTTCTATGCGGGTGGCTGAGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGTCTGCCCGTTTTT‑TTTACTCCGCACATTGTGTGAGCTTTTGCGATATTCTTAAAAAAATGAGAATTCAg                                                                                                                                                                                                <  2:587778/248‑1 (MQ=255)                                                       gTAACAGCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATAGTCTGAAAAAAATGAGAATTCAgg                                                                                                                                                                                               <  2:472807/249‑1 (MQ=255)                                                                 aaCAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCttt                                                                                                                                                                                   >  1:184347/1‑251 (MQ=255)                                                                      gCATCATTATCGGCTTCGCGTCTCTGGTGATTGGTCGTCTGGTGTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTTGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAg                                                                                                                                                                              <  1:684724/251‑1 (MQ=255)                                                                       catcatGGTCGGCGTAGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGGGAATGGCGCAAGTTGAACATCGGTTCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCGTTCTCAGa                                                                                                                                                                             <  2:155404/251‑1 (MQ=255)                                                                        atcatGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtctc                                                                                                                                                                                >  1:40623/1‑247 (MQ=255)                                                                          catGATCGGCTTCGCGACTTTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGATCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATc                                                                                                                                                                           <  2:403060/250‑1 (MQ=255)                                                                                cGGCTTCGCGACTCTGGTGATTGGTCGGCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGGTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaaaa                                                                                                                                                                    <  1:219217/251‑1 (MQ=255)                                                                                     tCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACa                                                                                                                                                                 <  2:184347/249‑1 (MQ=255)                                                                                       ccgACGCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACAGCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAATAAACATCGc                                                                                                                                                             <  2:119682/250‑1 (MQ=255)                                                                                        cgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCa                                                                                                                                                            <  2:201509/251‑1 (MQ=255)                                                                                          aCTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCaa                                                                                                                                                           >  2:590452/1‑250 (MQ=255)                                                                                               ggTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAg                                                                                                                                                       <  1:381201/249‑1 (MQ=255)                                                                                               ggTGATTGGTCGTCTGGTTGCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGTGCTGGTGACCTTCGATGCGGGTGGCTGGGCTAACTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGAAAAccgg                                                                                                                                                      <  2:17405/250‑4 (MQ=255)                                                                                                     ggggagtctggggtCCGGCGAATGGCGCAAGTTGAACATCGGGACGGTCGGGATCGCCGTGGCGCGGGTGACCTTCTAGGCGGGTGGCTGGGCTATCTAATCTGTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGGGAGCGTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTGTCGCAGATCATAAAAAACATCGCAAACAGGGAACgaa                                                                                                                                               <  2:115740/238‑3 (MQ=255)                                                                                                          gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGa                                                                                                                                           >  1:25677/1‑250 (MQ=255)                                                                                                              ggTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCAGGGaaa                                                                                                                                         <  1:580223/248‑1 (MQ=255)                                                                                                                         aaTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCAt                                                                                                                           >  1:324834/1‑250 (MQ=255)                                                                                                                          aTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATAc                                                                                                                         >  1:35182/1‑251 (MQ=255)                                                                                                                           tGGCGCACGTTGACCATCGGTACGGTCGTTATAGCCCTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGCAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGACAACAATGAGAATTCAGGCATAACGTCTTTCTCAGTTCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CCCCATACt                                                                                                                        >  2:704145/1‑251 (MQ=255)                                                                                                                              cgcAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGAt                                                                                                                     <  1:473775/251‑1 (MQ=255)                                                                                                                                  aGTTGAACATCGGTTCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATg                                                                                                                 <  2:537795/251‑1 (MQ=255)                                                                                                                                     tGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGAc                                                                                                               <  2:85094/250‑1 (MQ=255)                                                                                                                                          aTCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGAccgtacc                                                                                                         >  1:173695/1‑246 (MQ=255)                                                                                                                                             ggTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGAccgtaccgat                                                                                                      <  2:173695/251‑9 (MQ=255)                                                                                                                                                  ggTCGTTATCGCCGTGGCGCTGGTGGCCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAACAACATCGCAAACAGGAAACGAATGGAACTCTTCTT‑CACCATACTGATAATGACCCTCg                                                                                                          >  2:63668/1‑242 (MQ=255)                                                                                                                                                  ggTCGTTAGCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCg                                                                                                          <  1:63668/242‑1 (MQ=255)                                                                                                                                                        tATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTT‑GGGGTTTTTTGGAAAAAATTAGAATTCAGGGATTAAGGCCTTCTTAGATCATTAAAAAAATCGCAAACAAGGAACCCATGGAAATCTTCTTTAACCATAC‑GATAATTACCCTccgggtcttcgttgtccg                                                                                          >  1:649101/1‑234 (MQ=255)                                                                                                                                                             ccGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGggtg                                                                                      <  2:676419/251‑1 (MQ=255)                                                                                                                                                                               tCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATTGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACTCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGc                                                                      <  1:437713/249‑1 (MQ=255)                                                                                                                                                                                cTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGTGTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTTCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATACAAAACATCGCAAACAGAGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGCTCACTCGTGTCATGcact                                                                   >  1:69696/1‑248 (MQ=255)                                                                                                                                                                                  aTGCGGGTGGCTGGGCTATCTACTCTTTCATCGGATTCTGCAAACGGGTGGCCATGGCTGACCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGCCCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCttt                                                                 >  2:455966/1‑251 (MQ=255)                                                                                                                                                                                       ggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGAt                                                            <  2:542404/251‑1 (MQ=255)                                                                                                                                                                                               ggCTAACTAATCTTTCATCCGATTCTGAAAACGGGTGGCAATGGACGCCCGTTTTT‑TTTTCTCCGCAACTTGCGTGAGCTTTTGCGATATTCTGAAAAACATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGACTGGAAATCTTCTT‑CACAATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCccgcttc                                                    >  2:707549/1‑251 (MQ=255)                                                                                                                                                                                                gCTATCTTATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CTCCATACTGATAATGACCCTCGTTGTCTCGCGGTCCGGGGGGGTTACTCGTTTCATTCCCTGTCAGATCccgcttcc                                                   <  2:407039/251‑1 (MQ=255)                                                                                                                                                                                                        aTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGc                                           >  2:177236/1‑251 (MQ=255)                                                                                                                                                                                                             tCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATGCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCTTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAAt                                       <  2:384973/250‑1 (MQ=255)                                                                                                                                                                                                              cATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTAAGGCCTAACGTATTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCg                                     >  1:98528/1‑251 (MQ=255)                                                                                                                                                                                                                     ttCTGAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGAGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGACATCTTCTT‑CACCATACTGATAATGACACTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCgg                              >  2:614099/1‑251 (MQ=255)                                                                                                                                                                                                                         gAAAACGGGTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTgc                           <  1:294292/250‑1 (MQ=255)                                                                                                                                                                                                                                ggTGGCAATGGCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTgg                   >  1:179624/1‑251 (MQ=255)                                                                                                                                                                                                                                        gggcggCCCGTGTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATCATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGcc            <  2:582413/245‑1 (MQ=255)                                                                                                                                                                                                                                         ggCTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCgc          >  2:381377/1‑250 (MQ=255)                                                                                                                                                                                                                                          gCTGCCCGGTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGAc         <  1:462740/251‑1 (MQ=255)                                                                                                                                                                                                                                           cTGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCgccg        >  2:354830/1‑248 (MQ=255)                                                                                                                                                                                                                                            tGCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACgt       >  1:581957/1‑251 (MQ=255)                                                                                                                                                                                                                                             gCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACgtt      <  1:14728/251‑1 (MQ=255)                                                                                                                                                                                                                                             gCCCGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACgt       <  1:387403/250‑1 (MQ=255)                                                                                                                                                                                                                                                cGTTTTT‑TTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACgtttgg   >  2:533757/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gttttttttTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACgtttggg  >  2:580144/2‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ATTTTGTTGATGCGATGGCGGGTCTGGTTACGGCTGTATTCATCGTGCTGACCTGTAACATCGTAACAGGCATCATGATCGGCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTT‑CACCATACTGATAATGACCCTCGTGGTCTCGCTGTCCGGGGTGGTCACTCGTGTCATGCCCTTTCAGATCCCGCTTCCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGT  >  NC_000913/4279629‑4280120

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