Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,455	+G	intergenic (+589/+167)	rtcR → / ← glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,560,455	1	.	G	100.0%	218.4 / NA	60	intergenic (+589/+167)	rtcR/glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (28/32);  total (28/32)

ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTG  >  NC_000913/3560205‑3560693                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               aTCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGc                                                                                                                                                                                                                                                 >  1:591830/1‑251 (MQ=255)        tttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAt                                                                                                                                                                                                                                          >  1:530082/1‑251 (MQ=255)          tACCATCACCGTTTTAGCACCTTTGAGTTTGTTGGTGCCAATGCCCAATTTCCTTTGATTGAGCATTGCGTGCCCTCCCGCTTCCGGGTTGGTGCCTATATCGATAAACGGCGTCGAGCCATTGGGTATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                         <  1:330488/250‑1 (MQ=255)            ccATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                                                                                                                                                                                                      >  1:117363/1‑251 (MQ=255)            ccATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                                                                                                                                                                                                      >  2:376843/1‑251 (MQ=255)               tCAACGTGTTAGCAACGATGAGATTGTTGGTGACAATGCGCAACTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGttt                                                                                                                                                                                                                                   >  1:439912/1‑251 (MQ=255)                    gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                              >  1:293180/1‑251 (MQ=255)                               gTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCAt                                                                                                                                                                                                                   >  2:451834/1‑251 (MQ=255)                                    gaTTGTTGGAGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                                                                                                                                              >  1:629084/1‑251 (MQ=255)                                     aTTGTTGGGGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                                                                                                                             >  1:168033/1‑251 (MQ=255)                                       ggttgGGGAAAATGCGCAAATGGCTGTGATTGAGAAGTGGGTGCGCTCCCGCTTCCGGCGGGGTGCCGATAGCGATAAACAGCGTCGAGCCATGGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                            <  1:225150/249‑1 (MQ=255)                                          tgGTGAAAACGCGCAAATTGTTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCGGGCGTGGTGCCGATAGCGAAAAACAGCGTCGAGCCATGGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCGCGGGGT‑GTTACCCGAACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                                                                                                                         <  1:343583/250‑1 (MQ=255)                                          gggtgaaagtgcgaaatttacTGTGTTTGAGCGGTGCTGGCGTTCCCGCTTCGGGCGGGGTGCGGATCGCGAAAAACAGCGTCGTGCCTTTGGGGATTTGCCCCGCCCCGTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                                        <  2:5310/231‑1 (MQ=255)                                          gggggACAATGCGCAAATGGCTGTGATTGAGCAGTGCTTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                                        <  1:493573/247‑1 (MQ=255)                                               aCAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                                   >  1:549367/1‑251 (MQ=255)                                                 aaTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                                   >  2:385624/1‑249 (MQ=255)                                                    gcgcAAATTGCTGTGATTGAGCGGGGCGTGCGCTACCGCTTCCGGGGTGGTGCCTATATCGATAAACGGCTTCGAGCCTTTGGGGATTTGCTCCGCCCCTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                                                                                                                                              <  1:626233/251‑1 (MQ=255)                                                       cAAATTGCTGTGATTGAGCAGTGAGTGCGCTACCGCATCCGGCGTGGTGCCGATATCGATAAACAGCGGAGCGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAAGCAGCGCCTCACAGCCATGATGGCGCAGGAGCAGGTTTTGCTCAGCCAGCTCATTCAGGACGCGGCGAATAGTCTGCg                                                                                                                                                                                           >  2:585836/1‑251 (MQ=255)                                                            tgttgtgATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGGGGTGCGGATATCGATAACCGGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGAGGCGCTCTTTTTCTTCGGTCTGGTGG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                       <  1:527753/247‑1 (MQ=255)                                                              gtgggtttGAGCAGTGCGGGCGCTCCCGTTTCCGGCGTGGGGCCGATATCGATAAACAGCGTCGAGCCTTTGGGGTTTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                                                                                                                    <  1:132318/245‑1 (MQ=255)                                                                      gAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGGGCCGATATCGTTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                                                                                                            <  1:103698/251‑1 (MQ=255)                                                                       aGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                                                                                                                           >  1:454066/1‑251 (MQ=255)                                                                             gcgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAACCGGCGTCGGGCCATTGGGGTTTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCAATCGTGCCACGGCGG‑GTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                                                                                                                     <  1:340418/251‑1 (MQ=255)                                                                              cgtgccCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCCCCGCCACTTTG‑CGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCAGCCCGCTCATTGAGGTCGAGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                    >  1:425763/1‑251 (MQ=255)                                                                                           ttCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                                                                                                       >  1:270336/1‑251 (MQ=255)                                                                                                      gcacaattgcgaaaaacgtggTCGAGCCTTGGGGGTTTTCCCCCGCCCCTTGG‑CGGGCGATGCGCTTTTTTTCTTCGGTTTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGAAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                            <  2:357397/231‑1 (MQ=255)                                                                                                                     aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                             >  1:188038/1‑251 (MQ=255)                                                                                                                      cAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                            >  1:189056/1‑251 (MQ=255)                                                                                                                        tgatcggccctttgggtttttGCCCCGCCGCTTTG‑CGGGCGATGCCCTTTTTTTCTTCGGTCTGGGTG‑CCCTTGCGATCGTCCCCCGGCTT‑GTTAACCGAACTGGAGGGCACCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                           <  2:44454/232‑1 (MQ=255)                                                                                                                          gTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                           >  1:44454/1‑248 (MQ=255)                                                                                                                          gTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGATCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCgct                                                                                                                        >  1:486476/1‑249 (MQ=255)                                                                                                                            gcatcccttgggggttttcttccgcACTTTGGCGGGCGGTGCGCCGTTTTTCTTCGTTCTGG‑TGTGCCTTGCGATCGGCCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCTCACCGCCATGATGGCGCAGGTTCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACTGATAAATGCTCTACCAGCTCTTCGGTACTGACATAACTCTGCGGTTTAACCAGTTCGata                                                                                                                       <  2:501434/227‑1 (MQ=255)                                                                                                                              aGCCATTGGGGCTTTGCTCAGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCGGAGGGCTGACGGAGAAATGCTCTACCAGCTCTTAGGTACTGACATAACCCTGCTGTTTTACCAGTTCGataat                                                                                                                     >  2:518528/1‑250 (MQ=255)                                                                                                                                  aTTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                                                                                                                >  1:327882/1‑251 (MQ=255)                                                                                                                                          ttggccccGCCTCTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg                                                                                                         <  2:454066/244‑1 (MQ=255)                                                                                                                                            gggcccGCCCCTTTG‑CGGGCGTGGCTCTTTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTACCCGACCTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg                                                                                                      <  2:591830/247‑1 (MQ=255)                                                                                                                                                cccgcctgttggggggggttGCGCTCTTTTGCGTCGGTCTGGGTG‑GCCTGGCGATCGTGCCACGGGGT‑TTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCATGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt                                                                                                   <  2:480797/232‑1 (MQ=255)                                                                                                                                                 cGCCTCTTTG‑CGGGCTATGCGCTCTTTTTCTTCGGTCTGGGGG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                                 <  2:87827/251‑1 (MQ=255)                                                                                                                                                    cACTTTG‑CGGGTGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCTATCGTCCCACGGCGT‑GTTACCCGAACGGGAAGGCAGCGCTGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                              <  1:25587/251‑1 (MQ=255)                                                                                                                                                    cACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTTTGGGTG‑GCCTGGCGATCGTGCCACGGGGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                              <  1:479940/251‑1 (MQ=255)                                                                                                                                                             gggCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                      <  1:385624/251‑1 (MQ=255)                                                                                                                                                              ggCGCTGCGCTTTTTTTCTTCGGGCTGGGTG‑GCTTGGCGATCGTCCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTCTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                     <  1:340162/251‑1 (MQ=255)                                                                                                                                                                 gATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAt                                                                                  <  1:234254/251‑1 (MQ=255)                                                                                                                                                                           ttttCTTCGGTCTGGGTG‑GCCTGGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAat                                                                         <  1:490900/250‑1 (MQ=255)                                                                                                                                                                                 tcGGTCTGGTGG‑GACTTGCGATCGTGC‑ACCGTGTCGTTTACCGAACTTGATTGCGCCGCCGCACCGCCATGATGGGGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATATCCTGCGGGCTGACGGATAAATGCTCTACCATCTCTTCGGTACTGACATATCCCTGCTGTTTAACCAGGTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCtttaattctttt                                                                   <  2:579298/250‑12 (MQ=255)                                                                                                                                                                                               ccTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                                     >  1:13254/1‑251 (MQ=255)                                                                                                                                                                                                  ggCGACCGTGCCCCGGCGT‑GTGACCCAAACGGGAAGGCAGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAg                                                  <  2:592172/250‑1 (MQ=255)                                                                                                                                                                                                     gATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCg                                                >  1:429321/1‑250 (MQ=255)                                                                                                                                                                                                      aTCGTGCCACGGCGT‑GTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGaa                                              <  1:84791/251‑1 (MQ=255)                                                                                                                                                                                                        cGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGACAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                            >  2:396052/1‑251 (MQ=255)                                                                                                                                                                                                                ggCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTATACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                                      >  1:460898/1‑249 (MQ=255)                                                                                                                                                                                                                         aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCt                            <  1:130789/251‑1 (MQ=255)                                                                                                                                                                                                                           gcgatcgggAAGGGAGCGCAGCTCCGCCATGAGGGCGCAGGATCAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTGCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTaa                          <  2:74936/244‑1 (MQ=255)                                                                                                                                                                                                                             gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAc                        >  1:242471/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg                   >  2:231564/1‑251 (MQ=255)                                                                                                                                                                                                                                     gggtaGCGCCGCACCTCCATGATGGCGCAGAATGAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCTGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTACCCAGTTCGATAATCCCTTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGTAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCt                <  2:81917/247‑1 (MQ=255)                                                                                                                                                                                                                                     aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCt                >  2:424548/1‑251 (MQ=255)                                                                                                                                                                                                                                      gccAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa               <  2:117363/249‑1 (MQ=255)                                                                                                                                                                                                                                        ccgcgcCGCACCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGGTTCATTTATAAATCCCTGGAATTATGTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAc             <  2:439912/249‑1 (MQ=255)                                                                                                                                                                                                                                                  ccGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGTCTAACCCGGCGAt     <  2:530082/249‑1 (MQ=255)                                                                                                                                                                                                                                                   cGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCAAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtg  >  1:623091/1‑251 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTG‑CGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTG‑GCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTG  >  NC_000913/3560205‑3560693

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