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.4%	113.5 / ‑3.9	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 (1/0);  new base A (29/9);  total (30/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

TATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAAC  >  NC_000913/4567027‑4567515                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     gaTGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATc                                                                                                                                                                                                                                                 <  1:108140/250‑1 (MQ=255)  aTGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATcg                                                                                                                                                                                                                                                >  1:251111/1‑251 (MQ=255)    ggCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATcgcg                                                                                                                                                                                                                                              >  1:270873/1‑251 (MQ=255)        agaAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCaa                                                                                                                                                                                                                                           >  1:490078/1‑250 (MQ=255)                           aCGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGCTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATAACTGCTACTCCACtc                                                                                                                                                                                                                       >  2:375510/1‑250 (MQ=255)                            cgcgGCAAAACAGCTCCCCCCCCGCAAACCATTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAg                                                                                                                                                                                                                      <  2:38312/251‑1 (MQ=255)                               ggcgAAACAGCGCCACCACCGCAAACCAGGCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTc                                                                                                                                                                                                                   <  1:88490/251‑1 (MQ=255)                                     aCAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCACAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGc                                                                                                                                                                                                             >  1:37706/1‑251 (MQ=255)                                            caccacCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGc                                                                                                                                                                                                      >  1:115427/1‑251 (MQ=255)                                                accGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGgc                                                                                                                                                                                                  >  1:18391/1‑251 (MQ=255)                                                                ccAGCTCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCaa                                                                                                                                                                                  <  2:270873/251‑1 (MQ=255)                                                                    cgcgcCGACCATCGCCGCGTCAGCAATCGCCTACACGCCGCGCACCCAAAAATTGGGCGTCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGAGTAGTTCATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACGGCTCCTCCACTAGCTCGGCTGCCTGATGCGGGCGCGCGCGATGGAGCAATccc                                                                                                                                                                              <  2:29920/251‑1 (MQ=255)                                                                       gcCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGcc                                                                                                                                                                           >  1:375816/1‑251 (MQ=255)                                                                         cgacacccgccgtgtCAGCAAGCACCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCGTTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGAGCGCGCAAAAAACACGGCTCCTCCGCTAGCTCGGTTGCCTGATGCTGGCGCACGCGTTTGAGCAATCCCGCCAg                                                                                                                                                                         <  2:67777/237‑1 (MQ=255)                                                                          accctcggccttttcacaaaccGCCTCCACGCCGCTAACCCAAAAATTGGTCGGAAAAAACAGCGTAACGACAAAGTGAGCAGGGGCGATATGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGGGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCGGGATCGCGAAAAAATCACTGCCCCTCCTCTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAgg                                                                                                                                                                        <  2:394435/230‑1 (MQ=255)                                                                           aCCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAggg                                                                                                                                                                       >  1:195303/1‑251 (MQ=255)                                                                            ccttctccgcTACAGCAACCTCCTACACGCCGCTCACCCAAAAAATGGTCGGCAAAAACACCGATACGACAAAGGTAGCAGCGGCGATTAGTATTAAAAAGAACTCAAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATACATAACGTAAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGTCGCACGCGAATGTTCAATCCCGCCAGGGt                                                                                                                                                                      <  2:425739/245‑1 (MQ=255)                                                                                      ttCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATACCGCCAGGGTGAAAacgacg                                                                                                                                                            >  1:130471/1‑251 (MQ=255)                                                                                         aGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCCCTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCTATTGAGCAATCCCGCCAGGGTGAAAACGACGATa                                                                                                                                                         <  1:162565/251‑1 (MQ=255)                                                                                            aaTCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAAcg                                                                                                                                                      >  2:184496/1‑251 (MQ=255)                                                                                              tCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAAcgcc                                                                                                                                                    >  2:64560/1‑251 (MQ=255)                                                                                                gCCTTCACGACGCTCATCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTCAAGAGAGCGCCAACCTTTTGGCACGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACCGCGCCTCCACTAGCTCGGATGCCTGATGCTGGCACAAGCGATCGGCCAAACCCGCCAGGGCGGAAACGACGCTAAagcccc                                                                                                                                                  >  2:447282/1‑250 (MQ=255)                                                                                                            cTCACCCCAAAATTGGGGGGCAAAAACAGCGTAACCACAAAGGTAGCAGCGGCGATAAGCAGTCAAGAGAGCGCCAACCTTTTGGCGCGGCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCAATGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAATACGGCGATAACGCCCGCCACAACGGCt                                                                                                                                       >  1:263728/1‑250 (MQ=255)                                                                                                                   aaaaaTTGGGCGCCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTCATGCTGGCGCACGCGATTGAGCAATCCCGCCCGGGGGAAAACGACGACAACGCCCGACAAAACGGCTAACAgct                                                                                                                                >  2:193344/1‑248 (MQ=255)                                                                                                                                  aaaaCAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTg                                                                                                                                                                                                    <  2:125448/167‑1 (MQ=255)                                                                                                                                  aaaaCAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTg                                                                                                                                                                                                    >  1:125448/1‑167 (MQ=255)                                                                                                                                           tAACGACAAAGGTAGAAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTaa                                                                                                       >  1:273038/1‑250 (MQ=255)                                                                                                                                               gACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCAt                                                                                                  >  1:503092/1‑251 (MQ=255)                                                                                                                                                aCAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATACCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCAt                                                                                                  >  1:109125/1‑250 (MQ=255)                                                                                                                                                                    aaGCAGTAAAGAAAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGAGAAAACGCCGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAACCCATAGCCAGCGGCCGATATCgccc                                                                             >  2:717/1‑248 (MQ=255)                                                                                                                                                                         gTAAAGAGAGGGCCAACCTTCTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGCATAGAGCCTGACGACAGGATCGCGCAAAAGTCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCACTCCCGCCAGGGTGAAAACGACGAGAACGCCCGCCACAACGGCTCACAG‑ATGACACGGCAAACGGCCGCCGTTAAACCATACCAAGCGGccaaattcgcctgagg                                                                         >  2:206153/1‑236 (MQ=255)                                                                                                                                                                             agagagCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTCAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTAGACCCAGACCGAc                                                                    >  2:92173/1‑251 (MQ=255)                                                                                                                                                                                         ccTTTTAGCGCGACTGAGTTCGATGAGTTTATTAATAACGTTAAGACTAGAGCATGACGGCAGGATCGCGCAAAAATCACTGAACCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGAGTGAGCAATCCCGCCAGGGTGAAAAAGACGATACAGCCCGCCACAACGGCTA‑CAGGAAGAAACGGCAAGCGCACGCCGTACACCCAAAGCCAGAGGCCGCATTCGACCGCGACAGACATCACCATCaa                                                         >  2:307660/1‑249 (MQ=255)                                                                                                                                                                                             ttGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACACCGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTCCGACCGAGACCGACATCACCATCAGGATc                                                     >  1:83175/1‑250 (MQ=255)                                                                                                                                                                                                               aTGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGCCCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCCTCAGGATCACCATATTTAGCGATGCgg                                  >  2:110372/1‑251 (MQ=255)                                                                                                                                                                                                                                  tAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCCCGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTCACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGCTATCGACCGGTACCGACATCACCATCAGGATCACCAGATATAGACATGCTGACACGGTCCAtttgagta                >  2:14297/1‑245 (MQ=255)                                                                                                                                                                                                                                     gAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACtt            >  1:218497/1‑251 (MQ=255)                                                                                                                                                                                                                                      aaTAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACttt           >  1:411191/1‑250 (MQ=255)                                                                                                                                                                                                                                            gCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGa      >  1:406296/1‑250 (MQ=255)                                                                                                                                                                                                                                                gACGGCAGGATCGCGAAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAAACATAGCCAGCGGCCGATTTCGACCGCGACCGCCATACCCATCCGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGCAACTTATTGGaaaaa  >  2:299208/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     TATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAG‑ATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAAC  >  NC_000913/4567027‑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: ‑