Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	177,662	C→T	*267* (TAG→TAA)	btuF ←	vitamin B12 ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	177,662	0	C	T	100.0%	188.8 / NA	58	*267* (TAG→TAA)	btuF	vitamin B12 ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/31);  total (27/31)

AACACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTA  >  NC_000913/177529‑177797                                                                                                                                      |                                                                                                                                        aaCACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATc                                                                                                                                    >  2:290271/1‑139 (MQ=255) aaCACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATc                                                                                                                                    >  1:906/1‑139 (MQ=255)          aGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                           <  2:906/139‑1 (MQ=255)          aGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                           <  1:290271/139‑1 (MQ=255)              tCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgc                                                                                                                       >  1:52537/1‑139 (MQ=255)              tCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGACAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgc                                                                                                                       >  1:342479/1‑139 (MQ=255)               ccATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcg                                                                                                                      <  1:1169513/139‑1 (MQ=255)                cATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                           >  1:245606/1‑133 (MQ=255)                cATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                           <  2:245606/133‑1 (MQ=255)                    cGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTa                                                                                                                 <  1:459697/139‑1 (MQ=255)                       aaTTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACag                                                                                                              <  1:920021/139‑1 (MQ=255)                       aaTTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACag                                                                                                              >  1:976890/1‑139 (MQ=255)                       aaTTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACag                                                                                                              <  1:270630/139‑1 (MQ=255)                       aaTTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACag                                                                                                              >  2:267254/1‑139 (MQ=255)                         ttttCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                            <  2:52537/139‑1 (MQ=255)                         ttttCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                            <  1:1154131/139‑1 (MQ=255)                               ggCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                            <  2:649586/133‑1 (MQ=255)                               ggCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                            >  1:649586/1‑133 (MQ=255)                                     cAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                                <  2:561623/139‑1 (MQ=255)                                        taCGCCGTATATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAg                                                                                             <  1:451245/138‑1 (MQ=255)                                                 gaTGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACg                                                                                    <  1:267254/139‑1 (MQ=255)                                                     gCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                         <  2:1021903/98‑1 (MQ=255)                                                     gCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                         >  1:1021903/1‑98 (MQ=255)                                                     gCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTggg                                                                                <  2:1123511/139‑1 (MQ=255)                                                               cTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGtt                                                                      >  1:48466/1‑139 (MQ=255)                                                               cTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGtt                                                                      >  1:834517/1‑139 (MQ=255)                                                                 gtgCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCa                                                                    <  2:342479/139‑1 (MQ=255)                                                                            atCCAGCCAATAGACGAGCATTGTTCAATCCCCACTTAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcg                                                                                                                      <  1:843101/78‑1 (MQ=255)                                                                            atCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgca                                                                                                                      >  2:843101/1‑77 (MQ=255)                                                                               cAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGt                                                      <  1:92519/139‑1 (MQ=255)                                                                               cAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGt                                                      <  2:834517/139‑1 (MQ=255)                                                                                 gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                     <  2:1135160/58‑1 (MQ=255)                                                                                 gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                     >  1:1135160/1‑58 (MQ=255)                                                                                  ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcg                                                                                                                      >  2:799822/1‑72 (MQ=255)                                                                                  ccAACAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcg                                                                                                                      <  1:799822/72‑1 (MQ=255)                                                                                         aCGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTgag                                                   >  2:475799/1‑132 (MQ=255)                                                                                         aCGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTgag                                                   <  1:475799/132‑1 (MQ=255)                                                                                             gCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                               <  2:239301/84‑1 (MQ=255)                                                                                             gCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                               >  1:239301/1‑84 (MQ=255)                                                                                              cATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                               >  1:413916/1‑83 (MQ=255)                                                                                              cATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                               <  2:413916/83‑1 (MQ=255)                                                                                                   ttCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCa                                                                          <  1:1207109/99‑1 (MQ=255)                                                                                                   ttCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCa                                                                          >  2:1207109/1‑99 (MQ=255)                                                                                                     cAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                     <  2:763542/38‑1 (MQ=255)                                                                                                     cAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                     >  1:763542/1‑38 (MQ=255)                                                                                                       aTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                              <  2:48466/139‑1 (MQ=255)                                                                                                        tCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             >  1:593390/1‑107 (MQ=255)                                                                                                        tCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             <  2:593390/107‑1 (MQ=255)                                                                                                        tCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGc                             >  2:579411/1‑139 (MQ=255)                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             <  1:1257280/102‑1 (MQ=255)                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             >  2:1257280/1‑102 (MQ=255)                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAAtttt                                 <  1:771430/130‑1 (MQ=255)                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAAtttt                                 >  2:771430/1‑130 (MQ=255)                                                                                                                 aGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAccc                    >  1:1082128/1‑139 (MQ=255)                                                                                                                   cACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCa                  >  1:754964/1‑139 (MQ=255)                                                                                                                      cTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTa               >  2:350579/1‑139 (MQ=255)                                                                                                                                cAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                            <  2:789065/36‑1 (MQ=38)                                                                                                                                cAGGCGTTAATCTACCTGTGAAAGCGCATAACagag                                                                                                            >  1:789065/1‑36 (MQ=255)                                                                                                                                   gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTa  >  1:483599/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        AACACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTA  >  NC_000913/177529‑177797

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