Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	446,852	A→G	V285A (GTA→GCA)	cyoE ←	heme O synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	446,852	0	A	G	100.0%	182.0 / NA	51	V285A (GTA→GCA)	cyoE	heme O synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/27);  total (24/27)

TAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTA  >  NC_000913/446719‑446986                                                                                                                                      |                                                                                                                                       tAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATaa                                                                                                                                  >  1:923889/1‑140 (MQ=255)  aaaCATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATaaa                                                                                                                                 <  1:281970/140‑1 (MQ=255)     cATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATc                                                                                                                              >  2:889180/1‑140 (MQ=255)       tGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCaa                                                                                                                            <  2:14488/140‑1 (MQ=255)         cGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACg                                                                                                                          >  2:37611/1‑140 (MQ=255)                  tAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACatcatc                                                                                                                 <  1:499169/140‑1 (MQ=255)                                 gagGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTg                                                                                                  >  2:877611/1‑140 (MQ=255)                                   gggATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGAt                                                                                                <  1:815534/140‑1 (MQ=255)                                     gATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATgg                                                                                              <  1:349483/140‑1 (MQ=255)                                     gATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATgg                                                                                              <  1:474707/140‑1 (MQ=255)                                      atatGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                                                                             <  1:1030055/140‑1 (MQ=255)                                              cTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                                                                             <  1:294614/132‑1 (MQ=255)                                              cTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                                                                             >  2:294614/1‑132 (MQ=255)                                              cTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATag                                                                                     >  2:852211/1‑140 (MQ=255)                                               tAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATaga                                                                                    >  2:46527/1‑140 (MQ=255)                                                       aGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                                                                             >  2:756026/1‑123 (MQ=255)                                                       aGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                                                                             <  1:756026/123‑1 (MQ=255)                                                       aGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGa                                                                            >  1:286975/1‑140 (MQ=255)                                                                gTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgc                                                                   <  2:886574/140‑1 (MQ=255)                                                                   ccGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcg                                                                >  2:44570/1‑140 (MQ=255)                                                                   ccGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcg                                                                >  1:359472/1‑140 (MQ=255)                                                                   ccGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTACGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcg                                                                >  2:942901/1‑140 (MQ=255)                                                                    cGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcgc                                                               <  2:129754/140‑1 (MQ=255)                                                                     gTAGCACCTTTTTAATAGAGAGGTTTTGTTCCCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCACCACCGAGAGTGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGcc                                                              <  2:540470/140‑1 (MQ=255)                                                                      tAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGccc                                                             <  2:928797/140‑1 (MQ=255)                                                                      tAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTACGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGccc                                                             <  1:942901/140‑1 (MQ=255)                                                                                      gagagGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGc                                             <  1:37611/140‑1 (MQ=255)                                                                                        gagGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCaa                                           >  1:519742/1‑140 (MQ=255)                                                                                         agGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAAc                                          >  2:791876/1‑140 (MQ=255)                                                                                              ttGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTtat                                     >  2:1034640/1‑140 (MQ=255)                                                                                                     cacacaGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACcgcg                              <  2:782448/140‑1 (MQ=255)                                                                                                     cacacaGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACcgcg                              <  2:89506/140‑1 (MQ=255)                                                                                                      acacaGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCgcgc                             <  1:176491/140‑1 (MQ=255)                                                                                                      acacaGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCgcgc                             <  2:394673/140‑1 (MQ=255)                                                                                                      acacaGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCgcgc                             <  2:430299/140‑1 (MQ=255)                                                                                                       cacaGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCa                            <  2:286975/140‑1 (MQ=255)                                                                                                       cacaGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCa                            <  1:44570/140‑1 (MQ=255)                                                                                                        acaGCAGCCAGCAGCGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCag                           <  1:877611/140‑1 (MQ=255)                                                                                                           gcagccagcagcGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACcgcg                              >  1:164527/1‑134 (MQ=255)                                                                                                           gcagccagcagcGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACcgcg                              <  2:164527/134‑1 (MQ=255)                                                                                                           gcagccagcagcGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGc                        >  1:409315/1‑140 (MQ=255)                                                                                                            cagccagcagcGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCagagtc                       >  1:293784/1‑138 (MQ=255)                                                                                                            cagccagcagcGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGcc                       >  1:185422/1‑140 (MQ=255)                                                                                                              gccagcagcGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCagag                         >  2:878432/1‑136 (MQ=255)                                                                                                              gccagcagcGTATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCagag                         <  1:878432/136‑1 (MQ=255)                                                                                                                        tATGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAaccac           >  2:448239/1‑140 (MQ=255)                                                                                                                         aTGCGAGTCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAaccacc          <  1:791876/140‑1 (MQ=255)                                                                                                                                tCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCt   <  1:46527/140‑1 (MQ=255)                                                                                                                                tCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCt   >  1:355773/1‑140 (MQ=255)                                                                                                                                tCCGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCt   >  1:555686/1‑140 (MQ=255)                                                                                                                                 ccGGTGCCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTa  >  1:91861/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       TAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTA  >  NC_000913/446719‑446986

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