Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	197,755	G→T	G404W (GGG→TGG)	rseP →	intramembrane zinc metalloprotease RseP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	197,755	0	G	T	100.0%	231.3 / NA	71	G404W (GGG→TGG)	rseP	intramembrane zinc metalloprotease RseP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (42/29);  total (42/29)

GCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGT  >  NC_000913/197511‑197997                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   gCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATc                                                                                                                                                                                                                                              <  1:337928/251‑1 (MQ=255)    gTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATctgc                                                                                                                                                                                                                                           >  1:153031/1‑251 (MQ=255)    gTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATctgc                                                                                                                                                                                                                                           >  1:350792/1‑251 (MQ=255)       cAACGCCATCGTCGAAGACACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGt                                                                                                                                                                                                                                        >  2:674591/1‑251 (MQ=255)             cATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCACGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTg                                                                                                                                                                                                                                  >  1:302365/1‑251 (MQ=255)             cATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTg                                                                                                                                                                                                                                  >  1:575624/1‑251 (MQ=255)                       gCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGaaaa                                                                                                                                                                                                                        <  2:399529/251‑1 (MQ=255)                        ccACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATACTTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGaaa                                                                                                                                                                                                                         >  2:379071/1‑249 (MQ=255)                                     gTGGCAGCTGATGAAGCGGACGGTTAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTGTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCCGTTTCCGTTGACCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCgtac                                                                                                                                                                                                          <  2:87991/251‑4 (MQ=255)                                         cAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGCGATGTGAAAATGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGt                                                                                                                                                                                                      >  2:536774/1‑251 (MQ=255)                                              gatgaAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCGCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCg                                                                                                                                                                                                 <  1:673038/251‑1 (MQ=255)                                              gatgaAGCTGACGGGCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGTACCGGTATcc                                                                                                                                                                                                  <  2:650659/250‑1 (MQ=255)                                                tgaAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCg                                                                                                                                                                                                 >  2:648020/1‑249 (MQ=255)                                                 gaAGCTGCCGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGCTCGAAAAGATCAAGGGCGGACCGGTATCCGAGc                                                                                                                                                                                              >  2:374072/1‑251 (MQ=255)                                                 gaAGCTGACGGTCAGTATGCTGGGAAAATTGATCACAGGTGATGTGAAACTGAACAACATCAGGGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCATGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCAGTATCCgtgc                                                                                                                                                                                              >  2:203483/1‑248 (MQ=255)                                                 gaAGCCGACGGTCAGTATGCCGGGAAAATTGACCACTGGTGATGAGAACCTGAACAACCTCAGTGGGCCGATCTCTATAGACAAGGGGGCTGGGATGACAGCGGGACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGATGGCATCTGCTGTTCCTTGCGATCGAACAGTTCAAGGGCGGACCGGTATCCGAGc                                                                                                                                                                                              >  2:48405/1‑251 (MQ=255)                                                            tCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGCAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGAc                                                                                                                                                                                   >  2:238018/1‑251 (MQ=255)                                                               gTATGCTGGGAACATTGAAAACCGGTGATGTGAAACTAAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGTGTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCATGAACTTCGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCGTCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCgg                                                                                                                                                                                            >  2:627186/1‑239 (MQ=255)                                                               gTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCGTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGAATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCgg                                                                                                                                                                                            <  1:627186/239‑1 (MQ=255)                                                                         aaaaTTGCTCACCGGTGATGTGCCACTGAACACCCTCCGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGCCAGCGGCACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACAGGAATCCGAGCGGGTTCAAGACTTTTGTTAtagc                                                                                                                                                                       >  2:406604/1‑247 (MQ=255)                                                                          aaaTTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCAt                                                                                                                                                                     >  2:612515/1‑251 (MQ=255)                                                                                  cACCGGTGATGGGAAACTGAACAACCTCAGGGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACATCGGAACTCGGGGTGGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGGACTTGACGGGTGGCATCTGCTGTGCCTTGCTATCTAAAATGTCAAGGGCGGGCCGGTATCCGAGCGGTTTCAAGACTTGTGTTATCGCATTGGCTCg                                                                                                                                                              <  2:672194/250‑1 (MQ=255)                                                                                   cccGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGAt                                                                                                                                                            <  2:456273/250‑1 (MQ=255)                                                                                             tgAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGCTCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTg                                                                                                                                                  >  1:181823/1‑251 (MQ=255)                                                                                              gAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTg                                                                                                                                                  <  2:350792/250‑1 (MQ=255)                                                                                                    gAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGtt                                                                                                                                             <  1:341351/249‑1 (MQ=255)                                                                                                     aacaacCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAAt                                                                                                                                          >  1:478683/1‑251 (MQ=255)                                                                                                      acaacCCCAGTGGGCCGATCTTTATCGCCAAGGGGGCTGGGATGACGTCGGAACTCGGGGTTTTGTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACTGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTGTTGTTATCGCATTGGCTCGATCCTGCTGGTGCTGTTAATg                                                                                                                                         <  2:74066/251‑1 (MQ=255)                                                                                                          ccTCAGTGGACCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGGGAACTTAGGGATAATTAACCTGTTTCCGTAGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGGGCTGTTAATGGGGc                                                                                                                                     >  1:553525/1‑251 (MQ=255)                                                                                                             cAGTGCGCCGACCTCTATCACCAAGGGGGCTGGGATGACCGCGGAACTCGGGGTTGTTTATTACCTGACGTTTCTTGCGCTTATTCGCGCGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGCTGGCATCTGATGTTCCTTGCGATCGAACAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTg                                                                                                                                  >  2:643754/1‑251 (MQ=255)                                                                                                                                 ccaggggggCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTtctct                                                                                                              <  2:444118/247‑1 (MQ=255)                                                                                                                                  cAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACGTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAAGGATTgctctc                                                                                                             <  2:478683/251‑1 (MQ=255)                                                                                                                                   aaGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCg                                                                                                            <  1:612515/251‑1 (MQ=255)                                                                                                                                        ggCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTAt                                                                                                       <  1:304525/251‑1 (MQ=255)                                                                                                                                            gggATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATga                                                                                                     <  1:544497/249‑1 (MQ=255)                                                                                                                                            gggATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACGTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATgaga                                                                                                   <  1:705058/251‑1 (MQ=255)                                                                                                                                               gtGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTGCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACGTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGagtt                                                                                                <  2:497574/250‑1 (MQ=255)                                                                                                                                               aTGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGagtt                                                                                                <  1:674591/251‑1 (MQ=255)                                                                                                                                                 gACAGCGGAACTCGTGGTTGGTTATTACCTGCCGGTTCTTGCGCTTATTAGCGTGAACGTAGGGATAATTAACCTGTGTCCGTTGCCCGTAATTGACGGGTGGCAGCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCTGATCGTTATCCGATCTGGTTCAATACTTTTGTTATCGCATTGTCTCGATTCTGCTGGTGCTTTGAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGagttag                                                                                              <  1:688297/251‑1 (MQ=255)                                                                                                                                                  aCAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGagttagt                                                                                             >  1:578856/1‑251 (MQ=255)                                                                                                                                                          aCTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGga                                                                                        >  2:633644/1‑248 (MQ=255)                                                                                                                                                              ggggTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTATTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGAGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACg                                                                                  >  1:208781/1‑250 (MQ=255)                                                                                                                                                                     tttATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGACGAACGCataata                                                                           >  1:411717/1‑250 (MQ=255)                                                                                                                                                                     tttATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCataata                                                                           >  2:517241/1‑250 (MQ=255)                                                                                                                                                                     tttATTACCTGCCGATTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGTTTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGGGAGTTAGTTAGGAAGAACGCataata                                                                           >  2:191693/1‑250 (MQ=255)                                                                                                                                                                       tattaCCTTCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCGCTCGGTTATGAGAGTTAGTTATGAAGAACGCATAATAAc                                                                         <  2:419368/250‑1 (MQ=255)                                                                                                                                                                        attaCCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCTTTGCCCGTTCTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTAGCGCATTGGCTCGTTTTTGCTGGTGCTGTTAATGGGGCTTGCACTTTGCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGa                                                                       <  1:447921/251‑1 (MQ=255)                                                                                                                                                                          taCCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGCACGCATAATAACGATg                                                                     >  2:406133/1‑251 (MQ=255)                                                                                                                                                                               gCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGAGTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCg                                                                  <  2:295506/249‑1 (MQ=255)                                                                                                                                                                                ccGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATg                                                               >  1:453009/1‑251 (MQ=255)                                                                                                                                                                                            cTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCt                                                    >  1:657487/1‑250 (MQ=255)                                                                                                                                                                                                  aGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCg                                             <  1:238018/251‑1 (MQ=255)                                                                                                                                                                                                    cGTGATCTTAGGGATACTTAACCTGTTTCCGTTGCCAGTACTTGACGGGTGGAATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGCCCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTa                                                                                                                                            >  2:673465/1‑154 (MQ=255)                                                                                                                                                                                                    cGTGATCTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTa                                                                                                                                            <  1:673465/154‑1 (MQ=255)                                                                                                                                                                                                       gAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCgct                                        >  2:600200/1‑251 (MQ=255)                                                                                                                                                                                                         aCTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCgctgc                                      >  1:131000/1‑251 (MQ=255)                                                                                                                                                                                                            tAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCTATGAAAAAGTTGCTCATAGCGTCgctgct                                     <  1:197411/249‑1 (MQ=255)                                                                                                                                                                                                             aGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTGATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGtt                                  <  2:302365/251‑1 (MQ=255)                                                                                                                                                                                                              gggATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGttt                                 >  1:597352/1‑251 (MQ=255)                                                                                                                                                                                                                  tAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTagc                              >  1:398128/1‑250 (MQ=255)                                                                                                                                                                                                                     ttAACCTGTTTCGGTTGCCAGAACTTGACGGGTGGGAACTGCTGTTCCTTGCGATCGACGAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGAATTTTGATATCGCATTGGCTCGATTCTGCTGGTGCTGATAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGGTAGGAGAGTTAGGTAGGAAGAACGAATAATAACGATGGGGATGAAAAAGTTGCTCATAGCGTCGCGGCTGTTTAGAAgcg                          >  2:58797/1‑251 (MQ=255)                                                                                                                                                                                                                       agccTGTTTCCGTTGCCCGTACTTGACGGGTGGCCTCTGCTGTTCCTTGCGCTCGAAAAGATCAAGGGCGGACCGGGATCCGAGCGGGTTCACGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCCCTTTTAAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGACGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCAGCTGTTTAGCAGCGcc                        >  2:91956/3‑251 (MQ=255)                                                                                                                                                                                                                       aaCCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGcc                        <  2:14883/251‑1 (MQ=255)                                                                                                                                                                                                                        aCCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGcc                        >  1:106884/1‑250 (MQ=255)                                                                                                                                                                                                                        aCCTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACAGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGcc                        >  2:254193/1‑250 (MQ=255)                                                                                                                                                                                                                         ccTGTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCAc                      >  2:588003/1‑251 (MQ=255)                                                                                                                                                                                                                            gTTTCCGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGt                   >  1:568791/1‑251 (MQ=255)                                                                                                                                                                                                                                 cGTTGCCCGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACg              >  2:691132/1‑251 (MQ=255)                                                                                                                                                                                                                                      cccGTACTTGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCt         <  1:257533/251‑1 (MQ=255)                                                                                                                                                                                                                                            gttGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTCTCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGGGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAggg   >  2:184113/2‑251 (MQ=255)                                                                                                                                                                                                                                             ttGACGGGTGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGt  >  2:396470/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   GCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGATCGAAAAGATCAAGGGCGGACCGGTATCCGAGCGGGTTCAAGACTTTTGTTATCGCATTGGCTCGATTCTGCTGGTGCTGTTAATGGGGCTTGCACTTTTCAATGATTTCTCTCGGTTATGAGAGTTAGTTAGGAAGAACGCATAATAACGATGGCGATGAAAAAGTTGCTCATAGCGTCGCTGCTGTTTAGCAGCGCCACCGTATACGGTGCTGAAGGGT  >  NC_000913/197511‑197997

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