Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	156,068	A→G	V45A (GTA→GCA)	yadV ←	putative fimbrial chaperone YadV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	156,068	0	A	G	100.0%	70.7 / NA	26	V45A (GTA→GCA)	yadV	putative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/18);  total (8/18)

GCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTAT  >  NC_000913/155828‑156268                                                                                                                                                                                                                                                 |                                                                                                                                                                                                         gCATCTGGTTTTGGTGGAACTTCCATTGCGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCGCAGCGTTGTCATCGCCAGTATCTAACCATCTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGaca                                                                                                                                                                                                  <  2:242122/249‑1 (MQ=255)            tggttgAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTCCCAGCCGTGCGTTGACACTGTGTTGATCg                                                                                                                                                                                      <  2:461736/250‑1 (MQ=255)                ggAACTTTCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTATGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTTTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCTTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCttt                                                                                                                                                                                  <  2:544863/250‑1 (MQ=255)                     ttCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGATTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTtatat                                                                                                                                                                             >  1:73515/1‑250 (MQ=255)                                                          ttAGGCAGTGATGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTAGCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAGTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCaac                                                                                                                                          <  2:60872/248‑2 (MQ=255)                                                                aGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATg                                                                                                                                  >  1:8617/1‑250 (MQ=255)                                                                 gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGa                                                                                                                                 >  1:86279/1‑250 (MQ=255)                                                                               gtgtACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACa                                                                                                                  <  2:183739/251‑1 (MQ=255)                                                                                 gtACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAg                                                                                                                 <  1:386855/250‑1 (MQ=255)                                                                                                   ttGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGTTGAGTCCCTTTATTTGCCAGACGTGCGTTGACACttttt                                                                                                                                                                                            <  2:314766/157‑1 (MQ=255)                                                                                                   ttGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTCTTCCCTTTATTTTCCAGACGTGCGTTGCCACttttt                                                                                                                                                                                            >  1:314766/1‑157 (MQ=255)                                                                                                             ttGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGtt                                                                                     <  1:311368/250‑1 (MQ=255)                                                                                                                               cctgtccgcttatcAGTAAAAAGGACTGTATTACTGCCATGCTCAGCGTTGTCATCTCCATCATCTAAACAACTCTGGACTATCTACGGGTTATTCCCTTTATTTTCCAGACTTGCGTTGACCCTTTTTTGATCTCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGTTGATGAACTAAAATCCATACAGGTTACGCAGCATAAATCCGTTGTATGTTCGGCGTTTACAAACATata                                                                  <  2:95583/238‑3 (MQ=255)                                                                                                                                  ggcggcGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGCCAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGCGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCAATACAGGTTACGAAGCATCAAGCTGTTGTCTGTTTGGTGTTAAAAACCATAAAcc                                                                >  2:454054/1‑250 (MQ=255)                                                                                                                                  ggcggcGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAAcc                                                                >  2:544915/1‑250 (MQ=255)                                                                                                                                      gcgTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACGTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACTATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATccc                                                           <  2:98544/251‑1 (MQ=255)                                                                                                                                       cgTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCt                                                          <  1:544915/251‑1 (MQ=255)                                                                                                                                             aGTAAAAGGTACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACGCTGGACAAGCAACGGGTTATTCCCTTTATTTACCAGACGTGCGTGGACCCTTTTTTGATCGCGTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCGGTTGTATGTCTGGTGTTTAAAAACATAAACCATCCCTGATaga                                                    <  2:191479/251‑1 (MQ=255)                                                                                                                                             aGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAg                                                                                                                                        >  1:109899/1‑167 (MQ=255)                                                                                                                                             aGTAAAAGGGACTGTAATACTGCCAGGCGCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTGATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAg                                                                                                                                        <  2:109899/167‑1 (MQ=255)                                                                                                                                                                          cAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATaaaa                        <  1:206679/250‑1 (MQ=255)                                                                                                                                                                            gagttctcctccccAGTTTCTAACCAACCGTGGACAAGCTTCGGGTTATGTCCTTTATTTTCCAGACGTGCGGTGACCCTTTTTTGAACGTTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGATCTAAAAGCCATACAGGTTACGAATCATAAAGCTGTTGTATGCTTGGGGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATaaac                        <  2:427310/236‑2 (MQ=255)                                                                                                                                                                            gCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCGGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAg                     <  1:371680/251‑1 (MQ=255)                                                                                                                                                                              gTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGa                    >  1:33453/1‑250 (MQ=255)                                                                                                                                                                                      cGCCAGTATCTAACCAAGTCTGGACAAGCAACGGGTTATTCTCTTTAATTTCCAGACGTGCGTTGACACTTTTTTGATCGCTGGAATATATTACGCGAGTACGCGATATGACAATGGCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAGAGCTGTTGTATGTTTGGGGTTAAAAAACATAAATAATCCCTGATAGATGGAAAGAAAACAGGGAAGAATATAAAATAGAAattatta            <  2:86279/250‑1 (MQ=255)                                                                                                                                                                                                 aaCCAACTCTGTACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACCCTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGATAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGtat  <  1:318418/249‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                         GCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTAT  >  NC_000913/155828‑156268

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