Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,224,989	C→T	R388C (CGT→TGT)	metH →	cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,224,989	0	C	T	100.0%	197.1 / NA	56	R388C (CGT→TGT)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/29);  total (27/29)

GTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGA  >  NC_000913/4224853‑4225118                                                                                                                                         |                                                                                                                                  gTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTc                                                                                                                                >  2:475505/1‑140 (MQ=255) gTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTc                                                                                                                                <  1:486874/140‑1 (MQ=255) gTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTc                                                                                                                                <  1:464039/140‑1 (MQ=255)    cGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAAc                                                                                                                             <  1:780381/140‑1 (MQ=255)    cGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAAc                                                                                                                             <  2:704308/140‑1 (MQ=255)      gCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAg                                                                                                                           >  2:775610/1‑140 (MQ=255)            aGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGaa                                                                                                                     >  2:976676/1‑140 (MQ=255)                 cTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGaaacggc                                                                                                                >  2:464426/1‑136 (MQ=255)                   gAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGgcg                                                                                                              >  2:1113762/1‑140 (MQ=255)                   gAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGgcg                                                                                                              <  1:985767/140‑1 (MQ=255)                     aCATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGgcgcg                                                                                                            <  1:31105/140‑1 (MQ=255)                     aCATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGgcgcg                                                                                                            <  2:627964/140‑1 (MQ=255)                         tGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGa                                                                                                        <  1:986175/140‑1 (MQ=255)                              aaGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGgc                                                                                                               <  1:390431/128‑1 (MQ=255)                              aaGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGgc                                                                                                               >  2:390431/1‑128 (MQ=255)                                gATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGG‑AAACGGCGCGCAGATTATCGat                                                                                                <  1:464426/140‑1 (MQ=255)                                   aGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGatat                                                                                              <  1:775610/140‑1 (MQ=255)                                      cTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCaa                                                                                           <  1:193599/140‑1 (MQ=255)                                      cTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCaa                                                                                           <  1:308216/140‑1 (MQ=255)                                       tGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAAc                                                                                          <  2:981752/140‑1 (MQ=255)                                         tttGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATc                                                                                             <  2:711296/135‑1 (MQ=255)                                         tttGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATc                                                                                             >  1:711296/1‑135 (MQ=255)                                           tgtgAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATgg                                                                                      >  2:886650/1‑140 (MQ=255)                                           tgtgAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATgg                                                                                      >  2:174269/1‑140 (MQ=255)                                             tgAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGAt                                                                                    <  2:1144591/140‑1 (MQ=255)                                                 cGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAg                                                                                >  2:270036/1‑140 (MQ=255)                                                   tGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAggg                                                                              >  2:914863/1‑140 (MQ=255)                                                     ggTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGAt                                                                            >  1:1063399/1‑140 (MQ=255)                                                              aCCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGc                                                                   >  1:50155/1‑140 (MQ=255)                                                              aCCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGc                                                                   >  1:205199/1‑140 (MQ=255)                                                                     tCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAgcg                                                            >  1:109472/1‑140 (MQ=255)                                                                        ccGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAgcggcg                                                         >  1:299585/1‑140 (MQ=255)                                                                        ccGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAgcggcg                                                         <  1:886650/140‑1 (MQ=255)                                                                        ccGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAgcggcg                                                         >  1:712154/1‑140 (MQ=255)                                                                             tCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGt                                                    >  1:1013915/1‑140 (MQ=255)                                                                              ccGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTg                                                   >  2:1072450/1‑140 (MQ=255)                                                                                gatAAGTTCAAGCGTCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGCCGCGCAGACTATCGATATCAATTTGCATAAAGAGGCGCTCGAAGGCGATGCGGCGATggggtg                                                 <  2:308380/138‑5 (MQ=255)                                                                                gCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCg                                                 <  2:1063399/140‑1 (MQ=255)                                                                                gCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCg                                                 <  1:1113762/140‑1 (MQ=255)                                                                                              ccTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGa                                   >  1:637790/1‑140 (MQ=255)                                                                                                tGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGAtt                                 >  1:253260/1‑140 (MQ=255)                                                                                                      aagaagAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCACCAGGTGTAATACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGt                           <  2:253260/140‑1 (MQ=255)                                                                                                      aagaagAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGt                           <  1:475505/140‑1 (MQ=255)                                                                                                      aagaagAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGt                           <  2:712154/140‑1 (MQ=255)                                                                                                      aagaagAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGt                           <  1:854398/140‑1 (MQ=255)                                                                                                      aagaagAGAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGt                           <  1:174269/140‑1 (MQ=255)                                                                                                          agagAAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAAc                       <  1:953428/140‑1 (MQ=255)                                                                                                            agaAATACAGCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCg                     >  2:234991/1‑140 (MQ=255)                                                                                                                     gCGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcgct            >  2:431726/1‑140 (MQ=255)                                                                                                                      cGAGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcgctc           >  2:299603/1‑140 (MQ=255)                                                                                                                        aGGCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcgctcgc         >  1:564132/1‑140 (MQ=255)                                                                                                                         ggCGCTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTcgcg        >  1:512663/1‑140 (MQ=255)                                                                                                                           cgcTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTg      <  2:205199/140‑1 (MQ=255)                                                                                                                           cgcTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTg      >  2:180026/1‑140 (MQ=255)                                                                                                                           cgcTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTg      <  1:270036/140‑1 (MQ=255)                                                                                                                           cgcTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTg      <  1:1072450/140‑1 (MQ=255)                                                                                                                             cTGGATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGcc    <  2:412976/140‑1 (MQ=255)                                                                                                                               ggATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGa  <  1:299603/140‑1 (MQ=255)                                                                                                                               ggATGTCGCGTGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGa  <  2:99782/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  GTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGA  >  NC_000913/4224853‑4225118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑