Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,600	A→C	P141P (CCA→CCC)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,600	0	A	C	100.0%	167.9 / NA	51	P141P (CCA→CCC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/24);  total (27/24)

CCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTG‑CCTGGCGGGCGTAATAGCGAAATG  >  NC_000913/4638364‑4638837                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                ccTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGA‑GAACGCAAGATAGGGATTTATACCTGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGc                                                                                                                                                                                                                                  <  1:171384/251‑1 (MQ=255)     tgACCGAGCTTTATACGGTGCAGGAAGAGGATACAACCGTGGAGCGGAAACGAAGTTTTCTCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAGTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCAATTTCGAGGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGCCAAGCCCTTTATTGTGATTAGCGTc                                                                                                                                                                                                                               >  1:435018/1‑251 (MQ=255)            gCTTTATACTGTGAAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCGTTATTGTGATTAGCGTAGGGGATg                                                                                                                                                                                                                        <  2:31723/251‑1 (MQ=255)                     ggTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGCAGTTTTATCCATTTTTGGTTACCTGAGTCATTGCTGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGCTGCGCG‑GGGt                                                                                                                                                                                                               >  2:466985/1‑250 (MQ=255)                       ggCAGTAAGAGGATAAAGCGGTGGATCGGAAACGAAGTTGTATCCAGTTTTGGTTAGCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGAAAGATAGGGATTTATCCCGGTCAGGTCTGGCACAGGGATTTAACGCTAAAAGCTGATTTCGATGTTTCGCGTCGTAGCGAGCTCAACGCGGCAAATAGCATCTTAGGCAATCCCTTTATTGTGATTAGCGTCGGGGAGGCGCT‑TGGTagt                                                                                                                                                                                                            <  2:412619/250‑3 (MQ=255)                         cAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATtgg                                                                                                                                                                                                          >  1:122133/1‑251 (MQ=255)                           ggAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGCGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATtggtg                                                                                                                                                                                                        >  2:3361/1‑251 (MQ=255)                                agGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATtggtgtggtg                                                                                                                                                                                                   <  1:362077/251‑1 (MQ=255)                                      aaaCCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAgcg                                                                                                                                                                                             <  1:339553/251‑1 (MQ=255)                                       aaCCGTGGAGCGGTAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGATGAAGGCAATCAGAACGTAGCAGAACGCAAGATAGGGATTTATAACGGTCAGGTCTGGCAAAGTGATTTACCGCTAAAAGCCAATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATAATCTTAGTCAAGACCTTTTTTGTGATTAGCGTCGGCGATGCTCG‑TGATTTTGGTGAGGAGAATgcgc                                                                                                                                                                                            >  2:266832/1‑251 (MQ=255)                                                  ggAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGGGTGGTGGAAGGGCCTggaggtaa                                                                                                                                                                                  >  2:305282/1‑243 (MQ=255)                                                   gAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGAt                                                                                                                                                                                                                         <  1:212420/211‑1 (MQ=255)                                                   gAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGAt                                                                                                                                                                                                                         >  2:212420/1‑211 (MQ=255)                                                     aaggAAGTTTTATCCAATTTTGGTTACATGAGTCATTGATGGTTGATGGCAATCAGAACGTGGTAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTGACGCTAAAAGCCGATTTCGATGTTTCGCGTCTAAGCGAGCTCAACGCGCCAACTATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGAGTTGGGGTGGTGAAGGCGCCTGCAGTTaacgg                                                                                                                                                                               >  2:28315/4‑250 (MQ=255)                                                      acggaGTTTTATACATTTTTGGTTACCGGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGCACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCAAAGTAATTTAAAGATAAAGGCCTATTTCGATGTTTAGTGTCTAAGAGTGCTCAACGCGACAAATATCATCTTAGGAAAGACCTTTATTGTGCTTAGAGTCGGGGGTGCGCGGTGG‑ATTGGTGAGGTAAAAgagacgagatttagctga                                                                                                                                                                              >  2:261412/5‑233 (MQ=255)                                                        gAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaacggaacg                                                                                                                                                                           <  1:332330/251‑1 (MQ=255)                                                              ttATCAATTTTTGGTTACATGAGTCATTGATGGTTGATGGCAATAAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGGAAGAGCCTGAAGTTACCGGAACGGCGCTg                                                                                                                                                                     >  2:238391/1‑251 (MQ=255)                                                               tATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTgg                                                                                                                                                                    >  1:183988/1‑250 (MQ=255)                                                                       tttGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGAGTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGa                                                                                                                                                             <  2:452112/250‑1 (MQ=255)                                                                            tcaCCTGAGACATTGATGGTTGATGGCAATCAGAACGTGGTATAACGAAAGATAGGGATTTATACCGGTCAGGTCTGTCACAGTGGTTTAAAGATTAAAGACGATTTCGATGTTTAGCGTCTTAGCGATCTCAACGCGCAAAATATCATCATAGGCAAGCCCTTTATTGAGATTAGCGTCGGGGGTGAGAG‑TGGTATTGGTGTGGTGAAAGCGACGGGAGTTaacggagcagcgcataccatttaaacggg                                                                                                                                                       >  2:33505/3‑228 (MQ=255)                                                                             tACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGGGGTGAAAGCGCCTGAAGTTAACGGAACGGCCCTGACCATTGGACCCGGc                                                                                                                                                      >  2:91092/1‑251 (MQ=255)                                                                               ccTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGTAGTTAACGGAACGGCGCTGACCATTGAACCCggaac                                                                                                                                                    >  2:2537/1‑248 (MQ=255)                                                                                      cATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑GGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGtt                                                                                                                                              >  2:474511/1‑250 (MQ=255)                                                                                        ttgatggttgatggCAATCAGCACGTGGAAGCACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGAGATTTAACTCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGCTTAGCGTCGGGGAATCGCG‑TGGGATTGGTGTGGTGAAAGAGCCTGAAGTTAACGGGACGGCGCTGGCCCATGAGCACGGCACCGGGtgtgg                                                                                                                                           >  2:302763/1‑247 (MQ=255)                                                                                         tgatggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTag                                                                                                                                            >  1:190724/1‑249 (MQ=255)                                                                                          gatggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAg                                                                                                                                                                                                                                   <  2:127159/162‑1 (MQ=255)                                                                                          gatggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAg                                                                                                                                                                                                                                   >  1:127159/1‑162 (MQ=255)                                                                                                   tggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCgggg                                                                                                                                >  2:115586/1‑250 (MQ=255)                                                                                                    ggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCa                                                                                                                               <  1:270342/251‑1 (MQ=255)                                                                                                                gTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATGGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCTCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCatat                                                                                                                   <  2:435018/251‑1 (MQ=255)                                                                                                                    aagaaCGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTTGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTGGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGCATTGGTGTGGTGCAAGCGCCTGAAGTTAACGGACCGGCGCTGACCATTGACCCCGGCACCGGGTTAGAGCAAGGCGGGCCGGGCgggccaatccct                                                                                                               >  1:354183/1‑240 (MQ=255)                                                                                                                     agaaCGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCtt                                                                                                              <  2:236827/251‑1 (MQ=255)                                                                                                                              gATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAg                                                                                                                                                                                               >  1:228257/1‑161 (MQ=255)                                                                                                                              gATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAg                                                                                                                                                                                               <  2:228257/161‑1 (MQ=255)                                                                                                                               aTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAgg                                                                                                    >  1:77590/1‑251 (MQ=255)                                                                                                                               aTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATACGCG‑TAGCCTTGGTGTGGTGACAGCCGCTGCAGTTAACGGCACGGCGCTGACCATTGAACCCGGCCACGGGGTACAGCAAGGTGGGGACGGCGTGGTTATCCCTTTAccccaa                                                                                                      >  2:146644/1‑245 (MQ=255)                                                                                                                                         ataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCgg                                                                                          <  2:277960/251‑1 (MQ=255)                                                                                                                                         ataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCATCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCgg                                                                                          <  1:331689/251‑1 (MQ=255)                                                                                                                                               gTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCAAATCCCTTTACCCGAAGGGGACTGGAGGAAGCa                                                                                     >  1:360391/1‑250 (MQ=255)                                                                                                                                                 cAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCAGTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGaa                                                                                  <  1:302763/251‑1 (MQ=255)                                                                                                                                                         gCACAGTGAGTTAACGCTAAAAGCCGAGTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCGTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGc                                                                          <  2:485858/251‑1 (MQ=255)                                                                                                                                                                     aaCGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGccc                                                              >  1:374090/1‑251 (MQ=255)                                                                                                                                                                                ccGAGTTCGATGTTGCGCGGCTTAGCGCGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGTTGAATATGGCCCTGAATTTacgc                                                   <  2:259949/251‑4 (MQ=255)                                                                                                                                                                                   aTTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGt                                                <  1:186780/251‑1 (MQ=255)                                                                                                                                                                                     ttCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGGGCATATCCCTTTCCCCGAAGGGGACTGGCGGGAGCCGAACCTGAAGCTGGAAATGGCCCTGAGTTTAAACggaac                                              >  2:38607/1‑248 (MQ=255)                                                                                                                                                                                                           agcTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTG‑CCt                       <  1:243174/249‑1 (MQ=255)                                                                                                                                                                                                                     gcCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTG‑CCTGCCGGGCGTa             <  2:299555/251‑1 (MQ=255)                                                                                                                                                                                                                        aaaTATCATCTTAGGCAAGCCCTTTATTGGGATTAGCGTCGGGGATGCGCG‑TGGTAGTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTG‑CCTGCCGGGCGTAATa          <  2:122133/251‑1 (MQ=255)                                                                                                                                                                                                                            atcatcTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTG‑CCTGCCGGGCGTAATAGCg       <  1:474511/250‑1 (MQ=255)                                                                                                                                                                                                                              catcTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCAGATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTTTGGTG‑CCTGCCGGGCGTAATAGCGaaa    <  1:270641/248‑1 (MQ=255)                                                                                                                                                                                                                                 cTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGAAGAACCTGGAGCTGAATATGGCCCTGAATTTAAGCGGGACCGGCGATCCTTCTGTGGTGACCTC‑CGGGAGTAATAGCGGAATg  >  1:202596/1‑250 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                CCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTG‑CCTGGCGGGCGTAATAGCGAAATG  >  NC_000913/4638364‑4638837

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