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%	194.7 / NA	60	*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 (31/29);  total (31/29)

CTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGG  >  NC_000913/3870305‑3870571                                                                                                                                         |                                                                                                                                   cTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAAcc                                                                                                                                  >  1:68545/1‑139 (MQ=255)         cgcgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCt                                                                                                                          <  2:1624303/139‑1 (MQ=255)              cACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTcaca                                                                                                                     >  2:214444/1‑139 (MQ=255)              cACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTcaca                                                                                                                     >  2:478015/1‑139 (MQ=255)                   gCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACgcgc                                                                                                                           <  1:570948/128‑1 (MQ=255)                   gCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACgcgc                                                                                                                           >  2:570948/1‑128 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGTGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  1:49494/139‑1 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  1:1441229/139‑1 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        >  2:675098/1‑139 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  2:68545/139‑1 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  1:478015/139‑1 (MQ=255)                                       ttCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAagaga                                                                                            <  1:675098/139‑1 (MQ=255)                                               gATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCa                                                                                    >  2:224498/1‑139 (MQ=255)                                                           ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATaa                                                                                                <  1:1172827/115‑1 (MQ=255)                                                           ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATaa                                                                                                >  2:1172827/1‑115 (MQ=255)                                                              aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAgg                                                                                                      >  2:577301/1‑106 (MQ=255)                                                              aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAgg                                                                                                      <  1:577301/106‑1 (MQ=255)                                                              aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcg                                                                     <  1:214444/139‑1 (MQ=255)                                                                         cACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          >  2:1239588/1‑91 (MQ=255)                                                                         cACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          <  1:1239588/91‑1 (MQ=255)                                                                         cACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  1:440769/129‑1 (MQ=255)                                                                         cACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    >  2:440769/1‑129 (MQ=255)                                                                          aCCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATaa                                                                                                >  1:329201/1‑100 (MQ=255)                                                                          aCCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATaa                                                                                                <  2:329201/100‑1 (MQ=255)                                                                              gCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcgg                                                                       >  1:544335/1‑121 (MQ=255)                                                                              gCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcgg                                                                       <  2:544335/121‑1 (MQ=255)                                                                               cTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGc                                                    >  1:1518458/1‑139 (MQ=255)                                                                                                aTGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAAc                                                          <  2:295981/116‑1 (MQ=255)                                                                                                aTGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAAc                                                          >  1:295981/1‑116 (MQ=255)                                                                                                       gtggtgTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGcc                                                                                                                >  2:1481347/1‑55 (MQ=255)                                                                                                       gtggtgTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGcc                                                                                                                <  1:1481347/55‑1 (MQ=255)                                                                                                           tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAgg                        >  2:1120305/1‑139 (MQ=255)                                                                                                                cacaTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGc                                                                                                                 >  1:244265/1‑45 (MQ=255)                                                                                                                cacaTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGc                                                                                                                 <  2:244265/45‑1 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          >  2:801137/1‑50 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          <  1:801137/50‑1 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATaa                                                                                                >  1:202832/1‑60 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATaa                                                                                                <  2:202832/60‑1 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGc                                                    >  2:426401/1‑104 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGc                                                    <  1:426401/104‑1 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCcca                 >  2:876681/1‑139 (MQ=255)                                                                                                                    tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          <  1:825991/48‑1 (MQ=255)                                                                                                                    tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          >  2:825991/1‑48 (MQ=255)                                                                                                                    tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  2:1370630/86‑1 (MQ=255)                                                                                                                    tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    >  1:1370630/1‑86 (MQ=255)                                                                                                                     atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCg                                                                                                                  >  1:719057/1‑39 (MQ=255)                                                                                                                     atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCg                                                                                                                  <  2:719057/39‑1 (MQ=255)                                                                                                                     atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  2:1329728/49‑1 (MQ=255)                                                                                                                     atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        >  1:1329728/1‑49 (MQ=255)                                                                                                                     atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  2:1299996/85‑1 (MQ=255)                                                                                                                     atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    >  1:1299996/1‑85 (MQ=255)                                                                                                                       tCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCccaccacc            >  1:129557/1‑139 (MQ=255)                                                                                                                              ttGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGaa     >  2:530483/1‑139 (MQ=255)                                                                                                                               tGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAAc    >  2:1683008/1‑139 (MQ=255)                                                                                                                                gAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  2:1440441/38‑1 (MQ=255)                                                                                                                                gAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        >  1:1440441/1‑38 (MQ=255)                                                                                                                                 aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAg          >  1:673858/1‑131 (MQ=255)                                                                                                                                 aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAg          <  2:673858/131‑1 (MQ=255)                                                                                                                                 aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACgg  <  1:224498/139‑1 (MQ=255)                                                                                                                                 aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACgg  <  1:876681/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                   CTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGG  >  NC_000913/3870305‑3870571

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