Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	752,258	T→G	N180H (AAC→CAC)	ybgD ←	putative fimbrial protein YbgD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	752,258	0	T	G	100.0%	90.4 / NA	28	N180H (AAC→CAC)	ybgD	putative fimbrial protein YbgD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/10);  total (18/10)

GGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATCTGGCGCGGCTGAACCT  >  NC_000913/752140‑752379                                                                                                                       |                                                                                                                          ggATAC‑AATGAGAGTCGACAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCAcc                                                                                                            <  1:220531/133‑1 (MQ=255)      cAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATCCTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCCCCTGcc                                                                                                        <  2:73907/133‑1 (MQ=255)            agagTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACa                                                                                                  >  1:619090/1‑133 (MQ=255)            agagTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACa                                                                                                  >  1:388768/1‑133 (MQ=255)             gagTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTACGTGAGCCGTTTCTTCACCTGGCGTGACAt                                                                                                 >  1:304546/1‑133 (MQ=255)                tCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGa                                                                                              >  2:396207/1‑133 (MQ=255)                          ttCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTACACCTGCCGTGACATCGACTGCATTATc                                                                                    >  2:140234/1‑133 (MQ=255)                                  gtCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGtt                                                                            >  1:326240/1‑133 (MQ=255)                                       tACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTtccatc                                                                       >  2:332992/1‑133 (MQ=255)                                           tGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAgg                                                                   >  1:345450/1‑133 (MQ=255)                                             aTAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCg                                                                 >  2:16446/1‑133 (MQ=255)                                             aTAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCg                                                                 >  2:564074/1‑133 (MQ=255)                                             aTAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCg                                                                 >  1:362631/1‑133 (MQ=255)                                             aTAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCg                                                                 >  1:521924/1‑133 (MQ=255)                                               aaaTATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGGGGTTACTTCACCTGCCGTGACATCGACTGCAATATCAATTTGTTCCATCCAGGCGaa                                                               <  1:53205/133‑1 (MQ=255)                                                      ttATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTc                                                        <  2:521924/133‑1 (MQ=255)                                                           aaaGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGt                                                   <  1:364196/133‑1 (MQ=255)                                                           aaaGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGt                                                   <  2:362631/133‑1 (MQ=255)                                                             aGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCt                                                 <  1:564074/133‑1 (MQ=255)                                                                      tgatgCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCtgatg                                        <  2:345450/133‑1 (MQ=255)                                                                      tgatgCCCCCAAAATAATTCTTATTTATAATCCAGCACGTAGGTTGCGTGCGCGTTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCtgatg                                        <  2:313029/133‑1 (MQ=255)                                                                                aaaaTAATACTTATTTATAATCCAGCACGTAGGTTGCGTGAGCGTTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATc                              <  1:396207/133‑1 (MQ=255)                                                                                            atttatAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATc                  >  2:364673/1‑133 (MQ=255)                                                                                            atttatAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATc                  >  1:325935/1‑133 (MQ=255)                                                                                              ttatAATCCAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATCTg                >  2:484136/1‑133 (MQ=255)                                                                                                      cAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATCTGGTGCGGCt        >  2:443027/1‑133 (MQ=255)                                                                                                      cAGCACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATCTGGCGCGGCt        <  1:364673/133‑1 (MQ=255)                                                                                                         cACGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAAGTCAAGATCTGGCTCGGCTGaa     >  1:246968/1‑133 (MQ=255)                                                                                                          aCGTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATCTGGCGCGGCTGAAc    <  2:325935/133‑1 (MQ=255)                                                                                                            gTAGGTTGCGTGAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATCTGGCGCGGCTGAACCt  >  1:374182/1‑133 (MQ=255)                                                                                                                       |                                                                                                                          GGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGTTCTGATGAGCTTGCATCCAGGTCAAGATCTGGCGCGGCTGAACCT  >  NC_000913/752140‑752379

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