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%	208.4 / NA	64	*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 (32/32);  total (32/32)

CTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAA  >  NC_000913/3870305‑3870568                                                                                                                                         |                                                                                                                                cTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAAcc                                                                                                                               >  2:920324/1‑139 (MQ=255) cTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAAcc                                                                                                                               >  1:929384/1‑139 (MQ=255) cTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAAcc                                                                                                                               >  1:607149/1‑139 (MQ=255)  ttACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCa                                                                                                                              <  1:752178/139‑1 (MQ=255)              cACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTcaca                                                                                                                  <  2:892445/139‑1 (MQ=255)                         tGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  2:929384/139‑1 (MQ=255)                         tGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  2:607149/139‑1 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGTACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       >  2:224492/1‑137 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGTACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  1:224492/137‑1 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAc                                                                                                                            <  2:670808/116‑1 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAc                                                                                                                            >  1:670808/1‑116 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     >  1:145200/1‑139 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     >  2:307345/1‑139 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  1:920324/139‑1 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  1:307345/139‑1 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     >  1:628115/1‑139 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     >  1:444627/1‑139 (MQ=255)                           cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     >  1:553715/1‑139 (MQ=255)                             tttAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAgg                                                                                                   <  2:145200/139‑1 (MQ=255)                                   aCTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATaa                                                                                             <  2:444627/139‑1 (MQ=255)                                                    tAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCaa                                                                            >  1:882952/1‑139 (MQ=255)                                                         aGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGa                                                                       <  1:1134414/139‑1 (MQ=255)                                                         aGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGa                                                                       >  2:1134414/1‑139 (MQ=255)                                                          gCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAc                                                                      >  1:472528/1‑139 (MQ=255)                                                              aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     >  2:138811/1‑104 (MQ=255)                                                              aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  1:138811/104‑1 (MQ=255)                                                              aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcg                                                                  >  2:1787448/1‑139 (MQ=255)                                                              aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcg                                                                  >  2:1716651/1‑139 (MQ=255)                                                              aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcg                                                                  >  1:612741/1‑139 (MQ=255)                                                               tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 <  2:553715/139‑1 (MQ=255)                                                                    gACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                            <  1:355822/139‑1 (MQ=255)                                                                         cACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAAc                                                       <  1:1787448/139‑1 (MQ=255)                                                                                       aaGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         <  2:612741/139‑1 (MQ=255)                                                                                       aaGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         <  1:1716651/139‑1 (MQ=255)                                                                                         gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  1:250826/75‑1 (MQ=255)                                                                                         gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       >  2:250826/1‑75 (MQ=255)                                                                                         gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAgg                                                                                                   <  1:1123780/79‑1 (MQ=255)                                                                                         gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAgg                                                                                                   >  2:1123780/1‑79 (MQ=255)                                                                                         gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGgag                                                               <  2:1229641/115‑1 (MQ=255)                                                                                         gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGgag                                                               >  1:1229641/1‑115 (MQ=255)                                                                                                 tGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 <  1:221636/105‑1 (MQ=255)                                                                                                 tGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 >  2:221636/1‑105 (MQ=255)                                                                                                      aGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTgcg                          <  2:882952/139‑1 (MQ=255)                                                                                                           tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  2:1716770/59‑1 (MQ=255)                                                                                                           tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     >  1:1716770/1‑59 (MQ=255)                                                                                                           tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAACCCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAgg                     >  2:180898/1‑139 (MQ=255)                                                                                                           tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACTGTAACCCTGCGCCAgg                     <  2:1759662/139‑1 (MQ=255)                                                                                                             tgGCACACATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       >  2:1301719/1‑55 (MQ=255)                                                                                                             tgGCACACATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  1:1301719/55‑1 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     <  1:807658/52‑1 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                     >  2:807658/1‑52 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         >  2:715143/1‑112 (MQ=255)                                                                                                                  caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                         <  1:715143/112‑1 (MQ=255)                                                                                                                    tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 <  2:1292963/86‑1 (MQ=255)                                                                                                                    tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 >  1:1292963/1‑86 (MQ=255)                                                                                                                     atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       >  2:723938/1‑47 (MQ=255)                                                                                                                     atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  1:723938/47‑1 (MQ=255)                                                                                                                         cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGt                    <  1:1055774/126‑1 (MQ=255)                                                                                                                         cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGt                    >  2:1055774/1‑126 (MQ=255)                                                                                                                             gTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       <  2:267839/39‑1 (MQ=255)                                                                                                                             gTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                       >  1:267839/1‑39 (MQ=255)                                                                                                                              ttGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATCTAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGaa  >  1:1582558/1‑139 (MQ=255)                                                                                                                                gAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 >  2:537296/1‑74 (MQ=255)                                                                                                                                gAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                 <  1:537296/74‑1 (MQ=255)                                                                                                                                         |                                                                                                                                CTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAA  >  NC_000913/3870305‑3870568

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: ‑