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	59	*267* (TAG→TAA)	btuF	vitamin B12 ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/30);  total (29/30)

CACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATT  >  NC_000913/177531‑177794                                                                                                                                    |                                                                                                                                     cacCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  <  2:1565997/136‑1 (MQ=255) cacCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  >  1:1565997/1‑136 (MQ=255)   ccAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  <  2:114201/134‑1 (MQ=255)   ccAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  >  1:114201/1‑134 (MQ=255)   ccAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTAcc                                                                                                                             >  2:642257/1‑139 (MQ=255)        aGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                        <  2:1256553/139‑1 (MQ=255)        aGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                        <  2:1433074/139‑1 (MQ=255)                tACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  >  2:315893/1‑121 (MQ=255)                tACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  <  1:315893/121‑1 (MQ=255)                    cAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACa                                                                                                            <  1:1404193/139‑1 (MQ=255)                       ttttCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                         >  2:1271292/1‑139 (MQ=255)                                             gagaTGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  <  1:1002726/92‑1 (MQ=255)                                             gagaTGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  >  2:1002726/1‑92 (MQ=255)                                                        tACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGc                                                                        >  1:103274/1‑139 (MQ=255)                                                             cTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                      <  2:480656/88‑1 (MQ=255)                                                             cTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                      >  1:480656/1‑88 (MQ=255)                                                                 gCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCaaa                                                               >  1:172586/1‑139 (MQ=255)                                                                 gCCGACTATATCCAGCCAATAGACGACCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGTTTGCACGTTCaaa                                                               <  1:1064674/139‑1 (MQ=255)                                                                     aCTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGCAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGCGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAg                                                           >  1:960953/1‑139 (MQ=255)                                                                     aCTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAg                                                           >  1:1155872/1‑139 (MQ=255)                                                                          atCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACt                                                      >  2:57878/1‑139 (MQ=255)                                                                               gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTAc                                                                                                                              <  1:1447183/62‑1 (MQ=255)                                                                               gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTAc                                                                                                                              >  2:1447183/1‑62 (MQ=255)                                                                               gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGTTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTga                                                 >  2:667470/1‑139 (MQ=255)                                                                               gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTga                                                 >  1:680366/1‑139 (MQ=255)                                                                                ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                         <  2:868599/82‑1 (MQ=255)                                                                                ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                         >  1:868599/1‑82 (MQ=255)                                                                                ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATGCGTGGGCTTGCACGTTCAAACCAGTCACTCGTgag                                                >  2:1009497/1‑139 (MQ=255)                                                                                ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTgag                                                >  2:1636849/1‑139 (MQ=255)                                                                                   aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                      <  1:748390/66‑1 (MQ=255)                                                                                   aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                      >  2:748390/1‑66 (MQ=255)                                                                                          aGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                            >  1:1264090/1‑85 (MQ=255)                                                                                          aGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                            <  2:1264090/85‑1 (MQ=255)                                                                                                     aTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           <  2:680366/139‑1 (MQ=255)                                                                                                     aTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           <  2:103274/139‑1 (MQ=255)                                                                                                     aTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           <  1:1271292/139‑1 (MQ=255)                                                                                                      tCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcg                                                                                                                   <  2:1410/50‑1 (MQ=255)                                                                                                      tCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcg                                                                                                                   >  1:1410/1‑50 (MQ=255)                                                                                                          cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                      <  2:960953/139‑1 (MQ=255)                                                                                                           aCTGAGCACCTGAAAAGGTCATGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                      <  1:915720/138‑1 (MQ=255)                                                                                                           aCTGAGCACCTGAAAAGGTCATGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                      >  2:915720/1‑138 (MQ=255)                                                                                                            cTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGTTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCa                    <  1:667470/139‑1 (MQ=255)                                                                                                             tGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAc                   <  1:1009497/139‑1 (MQ=255)                                                                                                              gAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAc                   <  2:114693/138‑1 (MQ=255)                                                                                                              gAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAc                   >  1:114693/1‑138 (MQ=255)                                                                                                               aGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAccc                 <  1:57878/139‑1 (MQ=255)                                                                                                                gCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGc                                                                        >  1:676294/1‑83 (MQ=255)                                                                                                                gCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGc                                                                        <  2:676294/83‑1 (MQ=255)                                                                                                                gCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           >  1:612632/1‑128 (MQ=255)                                                                                                                gCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           <  2:612632/128‑1 (MQ=255)                                                                                                                 cACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                             <  2:1179435/61‑1 (MQ=255)                                                                                                                 cACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTAAAGAGCTGTTGTGCAGc                                                                                             >  1:1179435/1‑61 (MQ=255)                                                                                                                   ccTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGt             <  1:1636849/139‑1 (MQ=255)                                                                                                                    cTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTg                                                                               <  2:391926/72‑1 (MQ=255)                                                                                                                    cTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTg                                                                               >  1:391926/1‑72 (MQ=255)                                                                                                                         aaGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGtt       <  2:1155872/139‑1 (MQ=255)                                                                                                                              cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           <  2:44732/114‑1 (MQ=255)                                                                                                                              cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           >  1:44732/1‑114 (MQ=255)                                                                                                                              cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGAtt  >  1:863062/1‑139 (MQ=255)                                                                                                                                 gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                          >  1:1551771/1‑80 (MQ=255)                                                                                                                                 gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                          <  2:1551771/80‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     CACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATT  >  NC_000913/177531‑177794

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