Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,925	G→C	R77P (CGC→CCC)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,925	0	G	C	100.0%	233.3 / NA	69	R77P (CGC→CCC)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (35/34);  total (35/34)

GCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACG  >  NC_000913/4636680‑4637164                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                gCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACTCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGt                                                                                                                                                                                                                                            <  1:45242/251‑1 (MQ=255)      aGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTcgc                                                                                                                                                                                                                                       >  1:455042/1‑251 (MQ=255)       ggggCCTGTAAGGCGTAGCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACGGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTcgcg                                                                                                                                                                                                                                      <  1:396333/251‑1 (MQ=255)           ccTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGCGGCGGTGGCAGCCTGGATCGTAATGGCCATTTTTGTCAAAGAAATTAAACAGGGCGTG‑GGAAGAGCAATGGAGGGGACGTTGATAGACACCACAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGAAACTGGCC‑CAGGCGTTTAATCCGCTACAACA‑TACCCCGTTTCGCGCCa                                                                                                                                                                                                                                   >  1:442470/1‑250 (MQ=255)              gTAATGCGTGTCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGTTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACTCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACACCA‑TCCCCCGTTTCGCGCCAat                                                                                                                                                                                                                                 <  2:270879/249‑1 (MQ=255)                    cGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTTTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGCGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGGGGGTCCCGATTTGTTCTCTGGGGACCCAACGCATGGGCAACTGGGC‑CAGGCCTTTAATCAGCT‑CCACACTCCCCCGTTTCGAGCCACAATCGGTgg                                                                                                                                                                                                                         >  2:315669/1‑251 (MQ=255)                     gTATCGGCATGCGGTTGGTGCTGGGCTATTTTTTACTGGGGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACCCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATCGGGAACTGGCGGCAGG‑GTTTAATCAGCTGAAACA‑TCCCACGTTTCGCCCCAATATCGGGGGc                                                                                                                                                                                                                        >  2:272019/1‑251 (MQ=255)                          ggtatggggTTGTGGCTGGGCTATTTTTTCCTGGTGTCGTTGCCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGGT‑CGAAGATCAACGGAGTGGACGTTTACCGACACCGAAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCGGGGTACTCAACACATGGGCAACGTGCG‑TAGGCTTTTAATCAGGTTCAATA‑TCCCCCGTTTCGCGCCACTAGCGGTGTCATTaa                                                                                                                                                                                                                   <  2:171884/244‑1 (MQ=255)                            cATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCCTTAAcc                                                                                                                                                                                                                 >  2:48254/1‑250 (MQ=255)                                ccttcccgtcgctcttccgatctACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGTGACCCAACGCATGGGCAACTGGCG‑CAGGCGTGTAATCAGCTACAACA‑TCCCCCGTTGCGCGCCAATATCGGTGGCATTAACAAAGt                                                                                                                                                                                                             <  2:72422/229‑1 (MQ=255)                                            ggctatttataccTGGCGGCGGTGGCAGTCTGGTTCGTACTGGCCAATTATGTCAAAGAAGTTAAACCGGGCGTG‑CGAATAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGTAGCTGGCGCGTCCATATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCT‑CAGGCGTTTAATCAGCTACTACA‑TCCCCCGTTTCGCGCCAAGACCTGTGGCATGAACTAAGTGCGCGATGAata                                                                                                                                                                                                 <  2:321035/239‑1 (MQ=255)                                             gCTATTTTTTACTGGGGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAatat                                                                                                                                                                                                <  1:48254/251‑1 (MQ=255)                                                   ttttACTTGTGGCGGTGGCAGCCTGGTTCGTACTGGCCCTTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGt                                                                                                                                                                                           <  1:381330/250‑1 (MQ=255)                                                   ttttACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGt                                                                                                                                                                                           <  1:407230/250‑1 (MQ=255)                                                       aCTGGTGGCGGTGGCCGCCAGGTTCGTACTGGCCATTTTTGTCCAAGAAGTTAAACCGGGCGTG‑CGACGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTACCGATTTGATCTCTGGGGACCCAACGCATGGGCCACTGGCG‑CAGGCGTTTAAGCAGCTACACCA‑TCCCCCGGTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCACTGAATATCATGTCtatt                                                                                                                                                                                      >  2:381979/1‑250 (MQ=255)                                                       aCTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCCACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGt                                                                                                                                                                                                             >  1:72422/1‑228 (MQ=255)                                                          ggtggcggtggcAGCCTGGTTCGTACTGGGCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGCAGAGCAACGGAGGGGACGTTGGTCGACACCGCAACGGTGCTGGCGGAGCTGGTGCGTGCCGTTTTGCTCTCTGGGGGCCCAAAGGATGGGCAGCTGGAG‑CAGGCGTTTTATCAGCTACAACC‑TCCCCCGTTTCGCGCCAAACACGGGGGCATTTACAAAGTGCCCGATGTATATCATGTCtatttg                                                                                                                                                                                    >  2:141729/1‑247 (MQ=255)                                                                ggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATg                                                                                                                                                                             >  1:335970/1‑251 (MQ=255)                                                                      aGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAgg                                                                                                                                                                       <  1:202922/251‑1 (MQ=255)                                                                         cTGGTTCGCACTGGCCATTTTTTTCAAATAAGTTAAACCGAGCGTG‑CGAATAGCAACGGAGGGGACGTTGATCGACACCGCAACGAGGCTGGCGGCGCTGGCGCGTCCCGATTTGCTATCTGGGGACCCAACGCATGGGCAACTGACG‑CAGGCGTTTAATCAGCGACAACA‑TCCCCCGATTCGCGCCAATATCGGGGGCATTCAAAAAGTGCCCAATGAATCTCATGTcactatgcacagtacgcagggc                                                                                                                                                                     >  2:393373/1‑229 (MQ=255)                                                                           ggTTCGTACTGGCCATTTTTGTCAAAGAAGTTAATCCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACGCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAGGCGCAGGGCaa                                                                                                                                                                   <  2:71384/250‑1 (MQ=255)                                                                             ttCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGt                                                                                                                                                                >  1:24847/1‑251 (MQ=255)                                                                                  aCTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGt                                                                                                                                                           >  1:218292/1‑251 (MQ=255)                                                                                                 cAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAAGATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATa                                                                                                                                            <  2:425165/251‑1 (MQ=255)                                                                                                   aagaagGTAAACCGCGC‑TGACGAAGAGCAACGGAGGCGACGGTGAACGACACCGCAACGTTGCTGGCGGAGCTGGCGCGGCCTGATTTGCTCTCTGGGGACCCAACGGATGGGCAACTGGCT‑CAGGCGTTTAATCAGCTTCAACA‑TCCCCCGTTTCGCGCCAATATCGTTGGCATTAACAAAGTGCTCATTGAATATCATGTCTATATGACCGATGCTCAGCCCAAAGTATTGTGCGATTCGGCatattaa                                                                                                                                          <  2:39935/251‑7 (MQ=255)                                                                                                           aaaCCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCCACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTgg                                                                                                                                  >  1:319186/1‑251 (MQ=255)                                                                                                           aaaCCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTgg                                                                                                                                  <  1:272019/251‑1 (MQ=255)                                                                                                            aaCCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTgg                                                                                                                                  >  1:293136/1‑250 (MQ=255)                                                                                                             aCCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCCGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGAc                                                                                                                                >  1:300661/1‑251 (MQ=255)                                                                                                              ccgggctgg‑CGAAGAGCTACGGAGGGGACGTTGATCGTCATCGCAACGTTGCGGGCGGAGCTGGCGCGGCCCGATTGGCTCTCTTGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑GCCCCCGTTTCGCTCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAGGTCTATAGGTCCGTTGCGCAGGGCAAAATATTGTTCGATTCGGCAAATAAAGCCGTTGGACa                                                                                                                               <  2:293663/243‑1 (MQ=255)                                                                                                              ccGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGAc                                                                                                                                >  1:53204/1‑250 (MQ=255)                                                                                                              ccGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACa                                                                                                                               <  2:335970/251‑1 (MQ=255)                                                                                                                      g‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGGCCCGATTTGCTCGCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAAGATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGattat                                                                                                                        <  2:220497/250‑1 (MQ=255)                                                                                                                           agagCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTcgcg                                                                                                                   >  1:307138/1‑251 (MQ=255)                                                                                                                           agagCAACGGAGAGGACGTTAATCCACACCGCAAAGCTGCTGGCGGAGCTGGCGAGTCCCGATTTGATCTCTGGGGACCCAACGCACAGGCACCAGGCG‑CAGGCGTTTAATCAGCTACACCA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAACCGTTGGACAGGataattagcc                                                                                                                   >  2:47178/1‑247 (MQ=255)                                                                                                                              gCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTg                                                                                                                <  1:49460/251‑1 (MQ=255)                                                                                                                                         taCGTTGATCGACACCGCACCGTTGCTGGCGGAGCTGGAGCGTCCCGATTTGCCCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAAACAAGTGCGCAATGAATATCATGTCTATATGACAGATGCGCAGGGCAAAGTATTTTTCGATTCGGCAAATAAAGCCGCTGGGCAGGCTTATTCGCGCTGGAATGACGTCt                                                                                                     >  2:71911/2‑251 (MQ=255)                                                                                                                                            gTTGATCGACACCGCAACGTTGCTGGCGGCGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGc                                                                                                  >  1:133042/1‑251 (MQ=255)                                                                                                                                               gATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGGGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTa                                                                                                <  1:345099/250‑1 (MQ=255)                                                                                                                                               gATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGGCCCGTTTTGCTCTCTGGGGACCCAACGCATGGGCCACTGGAG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCACTATCGGTGGCATTAAAAAAGAGGGCAATGAATTTCATGACTATATGACCGATGAGCAGGGCAAAGTATTGTTCGATTCGGCAAATACAGCGGTTGTACAGGTTTATTCGCGCTGGAATGACGTCTGTCTac                                                                                               >  2:368946/1‑250 (MQ=255)                                                                                                                                                aTCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGCTGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAAc                                                                                              >  2:303308/1‑251 (MQ=255)                                                                                                                                                           aaCGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCACCGCATGGGCAACTGGCG‑CCGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTc                                                                                                                                                         >  1:446484/1‑181 (MQ=255)                                                                                                                                                           aaCGTTGCTGGCGGAGCTGGCGCGTCCCGAGTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCTTTTAATCAGCTACAGCA‑TCCCCCGTTTCGCGCCAATAGCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTAGTGTTc                                                                                                                                                         <  2:446484/181‑1 (MQ=255)                                                                                                                                                           aaCGTTGCTGGCGGAGCTGGCACGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGGGGCATTAACAAAGTGCGCAATGAATCTCATGTCTATCTGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGAAAGGCTTATTCGCGCTGGAATGACGTCTTGCTAAAGTTGCggtgtc                                                                                   >  2:253471/1‑246 (MQ=255)                                                                                                                                                                    ggCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTg                                                                          >  1:353583/1‑251 (MQ=255)                                                                                                                                                                        gAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCgcgc                                                                                                                  >  1:123961/1‑207 (MQ=255)                                                                                                                                                                        gAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCgcgc                                                                                                                  <  2:123961/207‑1 (MQ=255)                                                                                                                                                                          gCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCGGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCgca                                                                    <  1:197048/251‑1 (MQ=255)                                                                                                                                                                              gcgcGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCAc                                                                >  1:429100/1‑251 (MQ=255)                                                                                                                                                                                  gTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTg                                                            >  1:159407/1‑251 (MQ=255)                                                                                                                                                                                   gccctttttGCTCTCTGGGGACCCAACGCATTGGCCACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑GCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGTCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGTATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTg                                                            <  2:379692/244‑1 (MQ=255)                                                                                                                                                                                            gCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCCTTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGCTTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTg                                                  >  1:433330/1‑251 (MQ=255)                                                                                                                                                                                              tctctGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGcc                                                >  1:430501/1‑251 (MQ=255)                                                                                                                                                                                               cgctGGGTACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCGGGAATGACGTCTGGCTAACGTTGCGTGGGCAGTATGGTGCGCGCAGCACGTTGCAAAAGCCTGcc                                                <  2:429100/248‑1 (MQ=255)                                                                                                                                                                                                     ggACCCAACGCATGGGCAACTGGCG‑CAGGCTTTTAATCATCTACAACA‑TCCCCCGTTTCTCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATTTCATGTCTATATGACCGATGCGCAGGGCAAATGATTGTTCGATTCGTCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTTGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCg                                         <  1:454752/251‑1 (MQ=255)                                                                                                                                                                                                     ggACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCg                                         <  2:159407/251‑1 (MQ=255)                                                                                                                                                                                                              gCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCa                                                          <  1:175003/225‑1 (MQ=255)                                                                                                                                                                                                              gCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCa                                                          >  2:175003/1‑225 (MQ=255)                                                                                                                                                                                                                    gCAACTGGTG‑CAGGCGTTTAATCAGCTACAACA‑CCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGatgt                          <  2:137264/251‑1 (MQ=255)                                                                                                                                                                                                                          ggCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCTTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGtt                     <  1:323576/250‑1 (MQ=255)                                                                                                                                                                                                                          ggCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGatgtatg                       >  1:18888/1‑248 (MQ=255)                                                                                                                                                                                                                          ggCG‑CAGGCGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTAGTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGatgtatg                       <  2:18888/248‑1 (MQ=255)                                                                                                                                                                                                                                   cGTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCg            <  1:83579/251‑1 (MQ=255)                                                                                                                                                                                                                                    gTTTAATCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCCAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCg            >  1:279499/1‑250 (MQ=255)                                                                                                                                                                                                                                        aaTCAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTa        <  1:187980/250‑1 (MQ=255)                                                                                                                                                                                                                                           cAGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGa    <  2:353583/251‑1 (MQ=255)                                                                                                                                                                                                                                            aGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGatgtatgt                      >  1:161023/1‑232 (MQ=255)                                                                                                                                                                                                                                            aGCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGatgtatgt                      <  2:161023/232‑1 (MQ=255)                                                                                                                                                                                                                                             gCTACAACA‑TCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACg  >  1:57411/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                GCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTG‑CGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCG‑CAGGCGTTTAATCAGCTACAACA‑TCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACG  >  NC_000913/4636680‑4637164

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: ‑