Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	154,182	C→T	L415L (CTG→CTA)	htrE ←	putative fimbrial usher protein HtrE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	154,182	0	C	T	100.0%	75.9 / NA	23	L415L (CTG→CTA)	htrE	putative fimbrial usher protein HtrE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/11);  total (12/11)

GATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTCAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCATCGG  >  NC_000913/153968‑154417                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            gATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGttatta                                                                                                                                                                                                         <  1:324993/251‑1 (MQ=255)     tttCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCCAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGCATAGTTATTATCgg                                                                                                                                                                                                     >  2:314703/1‑251 (MQ=255)      ttCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGt                                                                                                                                                                                                    >  1:401216/1‑251 (MQ=255)        cACTTCAGCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCGTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCAGCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTa                                                                                                                                                                                                  <  2:253708/251‑1 (MQ=255)                  aTTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATAcc                                                                                                                                                                                        <  2:412175/251‑1 (MQ=255)                                                tAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGCAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGcc                                                                                                                                                          >  1:365096/1‑251 (MQ=255)                                                 aaTTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGcc                                                                                                                                                          <  2:139052/250‑1 (MQ=255)                                                                               tatTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGc                                                                                                                           >  1:171295/1‑251 (MQ=255)                                                                                atTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGc                                                                                                                           <  2:401216/250‑1 (MQ=255)                                                                                       gaagttgttGTTTCTTAGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTAGCATCCGGGATACGAACATGGGAATGAGTCAGATCGAAAGAAAATGCTCCAACAGAAGTATTTAGACAAAGACCTAAAAAAACAGCGGGATGGTTATTATCGGTTAGCTGAATACCTGTATAACCCGGCATTTAGTTATTCAGGCCGTAGTAGTAGATTGCTTGAAATAAATTAGGCGCACCCTg                                                                                                                   <  2:174007/244‑1 (MQ=255)                                                                                         aGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGTCCTAACAAACCAGCGGTAGAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGaa                                                                                                                 <  1:547746/251‑1 (MQ=255)                                                                                                             aCTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTg                                                                                              >  2:330050/1‑250 (MQ=255)                                                                                                                tggtCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCGTTTAAGACCTGAc                                                                                            <  2:436253/246‑1 (MQ=255)                                                                                                                      aGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCt                                                                                     >  2:135921/1‑250 (MQ=255)                                                                                                                                                   gTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGAATTTAGACCAAGACCTACCAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTa                                                        >  2:206231/1‑250 (MQ=255)                                                                                                                                                                                   aGTCACATCGAAAGAAAATGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCg                          >  2:441674/1‑248 (MQ=255)                                                                                                                                                                                           agaaagaaaATGACCCAACTGAATTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGTTTATCTGAATACCGGTATAACCCGTCAGATATTTATTCAGGCCGTAGTAGGAGCTGGCTTGAAATAAATTAGGCTCATCCGGAATATCATCTTTAAAGCCCGGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTACCATTTTAACAACGGATGAGAAAGGTTGCGAGAATGTCCGc               <  1:175789/250‑1 (MQ=255)                                                                                                                                                                                              aagaaaAGGCCCCAACGGAAGTATTTAGACCAAGACCTAACAACCCAGCGTTATAGTTATTATCGGTTATCTGAATCCCGTTATAACCCGTAAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTGGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCttt            <  1:408390/251‑1 (MQ=255)                                                                                                                                                                                              aaaaaaaTGCACAAACTGAAGTATTTAGCCCAAGCCCAAACAAACCAGCGTTATATTTATTAGCGGTTATATGAATACCGGTATATCCCGTCAGATATTTATTCAGGCCGTAGTAGTAGCTGGCATAAAATAAATTAGGCTCACCCTAAATACCTTCTTTAAAGCCCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCttt            <  2:261942/248‑1 (MQ=255)                                                                                                                                                                                                  aaaTGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGc        >  1:136637/1‑251 (MQ=255)                                                                                                                                                                                                  aaaTGCACCAACTGAAGGATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGACCGCTTTGAGc        >  2:448231/1‑251 (MQ=255)                                                                                                                                                                                                     tGCACCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCAt     >  1:32929/1‑251 (MQ=255)                                                                                                                                                                                                        aCCAACTGAAGTATTTAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCATCgg  >  1:14265/1‑251 (MQ=255)                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            GATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATGTTTTATCATCCGGGATACGAACATTGGAATGAGTCACATCGAAAGAAAATGCACCAACTGAAGTATTCAGACCAAGACCTAACAAACCAGCGGTATAGTTATTATCGGTTATCTGAATACCGGTATAACCCGTCAGATAGTTATTCAGGCCGTAGTAGTAGCTTGCTTGAAATAAATTAGGCTCATCCTGAATATCATCTTTTAAGACCTGACCGCCGCTAATATCCCAACGTCCAACGCCAGGGCGTAACATTTGAACAACGGATGAGAAAGGTTGCGAGAATGTCCGCTTTGAGCCATCGG  >  NC_000913/153968‑154417

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