Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,567,272	C→A	intergenic (‑26/+15)	yjiN ← / ← mdtM	DUF445 domain‑containing protein YjiN/multidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,567,272	0	C	A	97.5%	111.5 / ‑6.4	40	intergenic (‑26/+15)	yjiN/mdtM	DUF445 domain‑containing protein YjiN/multidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (28/11);  total (28/12)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAAC  >  NC_000913/4567026‑4567515                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ggagggCGAGAAATGACCGGAACCGTCACGCGGCGAAACAGCGCCCCCACCACAATCCAGTCCGCCAGCGCGCCGACCAGCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATGGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGCCGGCCGGAt                                                                                                                                                                                                                                                 <  2:99564/247‑1 (MQ=255)   aTGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATcg                                                                                                                                                                                                                                               >  1:86058/1‑251 (MQ=255)   aTGGCGAGAAAGGATCGGAAGCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATc                                                                                                                                                                                                                                                <  2:231734/250‑1 (MQ=255)    tGGCGAGAAATGATCGGAATAGGCACGCGGCGAAACAGCGTAACCACCGTAAACGAGTCCGCAAGCGCGCCGCCCATCGCCGATTCAGCAATCGACTTCACGCCGCTCACCCAAAAATTGGGCGGCAACAACAGCGTAACGACAAAGGCAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAGGAATAGAGGCTGAAGGCAGGAtggg                                                                                                                                                                                                                                              >  2:216446/1‑248 (MQ=255)        gaaaaaTGACCGGAACCGGCACGCGGCGAAACAGCGCCACCACCGCAACCCATGCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTAATCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCaa                                                                                                                                                                                                                                          <  2:68582/248‑1 (MQ=255)           aaaTGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCaaaa                                                                                                                                                                                                                                        >  1:58300/1‑250 (MQ=255)                 ccGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTg                                                                                                                                                                                                                                 <  2:233853/250‑1 (MQ=255)                   ggAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGct                                                                                                                                                                                                                               >  2:354952/1‑251 (MQ=255)                          ggaCGCGGCGAAACACCGCCCCCACCGCAAACCAGTCCCCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCGTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGTCGATAAGCAGTAAAGGGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCAc                                                                                                                                                                                                                        <  2:29917/249‑1 (MQ=255)                                ggcaaaaCAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTc                                                                                                                                                                                                                  <  2:7911/247‑1 (MQ=255)                                    aaaCAGCGCCACCACCGCAACCCATGCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCt                                                                                                                                                                                                              <  2:391027/251‑1 (MQ=255)                                     aaCAGCGCCACAACCGCAAACCAATCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGCTAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGCGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGCACGCGCAAAAATCCCTGCTCCTCCAATAGCtagcacg                                                                                                                                                                                                             >  2:39823/1‑245 (MQ=255)                                                 accGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGgc                                                                                                                                                                                                 >  1:96740/1‑251 (MQ=255)                                                  ccGCAAACCAGTCCGCCAGCGCGCCGACCAGCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGgcg                                                                                                                                                                                                <  1:360868/251‑1 (MQ=255)                                                     caacccAGACCGACAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCGTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACATCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTGTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGgcgc                                                                                                                                                                                               <  2:336667/245‑1 (MQ=255)                                                      aaaCAAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCCAAAATTGGGCGGCAAAAACAGCGTAACGCCTAAGGCAGCAGCGGCGATCAGCAGTAAAGAGAGCGCCACCCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTAATAACGTTAAGACTAGAGCCGGCAGGCAGCATCGCGCCAAAATCACGGCTCCTCCACTATCTCGGCTGCCTGATGCAGgcgcccg                                                                                                                                                                                            >  2:322845/1‑248 (MQ=255)                                                           aGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGAtt                                                                                                                                                                                       <  1:29242/251‑1 (MQ=255)                                                            gTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAAAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGGTCGCGCAAAAATCACTGCTCCTACCCTAGCTCGGCTGACAGATACTGGCGAAcgacattg                                                                                                                                                                                      >  2:335852/1‑245 (MQ=255)                                                                        gcCGACCATCGCCGCTTCAGCAATCGCCTTCTCGACGCTCACACAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTCGCAGCGGCGATAAGCAGGAACGAGAGCGCCAACCTTTTGGCGCGTATGAGTTCGATGAGTTTATTCATAACGTTAAGCATAGAGCCTGACGGCAGGTTCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCCATCCCgac                                                                                                                                                                          >  2:285545/1‑249 (MQ=255)                                                                         ccGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTg                                                                                                                                                                                                   <  1:43371/225‑1 (MQ=255)                                                                         ccGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTg                                                                                                                                                                                                   >  2:43371/1‑225 (MQ=255)                                                                                  cccgtgttcacatccGCCTGCACGCCGCGCACCCAAAAATTGGGCGGCAAAAACAGCGGAACGACAAAGGTTGCAGCTGCGATAATGAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCGGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGAGCAAAAATCACTGCCCCTCCACTAGCTCGGCGGCCTGATGCTGGCGCACGCGATTTAGCAATCCCGCCAGGGTGaaaa                                                                                                                                                                <  2:393869/237‑1 (MQ=255)                                                                                   ccgcTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTCTTCATAACGTTAAGAATAGAGCCTGCCGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAACCCCCCCAGGGTGAAAac                                                                                                                                                               >  2:306329/1‑251 (MQ=255)                                                                                       ttCAGCAATCGCCTACACGCCGCTCACCCACAAATTGGGCGGCAATAACAGCGTAACCACAAAGGTAGCCGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATAGCGCAAAAATCACGGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAacaacg                                                                                                                                                           >  2:207910/1‑251 (MQ=255)                                                                                                              tCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGCGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTa                                                                                                                                     >  2:234310/1‑250 (MQ=255)                                                                                                              tCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTa                                                                                                                                     >  1:83604/1‑250 (MQ=255)                                                                                                                                 cAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCgc                                                                                                                 >  1:328208/1‑250 (MQ=255)                                                                                                                                                 aCAAAGGTAGCAGCGGCGATAAGCAGTAACGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCAt                                                                                                  >  2:234352/1‑250 (MQ=255)                                                                                                                                                    acggTAGCAGCGGCGCTAAGCAGTAAAGAGAGCGCCACCCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTCAGAATAGAGCCTGACGGCAGGACCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGCCGATAACGCCCGCCACCACGCCTAACAGAAGAAACTGCCAGCGACCGCCGTTCAAcccaagcc                                                                                              >  2:271356/3‑245 (MQ=255)                                                                                                                                                                             acgagagCGCCAACCTTTTGGCGCGTCTGAGTTCTCTGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGCTTGACGAACACCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgaccgagaccgc                                                                     >  2:1872/3‑250 (MQ=255)                                                                                                                                                                               gagagCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCCAGCGACCGCCGTTAAACCATAGCCACCGGCCGATTTCGACCGAGACCGAcc                                                                   >  1:25426/1‑250 (MQ=255)                                                                                                                                                                                agagCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCGAACAGATGAAAAGGCAAGCGACCGCCGTTCAACCATAGCCAGCGGCCGATTTCAACCGAAACCGACac                                                                  >  2:9868/1‑250 (MQ=255)                                                                                                                                                                                                  cgcgTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGCTTTCGACCGAGACCGCCATCACCATCAGGATCACCa                                                 >  2:35132/1‑250 (MQ=255)                                                                                                                                                                                                                      ttATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGGGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGGAACGGCAAGCGACCGCCGTTAAACCATAACCAGCGGCCGATTTCGACCGAGACCGACAACACCATCAGGAGCACCATATTTAGCGATGCGGAGACgg                            >  2:246149/1‑251 (MQ=255)                                                                                                                                                                                                                             tAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCttt                     >  1:330807/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCATCCACTAGCTCGGCTGCCAGATGCTGGCGCAAGCGATTGAGCAATCCCGCCAGGGAGAAAACGCCGATAACGCCCGCCACAACGGCTAACAGATGAACCGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGCTTTCTACCGAGACCGACATCACCATCAGGAACACCAACTTTAGACATGCGGTGACGGTcctatacggcacattat          >  2:174855/1‑236 (MQ=255)                                                                                                                                                                                                                                          agagCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGGGACGGTCCCTTTCGGTAACTTATTg        >  1:310251/1‑251 (MQ=255)                                                                                                                                                                                                                                           gagCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCgaccgagac                                                                        >  2:346182/1‑186 (MQ=255)                                                                                                                                                                                                                                                tGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCCAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAAc  >  2:455144/1‑251 (MQ=255)                                                                                                                                                                                                                                                tGACGGCAGGATCGCGCAAAAATCAATGCTCCTCCGCAAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGACCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTACCTTTCGGTAACTTATTGGAAAAc  >  1:122627/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAAC  >  NC_000913/4567026‑4567515

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