Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	309,381	C→T	D216N (GAT→AAT)	ecpB ←	putative fimbrial chaperone EcpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	309,381	0	C	T	100.0%	264.6 / NA	71	D216N (GAT→AAT)	ecpB	putative fimbrial chaperone EcpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/35);  total (36/35)

GGCTACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCG  >  NC_000913/309137‑309625                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ggCTACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCt                                                                                                                                                                                                                                                 <  1:74841/250‑1 (MQ=255) ggCTACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTg                                                                                                                                                                                                                                                <  2:378263/251‑1 (MQ=255) ggCTACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAAGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTAccc                                                                                                                                                                                                                                                  <  2:474347/249‑1 (MQ=255)            agaTGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCa                                                                                                                                                                                                                                      >  1:555556/1‑250 (MQ=255)                 aaTATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGt                                                                                                                                                                                                                                <  2:391854/251‑1 (MQ=255)                        agCGGGACGCCAATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGt                                                                                                                                                                                                                         >  2:69512/1‑251 (MQ=255)                          cGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTa                                                                                                                                                                                                                       <  1:352635/251‑1 (MQ=255)                             gACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATcc                                                                                                                                                                                                                    >  2:229490/1‑251 (MQ=255)                                    aTGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTTCAATCCGTCCTTTGTTATCCGCAgtgt                                                                                                                                                                                                             <  2:28134/251‑1 (MQ=255)                                    aTGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCACACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGCACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAgtgt                                                                                                                                                                                                             >  2:497820/1‑251 (MQ=255)                                         gTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCAcg                                                                                                                                                                                                        <  1:68232/251‑1 (MQ=255)                                           ccTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCAcgcg                                                                                                                                                                                                      >  2:230667/1‑251 (MQ=255)                                            cTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCgg                                                                                                                                                                                                     >  2:456314/1‑251 (MQ=255)                                                   gcCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCCCGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAac                                                                                                                                                                                              <  2:357229/251‑1 (MQ=255)                                                        actgaAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGAcgc                                                                                                                                                                                         <  2:155954/251‑1 (MQ=255)                                                           gaAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGAcgcg                                                                                                                                                                                        <  2:62736/249‑1 (MQ=255)                                                           gaAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGAc                                                                                                                                                                                      <  1:188964/251‑1 (MQ=255)                                                           gaAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACGTATTACCCTGCCAAAGTGCAACCCGGCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGAcgcg                                                                                                                                                                                        <  2:217110/249‑1 (MQ=255)                                                                     ggAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCa                                                                                                                                                                            <  1:460819/251‑1 (MQ=255)                                                                        ataatCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCa                                                                                                                                                                         <  1:397651/251‑1 (MQ=255)                                                                                  cACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAAtt                                                                                                                                                               <  2:111485/251‑1 (MQ=255)                                                                                   cccTAGTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtc                                                                                                                                                              <  2:583727/250‑1 (MQ=255)                                                                                     cTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtctc                                                                                                                                                            >  2:144421/1‑251 (MQ=255)                                                                                                  aCTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCtt                                                                                                                                                >  2:422359/1‑250 (MQ=255)                                                                                                      aTGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGt                                                                                                                                           >  1:338646/1‑251 (MQ=255)                                                                                                      aTGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGGAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAAGGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTGTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAAGTCTCTTGACACTCCTTACCGt                                                                                                                                           <  2:252364/251‑1 (MQ=255)                                                                                                        gTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                                                                                                                                          <  2:561866/250‑1 (MQ=255)                                                                                                        gTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGGGCAACCCGGCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACGCCTTACCGtt                                                                                                                                          <  2:525877/250‑1 (MQ=255)                                                                                                            gcagcGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcgg                                                                                                                                     <  1:50653/251‑1 (MQ=255)                                                                                                                 gcTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTGAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTc                                                                                                                                >  2:299162/1‑251 (MQ=255)                                                                                                                 gcTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTc                                                                                                                                >  2:128066/1‑251 (MQ=255)                                                                                                                      cATCGGGCTGAACGCACAAAAAGACCCTGCCGAAGGCACACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGGAGTAATTCTCTTTACACTCCTTACCGTTGGCGGATTTCAGGc                                                                                                                            >  2:390103/1‑250 (MQ=255)                                                                                                                         cGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAg                                                                                                                         <  1:556232/250‑1 (MQ=255)                                                                                                                               gaacgaacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGt                                                                                                                   >  2:289933/1‑250 (MQ=255)                                                                                                                                    aacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCg                                                                                                              >  1:292837/1‑250 (MQ=255)                                                                                                                                       aaaaaGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAgg                                                                                                           >  1:32321/1‑250 (MQ=255)                                                                                                                                        aaaaGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTAGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGAt                                                                                                         <  1:282681/251‑1 (MQ=255)                                                                                                                                          aaGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATcc                                                                                                       <  2:132726/251‑1 (MQ=255)                                                                                                                                                  gCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTc                                                                                               <  1:118761/251‑1 (MQ=255)                                                                                                                                                   ccGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCg                                                                                              <  1:515279/251‑1 (MQ=255)                                                                                                                                                    cGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGc                                                                                             >  1:300503/1‑251 (MQ=255)                                                                                                                                                    cGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGc                                                                                             >  2:193716/1‑251 (MQ=255)                                                                                                                                                     gAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGcc                                                                                            >  2:11337/1‑250 (MQ=255)                                                                                                                                                            aaCTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGGAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTg                                                                                     >  1:249157/1‑251 (MQ=255)                                                                                                                                                                  gCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTg                                                                               <  1:542314/251‑1 (MQ=255)                                                                                                                                                                  gCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTg                                                                               >  2:376289/1‑251 (MQ=255)                                                                                                                                                                   cTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCCTTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGCAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCCGCGTCGCATTTCCTGTATTTGt                                                                              >  2:555531/1‑251 (MQ=255)                                                                                                                                                                     ttCAGTCCGGGGGAGAACCGTCGTAAAGGCAGTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAATTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGTATCCGCAGCGTCGCATTTCCTGTATTTGTc                                                                             <  2:529394/250‑1 (MQ=255)                                                                                                                                                                       cAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAgg                                                                          >  1:403910/1‑251 (MQ=255)                                                                                                                                                                          tCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAggg                                                                         >  1:420764/1‑249 (MQ=255)                                                                                                                                                                              ggggAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGcc                                                                     <  1:229490/249‑1 (MQ=255)                                                                                                                                                                              ggggAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCg                                                                    <  1:130364/250‑1 (MQ=255)                                                                                                                                                                               gggAGAACCGTCGTAACGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGtt                                                                  >  1:212978/1‑251 (MQ=255)                                                                                                                                                                                                  cATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGt                                                 >  1:207911/1‑249 (MQ=255)                                                                                                                                                                                                      gTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCg                                           >  2:288537/1‑251 (MQ=255)                                                                                                                                                                                                       tCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGc                                          >  1:392311/1‑251 (MQ=255)                                                                                                                                                                                                        cGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGcc                                         <  1:127470/251‑1 (MQ=255)                                                                                                                                                                                                          ttATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTg                                       >  1:462134/1‑251 (MQ=255)                                                                                                                                                                                                                  cTGCACAGCTATCTATTTCACGGGACTGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGAAGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTTAAAGTGGTAGTTAGCCTGACGAgggg                               >  2:251416/1‑251 (MQ=255)                                                                                                                                                                                                                    gCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAgggg                               <  2:238542/249‑1 (MQ=255)                                                                                                                                                                                                                    gCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAgggg                               <  1:468398/249‑1 (MQ=255)                                                                                                                                                                                                                     cacaGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGcg                             >  1:449765/1‑250 (MQ=255)                                                                                                                                                                                                                      acaGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGcgac                           <  1:208146/251‑1 (MQ=255)                                                                                                                                                                                                                           tatctatTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGAcc                       >  2:436351/1‑250 (MQ=255)                                                                                                                                                                                                                                atTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTGGCCCGGCATCAGGTAGTAATTCTCGTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAAt                  <  2:520315/250‑1 (MQ=255)                                                                                                                                                                                                                                    cACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGGAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGc             >  2:272990/1‑251 (MQ=255)                                                                                                                                                                                                                                    cACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCAGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGc             >  1:374415/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaCGGGCATGAACTTATTACCCTGCCAAAGTGCAACCAGTCCTTTGTTATCCGAAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCCTCAGGTAGCAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCACGGTCCGTAGGCGAGGATCAGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGc             >  1:374412/2‑251 (MQ=255)                                                                                                                                                                                                                                     aCGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGcc            <  2:207911/251‑1 (MQ=255)                                                                                                                                                                                                                                          aaTGAACTTATTACACTCCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATgc       >  1:476158/1‑251 (MQ=255)                                                                                                                                                                                                                                               aCTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATgcgggcg  >  2:574949/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GGCTACCGGCGAGATGAATATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCG  >  NC_000913/309137‑309625

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