Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,794,130	A→G	A48A (GCA→GCG)	rfaD →	ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,794,130	0	A	G	100.0%	154.1 / NA	46	A48A (GCA→GCG)	rfaD	ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/17);  total (29/17)

TTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGCCTGGAGCGTGAAATCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTG  >  NC_000913/3793883‑3794371                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   aatctcacttttggcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGCAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCCTCACCGATATTCTGGTGGGGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGAc                                                                                                                                                                                                                                                 >  1:220194‑M2/57‑251 (MQ=255)              gcgaatggctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAgg                                                                                                                                                                                                                                    <  1:355643‑M2/208‑1 (MQ=255)                     gctgaagccaacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTTTATGGATAAGGAAGACtt                                                                                                                                                                                                                             >  1:436916‑M2/37‑251 (MQ=255)                         aagccaacgccgtgaagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCctt                                                                                                                                                                                                                         >  1:138367‑M2/33‑250 (MQ=255)                             caacgccgtgtagtgatatgatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATcc                                                                                                                                                                                                                     >  1:238543‑M2/29‑251 (MQ=255)                                                gatagaagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAaga                                                                                                                                                                                                  <  1:806‑M2/242‑1 (MQ=255)                                                     aagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAtgg                                                                                                                                                                                                            <  1:304495‑M2/232‑1 (MQ=255)                                                     aagtATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCAACAAGGTTGTGAACCTGGGGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTtcctgatccagattatggggatcggaagggggt                                                                                                                                                                                             >  2:304495‑M2/5‑236 (MQ=255)                                                       gtATTCGCGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTCAGGCCCTGAATGCTAAAGGCATCACCGGTATTCTGGTGGGTGACAACCCGAAAGGCGGCACCAAGTTTGTGACCCTGGTGGGTCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCCGAttatggctggcggaaggtttggg                                                                                                                                                                                           >  2:171138‑M2/3‑240 (MQ=255)                                                                   gagaTTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGGCAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGTTCGGCGCTGTCgcaacg                                                                                                                                                                               >  2:289460/1‑246 (MQ=255)                                                                      aTTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGTTCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGTTCGGCGAAGTCGAAGCggtt                                                                                                                                                                            >  2:4637/1‑248 (MQ=255)                                                                      aTTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAtt                                                                                                                                                                            <  2:312104/251‑1 (MQ=255)                                                                                 cTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGGCTATATTGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAgg                                                                                                                                                                >  1:34230/1‑251 (MQ=255)                                                                                        aTTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAAGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGc                                                                                                                                                         <  2:73673/251‑1 (MQ=255)                                                                                             aaGGTTACAGTTATGATCATCGTTACCGGCGGCGTGGGCTTTATCGGTCGCAACATCGTTAAGGCCCTGGATGATAAAGGCATCACCGATATTCTGGTGGGGGACAAACTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGAATATATGGATAAGGAAGACTTCATGATCCAGATTATGGCTGGCGAAAGGTTTGGAGATGTAGGAACGGTTTT‑CCACGAAgtgggatgctcttt                                                                                                                                                    >  2:20724/1‑238 (MQ=255)                                                                                                gTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATCAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGACGCGATTTT‑CCACGAAGGCGCATGCTCTTcccc                                                                                                                                                 >  1:280302/1‑249 (MQ=255)                                                                                                     aGTTAGGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACAGCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCg                                                                                                                                            <  2:166820/251‑1 (MQ=255)                                                                                                       ttATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGACCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATTGCGAAAGGATTTTTCCA‑GAAGGCGCATGCTCTTccaaccacgcg                                                                                                                                          >  2:353819/1‑243 (MQ=255)                                                                                                                  gTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGGGCATGCTCTTCCACCACCGAGTGGGATGGCa                                                                                                                                >  1:205347/1‑250 (MQ=255)                                                                                                                     aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTTAATGATAACGGCCTCCCCGATATTCTGGGGGTGGAACACCTGAAAGAAGGCACCAAGTTTGTGAAACTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGGGATGTCGGAGAGCTTTT‑TCACCAAAGCGCCTGCCCTTCCACCCCCgcgggggatggccagt                                                                                                                             >  2:256821/1‑235 (MQ=255)                                                                                                                     aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAATACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGAGGGCAAGta                                                                                                                            <  2:248359/251‑1 (MQ=255)                                                                                                                     aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt                                                                                                                             >  1:103130/1‑250 (MQ=255)                                                                                                                                gggCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGAt                                                                                                                  >  1:14095/1‑250 (MQ=255)                                                                                                                                    tttATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATaaca                                                                                                              >  1:313521/1‑250 (MQ=255)                                                                                                                                          ggCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCGCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCa                                                                                                       <  2:153269/251‑1 (MQ=255)                                                                                                                                                aaCATCGTTAAGGCACTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAATGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTc                                                                                                 <  1:345626/251‑1 (MQ=255)                                                                                                                                                       ttAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGCGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCaaag                                                                                            >  1:14012/1‑247 (MQ=255)                                                                                                                                                                   aTGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCactac                                                                              >  2:301444/1‑251 (MQ=255)                                                                                                                                                                     gATAAAGGCATCTCCGATATTCTGGTGGTGAACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATAGGGATAAGGAAGACTGCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCaaatagcttccgcaatact                                                                             <  2:385317/250‑18 (MQ=255)                                                                                                                                                                            gCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAg                                                                     <  1:206229/251‑1 (MQ=255)                                                                                                                                                                               tCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCACGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCgc                                                                  >  1:234911/1‑250 (MQ=255)                                                                                                                                                                                   cGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGaa                                                               >  2:209972/1‑250 (MQ=255)                                                                                                                                                                                    gATATTCTGGTGGTGGACAACCTGAAAGACGGCAACAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAAt                                                             <  2:436916/251‑1 (MQ=255)                                                                                                                                                                                                  ggACAACCACAAAGCCAGCACCACCTTTGTGACCCTGGTGGTTCTCGATCTCGCGGACTATATCGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATTTCGAAGCGATTTT‑CCACGAAGCCGCATGCTCTTCCACCACCGAGAGGGATGTAAAGTATATGATGGGAAACACCTATCAATACTCCAAACACCTGCTGCAATACTGTCTGGAGCACGAAATCCCATTCCTGTa                                                 >  2:146211/1‑249 (MQ=255)                                                                                                                                                                                                               gACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGcggc                                  <  1:171138/251‑4 (MQ=255)                                                                                                                                                                                                                  ggCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCAc                               >  2:330023/1‑251 (MQ=255)                                                                                                                                                                                                                    cACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTG‑CCACGAAGGCGCATGCTCTTCCACCAGCGAGGGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCGGTTCCTGTATGCAGCGTCCGCGGCCACCt                             <  1:146211/251‑1 (MQ=255)                                                                                                                                                                                                                        aaGTTTGTGAACCTGGTGGATATGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CAACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGAGCGAAATCCCGTTCATGTAGGCATCTTACGCGGCCACCTAcgg                         >  1:165330/1‑251 (MQ=255)                                                                                                                                                                                                                             tgtgAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGAc                    >  1:256254/1‑251 (MQ=255)                                                                                                                                                                                                                                   ccTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCAccc              >  1:71251/1‑250 (MQ=255)                                                                                                                                                                                                                                      ggtggATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCg           >  1:108466/1‑251 (MQ=255)                                                                                                                                                                                                                                               ggatatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTAt    <  2:103130/247‑1 (MQ=255)                                                                                                                                                                                                                                               ggatatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGGCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCGCCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTAGCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAACCCCGTTCCTGGATGCATCTTACGCGGCCACCTACGGCGGACGCATCTCCGCCtt       <  2:385571/244‑1 (MQ=255)                                                                                                                                                                                                                                               ggatatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAAGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATTGCAAGTATATGATGGATAACAACTATCATTACTCCAAAGAGCTATTGCACTACTGTCTGGAGCGCGAAATACCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCAGACGCACCTCCGACTTTATTg  >  1:81556/3‑251 (MQ=255)                                                                                                                                                                                                                                                  tatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGCCGCACCCCCGGCttt      >  2:207623/1‑244 (MQ=255)                                                                                                                                                                                                                                                  tatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACttt      <  1:207623/244‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   TTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGCCTGGAGCGTGAAATCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTG  >  NC_000913/3793883‑3794371

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