Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,363,151	T→C	G115G (GGT→GGC)	yhcD →	putative outer membrane fimbrial subunit usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,363,151	0	T	C	100.0%	182.4 / NA	53	G115G (GGT→GGC)	yhcD	putative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/17);  total (36/17)

AATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAACTACTCTCTGATGAGC  >  NC_000913/3363005‑3363289                                                                                                                                                   |                                                                                                                                           aaTGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCAcc                                                                                                                                         >  2:770225/1‑150 (MQ=255) aaTGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACca                                                                                                                                        <  2:71020/151‑1 (MQ=255)      cGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATc                                                                                                                                   >  2:522621/1‑151 (MQ=255)       gAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATc                                                                                                                                   >  2:694090/1‑150 (MQ=255)       gAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATc                                                                                                                                   >  2:64477/1‑150 (MQ=255)           aaaaaGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACAcc                                                                                                                              >  2:493710/1‑151 (MQ=255)              aaGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCAtt                                                                                                                           >  1:69710/1‑151 (MQ=255)                  cAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCt                                                                                                                       >  2:769495/1‑151 (MQ=255)                   aaCATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTa                                                                                                                      <  2:190592/151‑1 (MQ=255)                    aCATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTaa                                                                                                                     <  1:549335/151‑1 (MQ=255)                    aCATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTa                                                                                                                      >  1:515225/1‑150 (MQ=255)                            aCGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGa                                                                                                               >  1:458184/1‑149 (MQ=255)                              gcgcTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCg                                                                                                           <  2:268587/151‑1 (MQ=255)                                 cTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGtt                                                                                                         <  2:661508/150‑1 (MQ=255)                                   ttGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCCCACCATTACCTAACCAGGAAGCGTTAg                                                                                                       <  2:295194/150‑1 (MQ=255)                                    tgatgaAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGa                                                                                                      <  2:11161/150‑1 (MQ=255)                                      atgaAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGatat                                                                                                   >  1:727481/1‑151 (MQ=255)                                      atgaAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGatat                                                                                                   >  1:62777/1‑151 (MQ=255)                                         aacaTGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGAtattat                                                                                                >  1:744802/4‑151 (MQ=255)                                               gCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTc                                                                                          >  2:625238/1‑151 (MQ=255)                                                 gACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTcac                                                                                        >  1:562243/1‑151 (MQ=255)                                                        gTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGc                                                                                 <  2:187340/151‑1 (MQ=255)                                                         tCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCg                                                                                >  2:752788/1‑151 (MQ=255)                                                         tCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCg                                                                                >  2:326932/1‑151 (MQ=255)                                                           ttGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGAt                                                                              >  1:808493/1‑151 (MQ=255)                                                           ttGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGAt                                                                              >  1:46100/1‑151 (MQ=255)                                                              aTCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCAt                                                                           >  2:735626/1‑151 (MQ=255)                                                               tCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCAtt                                                                          >  2:683072/1‑151 (MQ=255)                                                                 aGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTcc                                                                        <  1:194720/151‑1 (MQ=255)                                                                        tCTGCAACAAGCTGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAgg                                                                  <  2:385997/150‑1 (MQ=255)                                                                        tCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAggg                                                                 >  1:437616/1‑151 (MQ=255)                                                                         cTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAggg                                                                 >  2:220948/1‑150 (MQ=255)                                                                             aacaaGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTgg                                                             >  1:260044/1‑150 (MQ=255)                                                                              acaaGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAt                                                           >  1:748361/1‑151 (MQ=255)                                                                                 aGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtctc                                                        >  1:393468/1‑151 (MQ=255)                                                                                   ccGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCACATAGGGTTGGATCTca                                                       >  1:214092/1‑150 (MQ=255)                                                                                    cGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTcac                                                      >  1:838687/1‑150 (MQ=255)                                                                                           aaaaTTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACgc                                                >  1:139132/1‑149 (MQ=255)                                                                                                 ccTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtg                                        >  2:811808/1‑151 (MQ=255)                                                                                                 ccTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACt                                         <  1:522583/150‑1 (MQ=255)                                                                                                    tCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtggt                                      <  2:30057/150‑1 (MQ=255)                                                                                                          gaaAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAAc                                 >  1:816933/1‑149 (MQ=255)                                                                                                             aaaaaTGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAgctg                            >  1:237873/1‑151 (MQ=255)                                                                                                              aaaaTGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAgctgc                           >  2:255520/1‑151 (MQ=255)                                                                                                               aaaTGGCTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAgctgcg                          <  1:617020/151‑1 (MQ=255)                                                                                                                     cTGTTACGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCgctgct                     >  2:822310/1‑150 (MQ=255)                                                                                                                          aCGACTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAacta               <  1:694047/151‑1 (MQ=255)                                                                                                                             aCTATATATTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAACTActc            >  2:465606/1‑151 (MQ=255)                                                                                                                               tatataTTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAACTActctc          >  2:688041/1‑151 (MQ=255)                                                                                                                                atatatTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAACTActctct         <  2:191236/151‑1 (MQ=255)                                                                                                                                  atatTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAACTACTCTCtg        >  2:748538/1‑150 (MQ=255)                                                                                                                                    atTGTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAACTACTCTCtgat      <  2:409391/150‑1 (MQ=255)                                                                                                                                       gTCCTACCCGGGCACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAACTACTCTCTGATGagc  <  2:243811/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                           AATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTGCGCTGCTAAACTACTCTCTGATGAGC  >  NC_000913/3363005‑3363289

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