Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	502,232	G→A	T336I (ACC→ATC)	ybaL ←	putative transporter YbaL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	502,232	0	G	A	94.6%	109.9 / NA	38	T336I (ACC→ATC)	ybaL	putative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (20/15);  minor base C (2/0);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGATCTGCTTCTCTTCTTCGATTGCCTCTTCCAGCGTCTGCTCTTCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCAGTACCATCCCGGC  >  NC_000913/502001‑502470                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                gggATCTGCTTCTCTTCTTCGATTGCCTCTTCCAGCGTCTGCTCTTCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTgg                                                                                                                                                                                                                                <  1:1987/249‑1 (MQ=255)     tCTGCTTCTCTTCTACGATTGCCTCTTCCAGCGTATGCTCTTCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGttaggcgtg                                                                                                                                                                                                                           >  2:347395/1‑243 (MQ=255)        gCTTCTCTTCTTCGATTGCCTCTTCCAGCGTCTGCTCTTCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGAcc                                                                                                                                                                                                                        >  1:402808/1‑250 (MQ=255)              ctttgtcGATTGCTTCTTCCTGCGTCTGCTCTTCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGTTTGAGGATAATCGACAGGAGCGCCCCGTCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGAGAAACGCGAACTCACCAATCGGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACa                                                                                                                                                                                                                   <  1:20578/244‑1 (MQ=255)                                                    gTTTCGGTCTTCGCCAGAATTTTCCCCAGTAGTGCGAACATTCCCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAg                                                                                                                                                                          <  2:320434/251‑1 (MQ=255)                                                            cTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTAACAAACAGAATAATCg                                                                                                                                                                  >  1:366961/1‑251 (MQ=255)                                                                gccggATATATCTCCAGGAGTGCGAAAAGGACCGGGGTAAGCATAATCTACAGGATGGACCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGGCCCGCCAGGATAAACGAGAACTCACTAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGGCACCAAACAGACGCACCAGGAAAAAGGCGGCTACAC‑ACTTACCAACCAGAAAAATcacacg                                                                                                                                                              >  2:99642/5‑246 (MQ=255)                                                                gCCAGATATTTCTCCAGTATTGCGAACAGTGCCGGGTTTAGCATAATCGACAGGATCGCCCCTGCCATTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCTCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAg                                                                                                                                                              <  1:385096/251‑1 (MQ=255)                                                                          tctcCAGTAGTACGACCAGTACCGGGGTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGCCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAg                                                                                                                                                     <  2:16214/250‑1 (MQ=255)                                                                             ccAGTAGTGCGAACAGTACGGGGTTGAGCATAAGCGACAGGATCGCCCCTGCCAGTCCCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACt                                                                                                                                                 <  2:110605/251‑1 (MQ=255)                                                                                gtagtGCGACCAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGAACCAGGTTTTGTCCGGCCTGCGGCAGTCAATTCAATCCCATTCCCAGTCCCCCCAGGATAAACGCGAACTCACCAATCTGCCCCAGGCTGACGGCGATGCTTAATGCCGTACGTTGGGAGTGCCCAACCAGTCGCACCAGAAAACATGCGGCTA‑ACGACTTACCAAACCGAATAATCTCCAGCGTCGCctgccctgcc                                                                                                                                              >  1:20304/1‑251 (MQ=255)                                                                                gtagtGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGCAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCCACACTGcc                                                                                                                                              >  1:133459/1‑251 (MQ=255)                                                                                 tagtGCGAACACTACCGGGTGGAGCATAACCGACAGGATCCCCCCTGCCAGTACCGGGTTTTGGCCGGCCTGCGGCAGTAAATCCAATGCTATTCCCAGCCCCGCCAGGATAAACGCGAACTCACCAAGTTGCGCCAGGCGGGCGGCGATGATTAATGCCGTGCGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAAGGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGcc                                                                                                                                              <  1:239896/250‑1 (MQ=255)                                                                                      cGAACCGTACCAGGTTGAGCATAATAGACAGGATCGCCCCTGCCAGCACCAGGCTTTGTCCGGCCTGCGGCAGTACATTCAACGCCATTCCCAGTCCCGACAGGATAAACGCGAACTCACCAATCTGCGCCAGCATGGCGGCGCTGCTTCATGCCGTCCGATGGGAGTGACCAAACAGTCGCCCCAGGAAAAATGCGGCTA‑CCGACTTACCAAACAGAATAAACGACCGCGTCGCCAGCACTGCCAGCgg                                                                                                                                         >  2:63302/1‑250 (MQ=255)                                                                                       gAACAGGACCGGTTTGAGAATAATCGACAGGACCGCCCCTGCCATTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGtt                                                                                                                                       <  2:298801/251‑1 (MQ=255)                                                                                        aaCAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGTCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTg                                                                                                                                      >  1:11362/1‑251 (MQ=255)                                                                                              ccggggTTGGGCATAATCAACAGGATCGCCCCTGCCAGTACCAGTTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGaa                                                                                                                                 <  2:372508/247‑1 (MQ=255)                                                                                                                 aCAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTAAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAaac                                                                                                             >  1:347039/1‑251 (MQ=255)                                                                                                                  cagtatcaccccTGCCATAACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGACTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACCGTCGAACCAGGAAAAAAGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGAACaa                                                                                                               >  1:90488/9‑248 (MQ=255)                                                                                                                       tCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCCATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCAAGCGGATGATGAATCAGAATTAACGGATCAaacaacctc                                                                                                       >  1:237216/1‑248 (MQ=255)                                                                                                                                   gTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGAc                                                                                            >  1:237345/1‑250 (MQ=255)                                                                                                                                                   ggCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCaa                                                                           <  1:391966/251‑1 (MQ=255)                                                                                                                                                           gCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAAcgcgccg                                                                   >  1:392298/1‑248 (MQ=255)                                                                                                                                                              gTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCgcgc                                                                >  1:34138/1‑251 (MQ=255)                                                                                                                                                              gTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCgcgc                                                                <  1:153710/251‑1 (MQ=255)                                                                                                                                                                       aTGCCATTCCCAGGCCCGCCAGGATAAACGCGAACTCCCCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCGGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACTCGTCGCGCAATGGCAGc                                                       <  1:338417/251‑1 (MQ=255)                                                                                                                                                                        tGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCAAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCg                                                      >  1:65432/1‑251 (MQ=255)                                                                                                                                                                                ccAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTTGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATATTCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTgg                                              <  1:19612/251‑1 (MQ=255)                                                                                                                                                                                  aGTCCCGCCAGGATAAACGAGAACGCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGCCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCAATGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAGCAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTAGCGTGggc                                            <  1:160834/251‑1 (MQ=255)                                                                                                                                                                                              aTAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGACTTAACGGATCAAACAAAATCCCGACGGAGACAAAACACAGCACCGCAAACGCGTCCCGCAATGGCAGCGTATCGTGGGCGGCACGGTGAcc                                >  2:333913/1‑250 (MQ=255)                                                                                                                                                                                              aTAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACt                                <  1:175371/251‑1 (MQ=255)                                                                                                                                                                                              aTAAACGCGAACTCACCAAGCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACt                                <  2:151700/251‑1 (MQ=255)                                                                                                                                                                                                   cgcgAACTCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGtt                           >  1:226063/1‑251 (MQ=255)                                                                                                                                                                                                    gcgAAATCACCAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGAAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGGGACTCAGTTc                          >  2:20569/1‑251 (MQ=255)                                                                                                                                                                                                          tCACCAATCTGCGCCATTATGCCGGCGATGATTAATGCCGTACGCTGGGCGTGCCCCAACAGTCGCACCAGGAAAAATGCGGCTA‑AAGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGGCAAAAAACAGCTCCGCAAACGCGTCGCGCAATGCCAGCGTATCGTGGTCGGCACGGTTACCCACTTCCGACTc                    >  2:289999/1‑251 (MQ=255)                                                                                                                                                                                                             ccAATCTGCGCCAGGCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCCCCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGtt                 >  1:342620/1‑251 (MQ=255)                                                                                                                                                                                                                       ccAGTCTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAACAGCGTCGCGCCATGGCAGCGTATCGTGGGCGGCACGGCGACTCAGTTCAGACTCGTTCCGTACCAt        >  2:107215/1‑250 (MQ=255)                                                                                                                                                                                                                            cTGGCGGCGATGATTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCAGTACCATCCCGGc  >  1:49194/1‑251 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                GGGATCTGCTTCTCTTCTTCGATTGCCTCTTCCAGCGTCTGCTCTTCCAGCGTTTCGGTCTTCGCCAGATATTTCTCCAGTAGTGCGAACAGTACCGGGTTGAGCATAATCGACAGGATCGCCCCTGCCAGTACCAGGTTTTGTCCGGCCTGCGGCAGTAAATTCAATGCCATTCCCAGTCCCGCCAGGATAAACGCGAACTCACCAATCTGCGCCAGGCTGGCGGCGATGGTTAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTA‑ACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCAGTACCATCCCGGC  >  NC_000913/502001‑502470

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