Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,870,441	G→A	*355* (TAG→TAA)  G430G (GGC→GGT)	cbrA → dgoT ←	colicin M resistance protein, FAD‑binding protein, putative oxidoreductase D‑galactonate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,870,441	0	G	A	100.0%	166.2 / NA	51	*355* (TAG→TAA)  G430G (GGC→GGT)	cbrA dgoT	colicin M resistance protein, FAD‑binding protein, putative oxidoreductase D‑galactonate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/29);  total (22/29)

TACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAA  >  NC_000913/3870307‑3870568                                                                                                                                       |                                                                                                                                tACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCaa                                                                                                                             <  2:816669/139‑1 (MQ=255)      gcgcgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACgcgc                                                                                                                        <  2:1256807/139‑1 (MQ=255)       cgcgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCt                                                                                                                       >  2:1122617/1‑139 (MQ=255)       cgcgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCt                                                                                                                       >  1:800681/1‑139 (MQ=255)              cccGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATc                                                                                                                <  1:95827/139‑1 (MQ=255)                       tGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  2:217439/139‑1 (MQ=255)                         cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  1:1104243/139‑1 (MQ=255)                         cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  2:800681/139‑1 (MQ=255)                         cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  1:1122617/139‑1 (MQ=255)                         cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  1:1251963/139‑1 (MQ=255)                                    ctTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCGCATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAagag                                                                                          >  1:1299363/1‑139 (MQ=255)                                                     aaaGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAAcgc                                                                         <  1:88139/139‑1 (MQ=255)                                                         ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcg                                                                     >  1:703666/1‑139 (MQ=255)                                                            aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcg                                                                  >  1:404325/1‑139 (MQ=255)                                                             tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 <  2:977119/139‑1 (MQ=255)                                                             tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 <  2:240339/139‑1 (MQ=255)                                                              gCCTGACCGCACCGGCTTT‑‑‑‑AAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                            <  1:366575/139‑1 (MQ=255)                                                                  gACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                            <  2:703666/139‑1 (MQ=255)                                                                                  aGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTg                                            >  2:297209/1‑139 (MQ=255)                                                                                     aaGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         <  1:284069/139‑1 (MQ=255)                                                                                       gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCa                                                                                                                              >  2:349159/1‑52 (MQ=255)                                                                                       gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCa                                                                                                                              <  1:349159/52‑1 (MQ=255)                                                                                       gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGaaa                                       >  1:1158086/1‑139 (MQ=255)                                                                                          gATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGa                                                                       <  1:1288822/104‑1 (MQ=255)                                                                                          gATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGa                                                                       >  2:1288822/1‑104 (MQ=255)                                                                                               tGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAAcc                               <  2:829276/139‑1 (MQ=255)                                                                                                         tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAc                                                                      <  1:1177838/90‑1 (MQ=255)                                                                                                         tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAc                                                                      >  2:1177838/1‑90 (MQ=255)                                                                                                           tgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAagaga                                                                                         <  2:356341/69‑1 (MQ=255)                                                                                                           tgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAagaga                                                                                         >  1:356341/1‑69 (MQ=255)                                                                                                                caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  2:1155365/50‑1 (MQ=255)                                                                                                                caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       >  1:1155365/1‑50 (MQ=255)                                                                                                                     tCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         <  2:206502/107‑1 (MQ=255)                                                                                                                     tCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         >  1:206502/1‑107 (MQ=255)                                                                                                                       cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         <  2:1148581/105‑1 (MQ=255)                                                                                                                       cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         >  1:1148581/1‑105 (MQ=255)                                                                                                                       cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAg       <  2:1299363/139‑1 (MQ=255)                                                                                                                        acaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCg                                    <  1:658458/109‑1 (MQ=255)                                                                                                                        acaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCg                                    >  2:658458/1‑109 (MQ=255)                                                                                                                            ttGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGaa  >  1:394161/1‑139 (MQ=255)                                                                                                                            ttGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGaa  >  2:624354/1‑139 (MQ=255)                                                                                                                            ttGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGaa  >  1:715587/1‑139 (MQ=255)                                                                                                                             tGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 >  1:1091248/1‑75 (MQ=255)                                                                                                                             tGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 <  2:1091248/75‑1 (MQ=255)                                                                                                                               aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcgg                                                                    >  1:904477/1‑70 (MQ=255)                                                                                                                               aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcgg                                                                    <  2:904477/70‑1 (MQ=255)                                                                                                                               aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAgg                     <  2:1299461/117‑1 (MQ=255)                                                                                                                               aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAgg                     >  1:1299461/1‑117 (MQ=255)                                                                                                                                    ttAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 <  1:728858/68‑1 (MQ=255)                                                                                                                                    ttAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAg                                               <  2:975319/86‑1 (MQ=255)                                                                                                                                    ttAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAg                                               >  1:975319/1‑86 (MQ=255)                                                                                                                                       |                                                                                                                                TACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAA  >  NC_000913/3870307‑3870568

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