Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,301,071	C→T	A670T (GCA→ACA)	mdtO ←	putative multidrug efflux pump subunit MdtO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,301,071	0	C	T	98.5%	188.9 / NA	68	A670T (GCA→ACA)	mdtO	putative multidrug efflux pump subunit MdtO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/46);  minor base A (1/0);  total (22/46)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

AGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTGAAG  >  NC_000913/4300825‑4301317                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         aGCTGCGGGTCATTGCGTTGTTTACACCACTGCGCCTGCGGCCCGCCCGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCggg                                                                                                                                                                                                                                                      >  2:237466/1‑251 (MQ=255) aGCTGCGGGACATTGAGATGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCACGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGATGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGGCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCggg                                                                                                                                                                                                                                                      >  2:527284/1‑251 (MQ=255)      cGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGt                                                                                                                                                                                                                                                 <  2:449018/251‑1 (MQ=255)        ggTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTcc                                                                                                                                                                                                                                               >  2:158679/1‑251 (MQ=255)                aGTTGTTTCCACCACTGCGACTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGAGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGGTGTCCCGGCTGGc                                                                                                                                                                                                                                       >  2:521968/1‑251 (MQ=255)                       gccgccacTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAAGCAGCTTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTTACGATTGAGCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGCCGTTAATTTCGGGGCTGTCCAGTCTGTCAAGCggg                                                                                                                                                                                                                                <  2:360659/247‑1 (MQ=255)                            cacTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGaa                                                                                                                                                                                                                            >  1:547165/1‑250 (MQ=255)                                  gcCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGGCGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCTGTCTGGAACGCGGGCGGAAAGCTgg                                                                                                                                                                                                                     >  2:349381/1‑251 (MQ=255)                                  gcCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTgg                                                                                                                                                                                                                     <  1:181965/251‑1 (MQ=255)                                   ccTGCGGCCAGCCGGAGCGGGCAAGATGAATATCGTCGGCCAGTTTGATGTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGAGCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGc                                                                                                                                                                                                                    <  2:176741/251‑1 (MQ=255)                                      tcGGCCAGTCGGAGCTGGCAAGATGAATATCGTCGGCCAGTGTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCGTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTTCTGCCGAGAATGCTGCACAACAGCATACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAg                                                                                                                                                                                                                  <  1:102221/249‑1 (MQ=255)                                        ggCCAGCCGGAGCTGGCAAGATGAATATCGATGGCCAGTTTGCTTTGTTCCGGTTTGAGCTGTTGATCAGGTGCAGATTCCTTACGTACTTGGGCCCAGCCGGAAATCAGCGTAGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGAGGTTCCTTGTGCCTGTTACGTGGCGGGAGTTAATGTCGGGGCTGTCCAGCCCGGCAAGCGGGCGGAAAGATGGCAggcg                                                                                                                                                                                                               >  1:255149/1‑248 (MQ=255)                                            aGCCGGAGCTGTCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTGAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGtt                                                                                                                                                                                                           <  2:577173/251‑1 (MQ=255)                                                ggAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGATGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTc                                                                                                                                                                                                       >  1:486941/1‑251 (MQ=255)                                                 gAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTGTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCATCGTCGTGCGGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAATCGGGCGAAAAGCTGGCAGACTTGTTGCTcc                                                                                                                                                                                                      <  2:243906/251‑1 (MQ=255)                                                  aGCTGGCAAGATGAATATCGTCGTCCAGTTTGATTTGTTCCGGTGTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGGTGCCGAGAATGCTGCGCAACAGCAGACGTGAAAGTTTACGATTGATCATTGCGTGGCTCCTTGTGCCGTTTCCGTGGCGGGCGTGAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAATCGGGCAGACGTGTTGCACCt                                                                                                                                                                                                     <  2:228265/251‑1 (MQ=255)                                                   gCTGGCAAGATTAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACTTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCCCCGTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGTTGGCAGACGGGTTGCTCCT‑g                                                                                                                                                                                                   <  2:55667/251‑1 (MQ=255)                                                    cTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCCGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCCGACGTGTTGCTCCT‑Gt                                                                                                                                                                                                  >  1:263275/1‑251 (MQ=255)                                                     tGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGGCCAATCTGGCAAGCGGGGGAAAAGCTGGCCGACGTGTTGCTCCT‑Gtt                                                                                                                                                                                                 >  2:280567/1‑251 (MQ=255)                                                             tGAATATCGTCGGCCAGTTTGATTTGTTCCGGATTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTCATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGCAAAGCTGGCAGACGGGTTGGTCCT‑GGTTttata                                                                                                                                                                                           >  2:422973/1‑249 (MQ=255)                                                                  atgttcgGCCAGTTTGATCTGTGCCGTTTGGAGCCGCTTATGACGTGCAGAATCCTTACGGACCAGGGCACAGCCGGAAATCAGCGTCGTGCTTCCGAGAATGCTGCACAGCAGGAGACGTGAAAGTTGACGATTGATCATTGCTTAGCTCATTGTGACTGTTCCGTGTCCGGCGTTAATGTGGTGGGTGTCCAGACTGACAAGCGGGCGATTAGCTGGTAGACGTGTTGCGGCT‑GGTTTAAAAAGGtcgg                                                                                                                                                                                    <  2:56142/247‑4 (MQ=255)                                                                                 gATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTAAGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACCGCAGCCGTGAACGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGACGGGGCTGTCCACTATGGCCAGCGGGCGGAAAGCTGGCAGACGTGTTGGACCTAGTTTCAA‑AAGGTGGGCGGGGTTGCCTgagg                                                                                                                                                                     >  2:234547/1‑250 (MQ=255)                                                                                    ttGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAgg                                                                                                                                                                   <  1:254583/250‑1 (MQ=255)                                                                                     tGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGc                                                                                                                                                                  <  1:119696/250‑1 (MQ=255)                                                                                         ccGGTTTGAGCTGTTGATGAGGTGCAGAATCCGTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTTAAAGTTGACGATTGATCAGTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGgtg                                                                                                                                                             <  2:483796/251‑1 (MQ=255)                                                                                               tGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAAGTTCGGGGCTGTCCAGTCTGTGAAGCGGGCGGAAAGCTGGCAGGCGGGTTGCTCCT‑GGTTTTATAAGGGGGGGAGGATGGCCTGAGAGGCTGgtgtgtact                                                                                                                                                       >  2:131003/1‑246 (MQ=255)                                                                                                 aGCTGTTGATGAGGTGCATAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGTAGACGTGAAAGTTGACGTTTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGTGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCa                                                                                                                                                     <  2:478672/251‑1 (MQ=255)                                                                                                 aGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCggggggttcttc                                                                                                                                                      >  1:125687/1‑241 (MQ=255)                                                                                                  gCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCaa                                                                                                                                                    <  1:508477/251‑1 (MQ=255)                                                                                                   cTGTTGATGCGGTGCAGCATCCTTACGTACCATGGCACAGCCGGAAATCGGCGTCGAGCTGCCGAGAATGCTGCACAACAGCAGACGTGAACGTTGACGATTGCTCATGGCGTGGCTCCTAGTGCCTGTTCCGTTGCGGGCGTTAATGTTGGGGCTGGCCAGGTTGGATAGCGGGCGGAAAGGCGGCCAACGTGTTTCACAT‑GTTTTTATAAGGGGGGGAGGCTGGCCTgagttgcagggggttttttcg                                                                                                                                                    >  2:336761/1‑232 (MQ=255)                                                                                                     gTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAg                                                                                                                                                   <  2:59652/249‑1 (MQ=255)                                                                                                        gatgaGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGtt                                                                                                                                               <  1:144394/250‑1 (MQ=255)                                                                                                                   gAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCttttagcggttttttagggga                                                                                                                                    >  1:221476/1‑231 (MQ=255)                                                                                                                   gAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCGGAGAGGCGGGTGTCTCTTCAAGCGTTAGTTGAGTAGaa                                                                                                                                   <  2:402205/251‑1 (MQ=255)                                                                                                                           aCGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCTTTAAGCGGTTTTTTAGGGGAAAGGCTggg                                                                                                                           >  1:253910/1‑249 (MQ=255)                                                                                                                             gTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAACGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTCATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCTTTAAGGGGTATTTTAGGAGGACgggttagtg                                                                                                                         >  1:578196/1‑244 (MQ=255)                                                                                                                                  aGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑tt                                                                                                                   <  1:551116/251‑1 (MQ=255)                                                                                                                                    ggCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGc                                                                                                                 <  1:550147/251‑1 (MQ=255)                                                                                                                                            cGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGTTCTTCAAAAGTTATTTTAGGAGaaaagggcggggggggtgccaggacgg                                                                                                          >  1:373040/1‑226 (MQ=255)                                                                                                                                                     gCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGGTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGTAGCGTAtt                                                                                                <  1:149364/251‑1 (MQ=255)                                                                                                                                                       gtcgtGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTtc                                                                                              <  1:38110/251‑1 (MQ=255)                                                                                                                                                         cgtGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTGGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTtctc                                                                                            <  2:120354/251‑1 (MQ=255)                                                                                                                                                               gcCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAggg                                                                                       <  1:8332/250‑1 (MQ=255)                                                                                                                                                                  gagaATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTc                                                                                   <  2:253910/251‑1 (MQ=255)                                                                                                                                                                   agaATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCCAGCGGTATTTTTGGGGAACGGCTGGGGGGGG‑GTGGCCGACCGGCAGCGTATTTTTcccaggcggcg                                                                                  >  1:469888/1‑242 (MQ=255)                                                                                                                                                                    gaATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCgg                                                                                 <  2:469888/251‑1 (MQ=255)                                                                                                                                                                     aaTGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGc                                                                                <  2:569766/251‑1 (MQ=255)                                                                                                                                                                      aTGCTGCACAACAGCCGGCGTGAAAGGTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCa                                                                               <  2:623/251‑1 (MQ=255)                                                                                                                                                                      aTGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTTTCTTTCAGCGTTTTTTTTGGGGAACGGCTGGGGGCGGGTTGCC‑GGCCGGCAGCGGATTTTTCCCGGGGGTCGGcc                                                                               >  1:206568/1‑250 (MQ=255)                                                                                                                                                                       tgctgcACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCg                                                                                     <  1:55544/244‑1 (MQ=255)                                                                                                                                                                       tgctgcACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTTTTTTTAAACGGTATTTTGGGGGAAAGGGTGGGGGGGG‑TTTCCAGACCGGCAGCGGATTTTTCCCGGGCggccgcc                                                                               >  1:194420/1‑244 (MQ=255)                                                                                                                                                                                       aCGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTGGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTTTCGCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTc                                                              <  2:366517/251‑1 (MQ=255)                                                                                                                                                                                             aaGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAAGAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGAc                                                                                                            <  1:48534/199‑1 (MQ=255)                                                                                                                                                                                              aGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGt                                                       <  1:99671/251‑1 (MQ=255)                                                                                                                                                                                                   aCGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGc                                                  <  1:294454/251‑1 (MQ=255)                                                                                                                                                                                                    cGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGTGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCTG‑TTTCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCg                                                 <  1:379228/251‑1 (MQ=255)                                                                                                                                                                                                            cATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCa                                          <  2:206364/250‑1 (MQ=255)                                                                                                                                                                                                                 cGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCt                                    <  1:280567/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGt                    <  1:349381/250‑1 (MQ=255)                                                                                                                                                                                                                                  gTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCGCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTg                   <  1:221875/251‑1 (MQ=255)                                                                                                                                                                                                                                  gTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGGCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGt                    <  2:547165/250‑1 (MQ=255)                                                                                                                                                                                                                                        gggcgggcgTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATcc             <  2:373732/250‑1 (MQ=255)                                                                                                                                                                                                                                        gggcgggcgTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATcc             <  1:10467/250‑1 (MQ=255)                                                                                                                                                                                                                                          gcgggcgTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGTGCCTGCGCCGAGTTCCAGGTGGCATCCCa           <  1:422973/251‑1 (MQ=255)                                                                                                                                                                                                                                            gggcgTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAgg         <  1:53469/251‑1 (MQ=255)                                                                                                                                                                                                                                             ggcgTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGc        <  1:237466/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aaTGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTgaag  <  1:521968/251‑1 (MQ=255)                                                                                                                                                                                                                                                   aaTGTCGGGGCCGTCCAGTCTGGCCAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTGAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGATTGCGG‑TTGCCAGACCGGCACCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCACTGGCCGGCTGTAGTGCGTTAACCAGCGCCTGCGCCGAGTTCCATGTGTCATCCCAGTCGTgaag  <  2:286533/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         AGCTGCGGGTCATTGAGTTGTTTCCACCACTGCGCCTGCGGCCAGCCGGAGCTGGCAAGATGAATATCGTCGGCCAGTTTGATTTGTTCCGGTTTGAGCTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCT‑GTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGG‑TTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTGAAG  >  NC_000913/4300825‑4301317

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