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%	198.9 / NA	55	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 (27/28);  total (27/28)

TTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTG  >  NC_000913/3560213‑3560693                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               ttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                         >  1:370753/1‑251 (MQ=255)  tACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                         >  1:64945/1‑250 (MQ=255)     cATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACACAGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGt                                                                                                                                                                                                                                     >  1:62392/1‑251 (MQ=255)     cATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGt                                                                                                                                                                                                                                     >  1:152717/1‑251 (MQ=255)         cgcccgaccccaaggtggagtttgttgGTGAAAAGGCGCAAATTGGTTTGATGGAGCAGTGCGTGCGCTCCCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGGGTTTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGGGGCCTTGCGATGGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGtttt                                                                                                                                                                                                                                  <  1:188642/229‑1 (MQ=255)            gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                              >  2:452725/1‑251 (MQ=255)             tgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTcc                                                                                                                                                                                                                             >  1:374013/1‑251 (MQ=255)                   cAACGTTGAGATTGTTGGTGACCATGCGCAAATTGCTGTGCTAGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGtttt                                                                                                                                                                                                                                  >  2:188642/1‑240 (MQ=255)                               tgttgGTGAAAATGCGAAAATTGCTGTGATGGAGCAGTGCGTGCGCTCCCGCTTCCGGGGTGGTGCGGATACCGATAAACAGCGTCGAGCCATGGGGGATTGGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                            <  1:174882/250‑1 (MQ=255)                                gttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                                                                                                                          >  1:125175/1‑251 (MQ=255)                                  tggtgatcaTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                                        >  1:129346/9‑251 (MQ=255)                                    gTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                                                                                                                      >  1:309094/1‑251 (MQ=255)                                                aaattgttgtgatggggcgttgcgggcgCTACCGCTCCGGGCGGGGTGCGGATTTCGATAAACAGCTTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                          <  1:338283/231‑1 (MQ=255)                                                                              cccGCTTCCGGCGTGGTGCCGATATCGATAACCAGCGTCGGGCCTTGGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                            <  1:97773/250‑1 (MQ=255)                                                                                         cGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                                 >  1:258200/1‑251 (MQ=255)                                                                                              tgCCGATATCGATAACCAGCGACGAGCCCTTGGGGATTTGCTCCGCCAC‑ATTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                            >  2:242745/1‑251 (MQ=255)                                                                                               cacgtgatcaaaaaatgccggtgcgCCATTGGGGATTTTCTCCGC‑ACGTTGGCGGGGGATGCGCTCTTTTTCTCCGGGCTGGGTGGCGTGGCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCTGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGTCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                           <  2:152717/227‑1 (MQ=255)                                                                                                       cgatacaaaGCGTCGAGCCATAGGGGATTTGCTCCGCCAA‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGACTGTGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTAATTAAGGTCGCGGCGAATACTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                                    >  2:484774/9‑250 (MQ=255)                                                                                                             aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                             >  1:113216/1‑251 (MQ=255)                                                                                                               ggCGTCGAGCCATGGGGAATTTGCTCCGCCCC‑TTTGCGGGCTATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                           <  1:400030/250‑1 (MQ=255)                                                                                                               gcggtggtgCCAGTGGGGATGTGCCCCGCCTC‑TTTGCGGGCGATGCGCTCTTTTTCGCCGGTCTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                           <  1:452725/243‑1 (MQ=255)                                                                                                                gggtcgagccgggggggTTTTGCTCCGCCAT‑GTTGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGTGCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCGGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                                                                                                                          <  2:213008/238‑1 (MQ=255)                                                                                                                 ggTCGACCCTTTGGGGTTTTGCCCCGCCAC‑TTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                                                                                                                         <  1:48355/250‑1 (MQ=255)                                                                                                                         cATTGGGGTTTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                                 <  1:281389/251‑1 (MQ=255)                                                                                                                         cATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                                 >  1:245386/1‑251 (MQ=255)                                                                                                                                gatgtGCTCCCCCAC‑TTTGTGGGCGATGCGCTCTTTTTCTGCGGTCGGGGTGGCCGTGCGATCGGGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTTCCAGCTCTGCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAc                                                                                                          <  2:374013/247‑1 (MQ=255)                                                                                                                                  tttcccccGCCCC‑TTTGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTGGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg                                                                                                         <  2:309094/244‑1 (MQ=255)                                                                                                                                     gCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGCGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCACGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGCAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGCTGTGACGTtgt                                                                                                     >  2:21262/1‑251 (MQ=255)                                                                                                                                          gCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCATGCTGTTTAACCAGTTCGATAATACCGTTGTGCCGTTGtgtttg                                                                                                >  1:140279/1‑251 (MQ=255)                                                                                                                                                        ggATGCGCTTTTGTTCTGCGGGCTGGGGGGCCTGGCGAGCGGCCCAGGGCGTGTTAACGGAACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                                   <  2:184417/250‑1 (MQ=255)                                                                                                                                                          atgctctgtgttgctcgggTCGGGGTGGCCGGGCGATGGTGCCCCGGCGTGTGAACCGAACGGGAGGGCAGCGCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATc                                                                                 <  2:385575/234‑1 (MQ=255)                                                                                                                                                                tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTggg                                                                           >  2:207364/1‑250 (MQ=255)                                                                                                                                                                 tttttttttCGGTCGGTGTGGCCGTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGa                                                                           <  2:161094/249‑1 (MQ=255)                                                                                                                                                                   tattcttcAGTCTGTGTGGCCTTGCGATCGTACCACGGCGTGTTCACCGAAGTGGAAGGCGGCGCCGCACCTCCATGATGGCGCAGGATCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTTCGTGCTGACTGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAat                                                                         <  1:155591/248‑1 (MQ=255)                                                                                                                                                                          cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttc                                                                  <  1:414054/250‑1 (MQ=255)                                                                                                                                                                           ggTCTGGGGGGCCTGGCGATCGTGCCACGGGGTGTTAGCCTAATGGGAAGGCAGCGCCGGACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCATTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                <  2:158075/251‑1 (MQ=255)                                                                                                                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACAGGAAGGCCGCGCCGCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtt                                                               >  1:235305/1‑251 (MQ=255)                                                                                                                                                                              cTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                             >  2:225175/1‑251 (MQ=255)                                                                                                                                                                                  gTGGCCTTGCGATCGGCCCACGGCGTGTTAACCGACCGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGAGAAATGCTCTATCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgcg                                                         <  1:473153/251‑1 (MQ=255)                                                                                                                                                                                   gggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                                        <  1:47900/250‑1 (MQ=255)                                                                                                                                                                                    ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCa                                                       <  1:439603/251‑1 (MQ=255)                                                                                                                                                                                          gCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGc                                                 <  2:125175/251‑1 (MQ=255)                                                                                                                                                                                                gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                            >  1:87883/1‑250 (MQ=255)                                                                                                                                                                                                 tGCCACGGCGTGTTAACCGAACTGGAAGGCAGAGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATca                                           >  2:431403/1‑250 (MQ=255)                                                                                                                                                                                                  gcctgggctggtTAACCGTACGGGACGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGGCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAGCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaa                                          <  2:181308/240‑1 (MQ=255)                                                                                                                                                                                                   ccatggtgttgtAACCTAACGGGAAGGCCGCGCCTCCCCGACATGATGGCGCAGGATATGGCTTCGCCTCGCCATCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTATCCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGTGCATTGAGCGAATcaaca                                        <  2:29/240‑1 (MQ=255)                                                                                                                                                                                                                cccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTa                           <  1:105966/250‑1 (MQ=255)                                                                                                                                                                                                                   gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAc                        >  2:331265/1‑251 (MQ=255)                                                                                                                                                                                                                    aaCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAACAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTTTTTTCGTTTTCGCGCCTTGAGCGAATCTACAAACGCCATCGCTAAAcc                       >  2:87317/1‑251 (MQ=255)                                                                                                                                                                                                                      cgggaaggaaggacctgtccGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCa                     <  1:459047/233‑1 (MQ=255)                                                                                                                                                                                                                            ggCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa               <  2:118528/251‑1 (MQ=255)                                                                                                                                                                                                                               aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc            >  1:138090/1‑251 (MQ=255)                                                                                                                                                                                                                                    cgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                                      >  1:476014/1‑220 (MQ=255)                                                                                                                                                                                                                                    cgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCg       >  2:178525/1‑251 (MQ=255)                                                                                                                                                                                                                                         cgtcATGATGGCGCAGGATCAGGTGTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCTCGGCTAACCCGGCGAtgtg  <  2:53898/248‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               TTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAC‑TTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTG  >  NC_000913/3560213‑3560693

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