Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,026,473	Δ1 bp	coding (1486/1494 nt)	ubiD →	3‑octaprenyl‑4‑hydroxybenzoate decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,026,473	0	A	.	100.0%	297.4 / NA	63	coding (1486/1494 nt)	ubiD	3‑octaprenyl‑4‑hydroxybenzoate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (27/36);  total (27/36)

TTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTG  >  NC_000913/4026229‑4026722                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          ttGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCt                                                                                                                                                                                                                                                     <  1:589298/250‑1 (MQ=255)  tGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTg                                                                                                                                                                                                                                                    <  2:678277/250‑1 (MQ=255)              accCGTATGGACCCGGCGCATGATACTGTTCTGGTAGAAAATACGCCTATTAATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCgtttgtt                                                                                                                                                                                                                                       <  2:682978/251‑1 (MQ=255)                  gTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGc                                                                                                                                                                                                                                   >  1:523478/1‑251 (MQ=255)                       gACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTAt                                                                                                                                                                                                                              >  2:609141/1‑251 (MQ=255)                         cccGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTAttt                                                                                                                                                                                                                            <  1:309916/251‑1 (MQ=255)                                 gggATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATcc                                                                                                                                                                                                                    <  1:92181/251‑1 (MQ=255)                                  ggATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCg                                                                                                                                                                                                                   >  2:6011/1‑251 (MQ=255)                                         gTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACag                                                                                                                                                                                                               >  1:209169/1‑248 (MQ=255)                                              ggTAGAAAATCCGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAgc                                                                                                                                                                                                        <  2:403289/250‑1 (MQ=255)                                              ggTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAgcg                                                                                                                                                                                                       >  1:98838/1‑251 (MQ=255)                                              ggTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAgc                                                                                                                                                                                                        <  2:43690/250‑1 (MQ=255)                                                  gAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATg                                                                                                                                                                                                   >  1:60365/1‑251 (MQ=255)                                                          gCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCtt                                                                                                                                                                                           <  2:588539/251‑1 (MQ=255)                                                             tATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATCCCACGAATAAATGGCCGGGGGAAACCCCGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCAATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGAGCATGACAACCTTa                                                                                                                                                                                          >  1:286107/1‑249 (MQ=255)                                                                tGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCt                                                                                                                                                                                      <  2:212673/250‑1 (MQ=255)                                                                       cTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTg                                                                                                                                                                              >  2:335369/1‑251 (MQ=255)                                                                              ttGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCg                                                                                                                                                                        <  1:567831/250‑1 (MQ=255)                                                                               tGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGt                                                                                                                                                                      >  1:271576/1‑251 (MQ=255)                                                                                           ctcCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCAc                                                                                                                                                          >  2:566779/1‑251 (MQ=255)                                                                                             ccGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCa                                                                                                                                                           <  2:514054/248‑1 (MQ=255)                                                                                               gggCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGAt                                                                                                                                                      <  2:293257/251‑1 (MQ=255)                                                                                                     ggTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGt                                                                                                                                                 >  1:303631/1‑250 (MQ=255)                                                                                                      gTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGtat                                                                                                                                               >  1:576226/1‑251 (MQ=255)                                                                                                       gagttgtatatgtccaagtctgccGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCGGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGtata                                                                                                                                              <  2:165682/228‑1 (MQ=255)                                                                                                          aaaaaTGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCg                                                                                                                                           >  1:633278/1‑251 (MQ=255)                                                                                                             aaTGGGGCTGGATGCCACGAATAAATGGCCGGGGGAATCCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCgtgt                                                                                                                                        <  1:581814/251‑1 (MQ=255)                                                                                                             aaTGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGTCCTCGGTAGAAGCTATCACGGATACCGTATATCgtgt                                                                                                                                        <  2:91685/251‑1 (MQ=255)                                                                                                             aaTGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCgtgt                                                                                                                                        <  2:206119/251‑1 (MQ=255)                                                                                                                  ggCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTcc                                                                                                                                      >  1:587947/1‑248 (MQ=255)                                                                                                                  ggCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTcc                                                                                                                                      <  2:587947/248‑1 (MQ=255)                                                                                                                      ggATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATc                                                                                                                                 <  1:335369/249‑1 (MQ=255)                                                                                                                      ggATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATc                                                                                                                                 <  2:209169/249‑1 (MQ=255)                                                                                                                        aTGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGc                                                                                                                             <  1:124553/251‑1 (MQ=255)                                                                                                                                       ggCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATc                                                                                                                                 <  1:552412/232‑1 (MQ=255)                                                                                                                                       ggCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATc                                                                                                                                 >  2:552412/1‑232 (MQ=255)                                                                                                                                           gggggAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCcttttcttt                                                                                                          <  2:136319/251‑1 (MQ=255)                                                                                                                                             gggAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCttttcttttc                                                                                                        <  2:531446/251‑1 (MQ=255)                                                                                                                                              ggAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCg                                                                                                       >  1:71908/1‑251 (MQ=255)                                                                                                                                               gAAACCCAGCGCGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGt                                                                                                      >  2:379114/1‑251 (MQ=255)                                                                                                                                                aaaCCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTg                                                                                                     <  2:60365/251‑1 (MQ=255)                                                                                                                                                     cAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGt                                                                                                >  1:145287/1‑251 (MQ=255)                                                                                                                                                         gTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGt                                                                                            >  2:218663/1‑251 (MQ=255)                                                                                                                                                           gAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTAGTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCTGTAt                                                                                          <  2:271192/251‑1 (MQ=255)                                                                                                                                                                   aCGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGGATTTGATggg                                                                                  >  1:127670/1‑250 (MQ=255)                                                                                                                                                                     gTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTGTCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGt                                                                                  >  2:806/1‑249 (MQ=255)                                                                                                                                                                     gTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTa                                                                                 >  2:151437/1‑250 (MQ=255)                                                                                                                                                                       cccATCAAAAAAGATCCAGATGTTGTCGCTCATATTGACGCCATCTGGGATTAACTGGCTATTTTTAACAACTGTAAA‑GCGGCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATTACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTg                                                                              <  2:624413/251‑1 (MQ=255)                                                                                                                                                                              aaaaaGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTgatggatn                                                                       <  2:303631/251‑2 (MQ=255)                                                                                                                                                                                aaaGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGGCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGTTAGTGATGGATGAg                                                                     <  1:505336/251‑1 (MQ=255)                                                                                                                                                                                          aTGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTGTCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAAc                                                           <  1:806/251‑1 (MQ=255)                                                                                                                                                                                          aTGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAAc                                                           >  1:458034/1‑251 (MQ=255)                                                                                                                                                                                           tgttgtCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACg                                                          >  2:483225/1‑251 (MQ=255)                                                                                                                                                                                                 cgcgCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGGCGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGtt                                                    <  2:572827/251‑1 (MQ=255)                                                                                                                                                                                                        atTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAAt                                              <  2:127670/250‑1 (MQ=255)                                                                                                                                                                                                        atTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGCCCGTTCTCAATg                                             <  1:483225/251‑1 (MQ=255)                                                                                                                                                                                                           gACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCt                                          <  1:218663/251‑1 (MQ=255)                                                                                                                                                                                                                     gggATGAACTGGCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCg                                 <  2:271576/250‑1 (MQ=255)                                                                                                                                                                                                                               ggCTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGaaaaa                        <  2:576226/249‑1 (MQ=255)                                                                                                                                                                                                                                 cTATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGt                    >  1:10768/1‑251 (MQ=255)                                                                                                                                                                                                                                  tATTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCACTGGCTTCGACGCCGGATGAAAAAGGGtt                   >  2:635213/1‑251 (MQ=255)                                                                                                                                                                                                                                   aTTTTTAACAACGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGttt                  <  2:145287/251‑1 (MQ=255)                                                                                                                                                                                                                                            aacGGTAAA‑GCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTg         >  2:188249/1‑251 (MQ=255)                                                                                                                                                                                                                                                    aaGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAACCGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTg  >  2:166948/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          TTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTG  >  NC_000913/4026229‑4026722

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