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%	160.8 / NA	50	*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 (23/27);  total (23/27)

TTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATG  >  NC_000913/3870306‑3870576                                                                                                                                        |                                                                                                                                        ttACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCa                                                                                                                                      <  2:397528/139‑1 (MQ=255)   aCTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAc                                                                                                                                    >  1:584302/1‑139 (MQ=255)                                                tAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCcca                                                                                                                   >  2:304169/1‑111 (MQ=255)                                                tAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCcca                                                                                                                   <  1:304169/111‑1 (MQ=255)                                                   tAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                               >  2:982049/1‑112 (MQ=255)                                                   tAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                               <  1:982049/112‑1 (MQ=255)                                                          ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCt                                                                                                                               <  1:1099590/89‑1 (MQ=255)                                                          ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCt                                                                                                                               >  2:1099590/1‑89 (MQ=255)                                                             aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCccac                                                                                                                  >  2:554431/1‑99 (MQ=255)                                                             aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCccac                                                                                                                  <  1:554431/99‑1 (MQ=255)                                                             aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcg                                                                          >  2:600937/1‑139 (MQ=255)                                                              tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                         <  2:584302/139‑1 (MQ=255)                                                              tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                         <  2:960696/139‑1 (MQ=255)                                                              tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                         <  1:745248/139‑1 (MQ=255)                                                              tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                         <  2:881967/139‑1 (MQ=255)                                                                    aCCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                               >  1:1190993/1‑95 (MQ=255)                                                                    aCCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                               <  2:1190993/95‑1 (MQ=255)                                                                        cACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAAc                                                               >  1:684225/1‑139 (MQ=255)                                                                         aCCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAAc                                                                                   >  1:186116/1‑118 (MQ=255)                                                                         aCCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAAc                                                                                   <  2:186116/118‑1 (MQ=255)                                                                               tttAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCa                                                        <  2:492068/139‑1 (MQ=255)                                                                               tttAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCa                                                        <  2:63522/139‑1 (MQ=255)                                                                                  aaGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGt                                                     >  2:219472/1‑139 (MQ=255)                                                                                   aGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAcg                                                                                                                                   <  1:780453/60‑1 (MQ=255)                                                                                   aGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAcg                                                                                                                                   >  2:780453/1‑60 (MQ=255)                                                                                                tGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAAcc                                       <  1:802901/139‑1 (MQ=255)                                                                                                tGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAAcc                                       <  1:600937/139‑1 (MQ=255)                                                                                                tGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAAcc                                       <  2:684225/139‑1 (MQ=255)                                                                                                  cgcAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCt                                     >  2:923006/1‑139 (MQ=255)                                                                                                      gtggtgTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcgcg                                                                                            >  1:322329/1‑80 (MQ=255)                                                                                                      gtggtgTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcgcg                                                                                            <  2:322329/80‑1 (MQ=255)                                                                                                          tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                               >  1:474825/1‑57 (MQ=255)                                                                                                          tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                               <  2:474825/57‑1 (MQ=255)                                                                                                                 caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGTGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCcca                      >  1:1084205/1‑139 (MQ=255)                                                                                                                 caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAgg                             <  2:670577/132‑1 (MQ=255)                                                                                                                 caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAgg                             >  1:670577/1‑132 (MQ=255)                                                                                                                   tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcg                                                                                              >  1:419827/1‑65 (MQ=255)                                                                                                                   tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcg                                                                                              <  2:953666/65‑1 (MQ=255)                                                                                                                   tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcg                                                                                              <  2:419827/65‑1 (MQ=255)                                                                                                                   tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcg                                                                                              >  1:953666/1‑65 (MQ=255)                                                                                                                    atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCccacca                   >  2:310285/1‑139 (MQ=255)                                                                                                                        cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                             <  2:276711/45‑1 (MQ=255)                                                                                                                        cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                             >  1:276711/1‑45 (MQ=255)                                                                                                                          caGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGt                            <  2:433872/124‑1 (MQ=255)                                                                                                                          caGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGt                            >  1:433872/1‑124 (MQ=255)                                                                                                                           aGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCgg            >  2:12767/1‑139 (MQ=255)                                                                                                                                aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACgg       <  1:12767/139‑1 (MQ=255)                                                                                                                                  ggATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                         <  1:827421/71‑1 (MQ=255)                                                                                                                                  ggAATAACCAACGCGCTTCAAATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                         >  2:827421/1‑71 (MQ=255)                                                                                                                                     ttAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATg  <  1:219472/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        TTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATG  >  NC_000913/3870306‑3870576

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