Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	954,813	(A)7→6	intergenic (‑347/+59)	focA ← / ← ycaO	formate channel FocA/ribosomal protein S12 methylthiotransferase accessory factor YcaO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,807	0	A	.	100.0%	188.1 / NA	44	intergenic (‑341/+65)	focA/ycaO	formate channel FocA/ribosomal protein S12 methylthiotransferase accessory factor YcaO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (21/23);  total (21/23)

AACAGTTAATTAAAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGATGTGCA  >  NC_000913/954673‑954941                                                                                                                                       |                                                                                                                                       aaCAGTTAATTAAAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaa                                                                                                                                   >  1:317208/1‑140 (MQ=255)  aCAGTTAATTAAAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaac                                                                                                                                  >  1:298883/1‑139 (MQ=255)   cAGTTAATTAAAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaacg                                                                                                                                 <  1:590190/140‑3 (MQ=255)   cAGTTAATTAAAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaacg                                                                                                                                 <  2:8222/140‑3 (MQ=255)   cAGTTAATTAAAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaacg                                                                                                                                 >  1:342017/1‑138 (MQ=255)      ttaattaaAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGccc                                                                                                                             <  2:728901/140‑1 (MQ=255)        aattaaAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCg                                                                                                                           <  2:379965/140‑1 (MQ=255)        aattaaAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCg                                                                                                                           <  2:342017/140‑1 (MQ=255)        aattaaAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCg                                                                                                                           <  2:298883/140‑1 (MQ=255)        aattaaAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCg                                                                                                                           <  1:1122838/140‑1 (MQ=255)              aaGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAatat                                                                                                                     >  2:172770/1‑140 (MQ=255)                  gAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACt                                                                                                                 <  2:317208/140‑1 (MQ=255)                  gAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACt                                                                                                                 <  1:967136/140‑1 (MQ=255)                    gCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCa                                                                                                               >  1:839914/1‑140 (MQ=255)                    gCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCa                                                                                                               >  1:768572/1‑140 (MQ=255)                    gCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCa                                                                                                               >  2:485550/1‑140 (MQ=255)                       tCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGAc                                                                                                            <  1:1098917/140‑1 (MQ=255)                              aaTAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTCTTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCAt                                                                                                     <  1:1022840/140‑1 (MQ=255)                              aaTAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCAt                                                                                                     <  1:693096/140‑1 (MQ=255)                              aaTAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCAt                                                                                                     <  1:172770/140‑1 (MQ=255)                                tAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTa                                                                                                          <  1:741909/133‑1 (MQ=255)                                tAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTa                                                                                                          >  2:741909/1‑133 (MQ=255)                                   aGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAACATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGc                                                                                                <  2:808900/140‑1 (MQ=255)                                       aCAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTAc                                                                                            <  1:328361/140‑1 (MQ=255)                                        cAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACa                                                                                           <  2:503853/140‑1 (MQ=255)                                        cAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACa                                                                                           <  2:625951/140‑1 (MQ=255)                                        cAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACa                                                                                           <  1:882168/140‑1 (MQ=255)                                             tcgatcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCg                                                                                      >  1:468707/1‑140 (MQ=255)                                             tcgatcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCg                                                                                      >  1:195305/1‑140 (MQ=255)                                                               tACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                                                    >  2:671807/1‑140 (MQ=255)                                                               tACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                                                    >  2:361022/1‑140 (MQ=255)                                                                aCAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatt                                                                   <  2:1126295/140‑1 (MQ=255)                                                                                              cAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCa                                     >  2:453406/1‑140 (MQ=255)                                                                                                 ttattttattAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCt                                  <  2:513342/140‑1 (MQ=255)                                                                                                        attaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTa                           >  2:354832/1‑140 (MQ=255)                                                                                                         ttaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTg                                       <  1:1643/127‑1 (MQ=255)                                                                                                         ttaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTg                                       >  2:1643/1‑127 (MQ=255)                                                                                                            aTGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGcc                       >  2:827683/1‑140 (MQ=255)                                                                                                                 aaTATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTcaa                  >  1:499145/1‑140 (MQ=255)                                                                                                                  atatatGTAAAGCGGGGCATT‑ATAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTcaac                 <  1:827683/140‑1 (MQ=255)                                                                                                                  atatatGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTcaac                 <  1:97224/140‑1 (MQ=255)                                                                                                                       tGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGa            >  2:1121956/1‑140 (MQ=255)                                                                                                                       tGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGa            >  1:273653/1‑140 (MQ=255)                                                                                                                          aaaGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTg         <  2:839914/140‑1 (MQ=255)                                                                                                                            aGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGAt       >  1:68781/1‑140 (MQ=255)                                                                                                                            aGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGAt       >  1:807628/1‑140 (MQ=255)                                                                                                                                gggCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGATGTgc   >  1:276273/1‑140 (MQ=255)                                                                                                                                gggCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGATGTgc   >  1:594258/1‑140 (MQ=255)                                                                                                                                 ggCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGATGTgca  >  1:755998/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       AACAGTTAATTAAAAGGGAGCATCAGGCGAATAAAGTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGATGTGCA  >  NC_000913/954673‑954941

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