Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	826,038	G→A	A24V (GCG→GTG)	ybhR ←	ABC exporter membrane subunit YbhR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	826,038	0	G	A	97.4%	118.3 / NA	39	A24V (GCG→GTG)	ybhR	ABC exporter membrane subunit YbhR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/24);  minor base C (0/1);  total (13/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCAGCCGGGAAACGCACCAGTAGTAACGCCTTTTGTGTGTCGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTAC  >  NC_000913/825796‑826278                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  aGTCAGCCGGGAAACGCACCAGTAGTAACGCCTTTTGTGTGTCGATGGTTGGGCGGATCTCCGGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGtt                                                                                                                                                                                                                                           <  2:377768/251‑1 (MQ=255)    cAGCCGGGAAACGCACCAGTAGTAACGCCTTTTGTGTGTCGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTg                                                                                                                                                                                                                                         <  1:27746/250‑1 (MQ=255)       ccGGGAAACGCACCAGTAGTAACGCCTTTTGTGTGTCGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGt                                                                                                                                                                                                                                     >  1:13865/1‑251 (MQ=255)         gggAAACGCACCAGTAGTAACGCCTTTTGTGTGTCGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTc                                                                                                                                                                                                                                   >  1:655/1‑251 (MQ=255)            aaaCGCACCAGTAGTAACGGCTTTTGTGTGTCGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGTGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTAACTTACAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGATTGCGGTTc                                                                                                                                                                                                                                   >  1:320358/1‑248 (MQ=255)            aaaCGCACCAGTAGTAACGCCTTTTGTGTGTCGATGGTGTGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAACCGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTc                                                                                                                                                                                                                                   <  2:320358/248‑1 (MQ=255)                         gtaACGCCTTTTGTGTGTCGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGc                                                                                                                                                                                                                   >  2:15647/1‑251 (MQ=255)                         gtaACGCCTTTTGTGTGTCGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGc                                                                                                                                                                                                                   <  1:165943/251‑1 (MQ=255)                                           cGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTGCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGTGGGTTTGCGGTTCGCGCAGCAACGACGGCAACTCTTTGCGGATTAAc                                                                                                                                                                                                 <  1:334456/251‑1 (MQ=255)                                                             ccTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAAc                                                                                                                                                                               <  1:410129/251‑1 (MQ=255)                                                             ccTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGACGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATg                                                                                                                                                                                   <  1:283808/247‑1 (MQ=255)                                                                gtgGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGGCCATAAGCGATGAAACAt                                                                                                                                                                             <  2:340223/250‑1 (MQ=255)                                                                                cagnccaTGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCTGCTCCACCGAATGTTCGCCGTTATCTTTATCGTAGATGGCGATTGTGGCGTTAGTCACTTCTAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCTCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATGAACGTCCATAAGCGATGAAACATGCTCGTCTCCCTAATcc                                                                                                                                                            <  2:404275/246‑1 (MQ=255)                                                                                    acaTGAGTAAAGGCGCGGGCGCGTGCAAAACGTTGGTTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACTGGTAAAATCAGAATCACGCGGGTTTGCGGGTCGCGCAGCACCGACTGCAGCTCTTTGCGGATTAACGTCCAGAAGCGATGAAACAGGCTCTGCTCCCTAATCCAGcc                                                                                                                                                        <  2:281065/251‑1 (MQ=255)                                                                                             aaGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGAGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGCAGTTTgg                                                                                                                                               >  2:20431/1‑251 (MQ=255)                                                                                                   cTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTc                                                                                                                                          >  2:10363/1‑250 (MQ=255)                                                                                                     ggCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑Gc                                                                                                                                      <  2:116928/251‑1 (MQ=255)                                                                                                                aaCGTTGGGTCAGCTCCTCCGAAGGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGcc                                                                                                                           <  2:187308/251‑1 (MQ=255)                                                                                                                  cGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGa                                                                                                                         <  1:450459/251‑1 (MQ=255)                                                                                                                   gTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGAt                                                                                                                        <  1:175743/251‑1 (MQ=255)                                                                                                                      gggTCAGCTCCACAGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGTCCGATaaa                                                                                                                     <  2:160124/251‑1 (MQ=255)                                                                                                                         tCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACAt                                                                                                                  >  1:144875/1‑251 (MQ=255)                                                                                                                                cACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCAccgcc                                                                                                           <  1:256371/251‑1 (MQ=255)                                                                                                                                   cGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCa                                                                                                                                         <  1:243592/219‑1 (MQ=255)                                                                                                                                   cGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCa                                                                                                                                         >  2:243592/1‑219 (MQ=255)                                                                                                                                        gCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGAt                                                                                                   <  1:15647/251‑1 (MQ=255)                                                                                                                                             ccGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGTTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCAAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCgaaa                                                                                              <  1:229041/251‑1 (MQ=255)                                                                                                                                                      tcatcggagagtgCGATGGTGGGGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGAAGCAACGACTGCAACGTTTTGCGTATTAACGTCCAGAAGCGAGAAAACATGCTCTTCTCCCTAATCCGGCCGTCGTTTGGTTTTCA‑GCCCCGTCAGCCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCAc                                                                                       <  1:58056/237‑1 (MQ=255)                                                                                                                                                                    gATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACATGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcacc                                                                       >  1:86265/1‑251 (MQ=255)                                                                                                                                                                        gtggCGTTAGTCACTTCCAGCGGGGAGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAAACACGCGGGTTTGCGGTTCGCGCAGCAACGACGGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTg                                                                   <  2:348064/251‑1 (MQ=255)                                                                                                                                                                            cGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAat                                                               >  1:132005/1‑251 (MQ=255)                                                                                                                                                                                             gTGGCGGCGACCGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCCCGCGGTTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATTCTCTTCTCCCTAATCCAGCCGACTTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACa                                              <  2:383482/251‑1 (MQ=255)                                                                                                                                                                                             gTGGCGGCGAACGGGAACAGGATCACCGGAATTAGCACGGGTAAAATCAGATCCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGAGGATCAACGTCCATAAGCGATGAAACATGCTCTTCCCCCTAATCCAGCCGACGTTTG‑TTTTCAAGCCCCGTCAGTCCGATAAACATCCCCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACa                                              <  2:132005/251‑1 (MQ=255)                                                                                                                                                                                                           gAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGTCTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGAACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCCGGGTg                                >  1:289253/1‑251 (MQ=255)                                                                                                                                                                                                                    cACCTTAATTAGCACGGGTAAAATCAAAATCACGCGGTTTTGCGGTTCGCGCAGCAACCACTGCAACTCTTTTCTAATTAACGTCCATAAGCAATGAAACATGCTCTTCTCCCTAATCCAGCCGACTTTTGGTTTTCA‑GCCACGTCAGGCTGATAAACTTCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCTGCGAGGAACAGGCTTTGCAGGGTGCTGACGaa                        <  1:178652/250‑1 (MQ=255)                                                                                                                                                                                                                     aCCTGAATTAGCACGGGTAAAATCAGAATCACGCGGGATTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCAGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGNAAt                      >  2:413448/1‑251 (MQ=255)                                                                                                                                                                                                                              aGCACGGGTAAAATCAGAATCACGGGGGTTTGCGGTTCGCGCAGCAACGCCGTCAACTCTTGGCGGATTAACGTCCAAAGGCGATGAAACATGCTCTTCTCCCTATTCCAGCCGACGTTTGGTTTTCA‑GCCAGGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCgg             <  2:163772/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaaTCAGAATCACGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGt    >  2:425275/1‑250 (MQ=255)                                                                                                                                                                                                                                         aaTCAGAATCACGCGGGTTGGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTAc  <  1:10363/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  AGTCAGCCGGGAAACGCACCAGTAGTAACGCCTTTTGTGTGTCGATGGTTGGGCGGATCTCCTGTGGGCTTTTCAGCAGCAGCACATGAGTAAAGGCGCTGGCGCGGGCAAAACGTTGGGTCAGCTCCACCGAATGCTCGCCGTTATCTTCATCGTAGATGGCGATGGTGGCGTTAGTCACTTCCAGCGTGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCA‑GCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGCTGACGAAATAACGAGCGGGAATAATGTAC  >  NC_000913/825796‑826278

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