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

CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATT  >  NC_000913/177533‑177794                                                                                                                                  |                                                                                                                                     ccAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTTAGCACCTGAAAAGGTCAGGCGTTAATCTAcc                                                                                                                             >  1:820160/1‑139 (MQ=255) ccAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTAcc                                                                                                                             >  2:555585/1‑139 (MQ=255)      aGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                        <  1:555585/139‑1 (MQ=255)        ggTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                      <  2:820160/139‑1 (MQ=255)          tCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTAcc                                                                                                                             <  2:1382146/130‑1 (MQ=255)          tCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTAcc                                                                                                                             >  1:1382146/1‑130 (MQ=255)                     ttttCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                         <  2:425681/139‑1 (MQ=255)                                 cAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                             <  1:351673/139‑1 (MQ=255)                                           gagaTGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGataata                                                                                   >  2:1271077/1‑139 (MQ=255)                                                 gCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTggg                                                                             <  1:524273/139‑1 (MQ=255)                                                      tACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                             <  2:160164/118‑1 (MQ=255)                                                      tACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                             >  1:160164/1‑118 (MQ=255)                                                           cTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCa                                                                                                                 >  2:1206670/1‑93 (MQ=255)                                                           cTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCa                                                                                                                 <  1:1206670/93‑1 (MQ=255)                                                                 cGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGCGGGCTTGCACGTTCAAAcc                                                             <  1:1022719/139‑1 (MQ=255)                                                                   aCTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTATCTGTGAAAGCGCATTACAGAGCTGTTGTGc                                                                                                <  1:532291/102‑1 (MQ=255)                                                                   aCTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGc                                                                                                >  2:532291/1‑102 (MQ=255)                                                                        atCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACt                                                      >  1:1195260/1‑139 (MQ=255)                                                                             gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                        <  2:1124540/68‑1 (MQ=255)                                                                             gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                        >  1:1124540/1‑68 (MQ=255)                                                                             gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTga                                                 >  2:193830/1‑139 (MQ=255)                                                                              ccAATAGACGAGCATTGTTCATTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACag                                                                                                           >  2:1131552/1‑80 (MQ=255)                                                                              ccAATAGACGAGCATTGTTCATTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACag                                                                                                           <  1:1131552/80‑1 (MQ=255)                                                                              ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTggg                                                                             <  2:75063/110‑1 (MQ=255)                                                                              ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTggg                                                                             >  1:75063/1‑110 (MQ=255)                                                                                 aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                      <  2:1089554/66‑1 (MQ=255)                                                                                 aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGaaa                                                                                                                      >  1:1089554/1‑66 (MQ=255)                                                                                 aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAgg                                             >  1:431283/1‑139 (MQ=255)                                                                                        aGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCt                                                                                                                                >  2:270745/1‑49 (MQ=255)                                                                                        aGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCt                                                                                                                                <  1:270745/49‑1 (MQ=255)                                                                                         gCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACgg                                     <  2:1195260/139‑1 (MQ=255)                                                                                          cATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                             <  1:50428/82‑1 (MQ=255)                                                                                          cATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                             >  2:50428/1‑82 (MQ=255)                                                                                            ttGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGaa                                  >  2:457513/1‑139 (MQ=255)                                                                                               ttCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  <  2:595816/40‑1 (MQ=255)                                                                                               ttCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  <  2:60102/40‑1 (MQ=255)                                                                                               ttCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  >  1:60102/1‑40 (MQ=255)                                                                                               ttCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                  >  1:595816/1‑40 (MQ=255)                                                                                                   aTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           <  1:529546/139‑1 (MQ=255)                                                                                                   aTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCAGTACAGAGCTGTTGTGCAGCGAGGATATTACGTGGGCTTGCACGTTCAAACCAGTCACGCGTGAGAGGAATAACGGGAATTTTGAg                           <  2:431283/139‑1 (MQ=255)                                                                                                        cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                      <  1:193830/139‑1 (MQ=255)                                                                                                         aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                         <  1:1241519/55‑1 (MQ=255)                                                                                                         aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                         >  2:1241519/1‑55 (MQ=255)                                                                                                         aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                            <  2:1309694/68‑1 (MQ=255)                                                                                                         aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                            >  1:1309694/1‑68 (MQ=255)                                                                                                              gCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAcccc                >  1:1294167/1‑139 (MQ=255)                                                                                                                  cTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTa            >  1:1020178/1‑139 (MQ=255)                                                                                                                          gTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           <  2:115963/116‑1 (MQ=255)                                                                                                                          gTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                           >  1:115963/1‑116 (MQ=255)                                                                                                                            cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGc                                                                                                <  2:1088499/45‑1 (MQ=255)                                                                                                                            cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGc                                                                                                >  1:1088499/1‑45 (MQ=255)                                                                                                                            cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                            <  2:232818/49‑1 (MQ=255)                                                                                                                            cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                            >  1:232818/1‑49 (MQ=255)                                                                                                                            cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTggg                                                                             <  2:1096706/64‑1 (MQ=255)                                                                                                                            cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTggg                                                                             >  1:1096706/1‑64 (MQ=255)                                                                                                                            cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGAtt  >  2:1311387/1‑139 (MQ=255)                                                                                                                               gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                            >  2:1298541/1‑46 (MQ=255)                                                                                                                                  |                                                                                                                                     CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATT  >  NC_000913/177533‑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: ‑