Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,794,122	A→G	N46D (AAT→GAT)	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,122	0	A	G	100.0%	141.0 / NA	41	N46D (AAT→GAT)	rfaD	ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/14);  total (27/14)

TTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGCCTGGAGCGTGAAATCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCG  >  NC_000913/3793883‑3794362                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  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)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  TTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGCCTGGAGCGTGAAATCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCG  >  NC_000913/3793883‑3794362

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: ‑