Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,377,330	G→A	C131C (TGC→TGT)	blc ←	outer membrane lipoprotein Blc

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,377,330	0	G	A	100.0%	64.1 / NA	20	C131C (TGC→TGT)	blc	outer membrane lipoprotein Blc
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/6);  total (14/6)

GCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCGCAAA‑CCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGC  >  NC_000913/4377200‑4377462                                                                                                                                   |                                                                                                                                      gctgctGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAgc                                                                                                                             >  2:154156/1‑140 (MQ=255)   tgctgTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAgcgc                                                                                                                           >  1:38791/1‑140 (MQ=255)     ctgTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCATTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCAt                                                                                                                         <  1:333271/140‑1 (MQ=255)     ctgTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGGGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCAt                                                                                                                         >  2:247296/1‑140 (MQ=255)     ctgTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCGGAAA‑CCAGCGCac                                                                                                                         >  2:625149/1‑139 (MQ=255)     ctgTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCAt                                                                                                                         <  2:496766/140‑1 (MQ=255)     ctgTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCAt                                                                                                                         <  1:588349/140‑1 (MQ=255)         aCCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCg                                                                                                                     <  2:38791/140‑1 (MQ=255)               ataaatTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGAATTc                                                                                                               >  2:468841/1‑140 (MQ=255)               ataaatTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGCCGCGGTCCGGGTCAAAAA‑CCACCGCATGGCGGTATTc                                                                                                               >  1:320398/1‑140 (MQ=255)                     ttACTGACATCAAGCCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTCCCGATc                                                                                                         <  1:597060/140‑1 (MQ=255)                           aCATCAAACCCTTCCCGGGGCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGAAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTTCCGATCGAGTGc                                                                                                   >  1:630614/1‑140 (MQ=255)                             aTCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTCCCGATCGAGTGCaa                                                                                                 <  1:151001/140‑1 (MQ=255)                             aTCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTCCCGATCGAGTGCaa                                                                                                 >  2:475372/1‑140 (MQ=255)                                       tCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGGATACCAGATCGACTGCAATAACGTtata                                                                                       >  1:308317/1‑140 (MQ=255)                                       tCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGG‑‑GAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTGTTCACGATCGAGTGCAATAACGTTATAAc                                                                                     >  1:211046/1‑140 (MQ=255)                                                   aCTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGTATGGCGGTCTTCCCGACCGAGTGCAATAACGTTATAACCGCCATAGaa                                                                           >  2:263256/1‑140 (MQ=255)                                                       ccAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGa                                                                       >  2:536545/1‑140 (MQ=255)                                                       ccAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTCCCGACCGACTGCAATAACGTTACAACCGCCATAGAAAGGa                                                                       >  2:296660/1‑140 (MQ=255)                                                                         cTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTc                                                     >  2:390277/1‑140 (MQ=255)                                                                              tCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTcagcgc                                                >  2:465892/1‑140 (MQ=255)                                                                                                 gggAGAGTACCCAAAGGTAGCCGCGGCCGGGCCCACAAC‑CCAGGGCAGGGGGTTTTCCCCGATCGAGTGCAAAAACTTTAAAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCg                             <  2:27145/140‑1 (MQ=255)                                                                                                                    gTCGCGGTCCGGGCCACAAA‑CCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCtt          >  1:654608/1‑140 (MQ=255)                                                                                                                             cGGGCCCAAAACCCAGCGCAGGGCGTTATCCCGGTGCGGGTGCAAAAGTTTAATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTctgc  <  1:475372/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                      GCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCGCAAA‑CCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGC  >  NC_000913/4377200‑4377462

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