Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	277,026	A→G	E104E (GAA→GAG)	mmuM →	homocysteine S‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	277,026	0	A	G	100.0%	99.7 / NA	30	E104E (GAA→GAG)	mmuM	homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/20);  total (10/20)

ACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTATCATTGTAGCGTTGAGG  >  NC_000913/276890‑277147                                                                                                                                         |                                                                                                                          acCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGc                                                                                                                         >  2:85930/1‑139 (MQ=255) acCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCg                                                                                                                        <  2:435496/140‑1 (MQ=255)    gggcgggGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTAt                                                                                                                     <  2:190802/140‑1 (MQ=255)    gggcgggGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTAt                                                                                                                     <  2:184176/140‑1 (MQ=255)         ggggCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGc                                                                                                                >  2:71435/1‑140 (MQ=255)            gcgcAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCaga                                                                                                             <  2:142278/140‑1 (MQ=255)                    cgcgATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGGCTGGATGAAGCGCAGGCGAAAGAGCTGATTTGCAAAAGCGCGGAACTGGCGCGTAAAACCCGTGAGGCGTCTCTGGCAGAGAACCCGc                                                                                                     >  2:362511/1‑140 (MQ=255)                        aTCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGc                                                                                                 <  1:635199/140‑1 (MQ=255)                             tGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGaa                                                                                            >  2:112938/1‑140 (MQ=255)                                   cTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTc                                                                                      >  1:36843/1‑140 (MQ=255)                                   cTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTc                                                                                      >  1:65079/1‑140 (MQ=255)                                     aTCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCtg                                                                                    >  1:276975/1‑140 (MQ=255)                                              cGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGGACCCTTCTGGTGGCgggg                                                                           >  2:368112/1‑139 (MQ=255)                                              cGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGa                                                                           <  1:71435/140‑1 (MQ=255)                                              cGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGa                                                                           <  1:726613/140‑1 (MQ=255)                                                           cgccgcGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCtt                                                              >  2:338292/1‑140 (MQ=255)                                                                         tGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTg                                                <  2:35820/140‑1 (MQ=255)                                                                              gAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGa                                           <  1:338292/140‑1 (MQ=255)                                                                                aGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATg                                         <  1:368112/140‑1 (MQ=255)                                                                                 gcgcAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGGTCCGTCGGGGCTTTCGGCGCGGATCTGgggggggg                                        >  2:380553/1‑136 (MQ=255)                                                                                 gcgcAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATgg                                        >  1:356377/1‑140 (MQ=255)                                                                                  cgcAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGc                                       <  1:289260/140‑1 (MQ=255)                                                                                            aGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGa                                           <  1:275322/126‑1 (MQ=255)                                                                                            aGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGa                                           >  2:275322/1‑126 (MQ=255)                                                                                               gcTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTg                          <  2:36843/140‑1 (MQ=255)                                                                                                cTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTgg                         <  2:276975/140‑1 (MQ=255)                                                                                                  gATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCg                       <  1:380553/140‑1 (MQ=255)                                                                                                  gATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCg                       <  1:362511/140‑1 (MQ=255)                                                                                                      ggCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTa                   <  1:112938/140‑1 (MQ=255)                                                                                                       gCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTAt                  <  2:65079/140‑1 (MQ=255)                                                                                                         aaaaGCGTGGAGCTGGCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTATCa                <  2:55096/140‑1 (MQ=255)                                                                                                                       ggCGCGTAAAGCCCGTGAGGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTATCATTGTAGCGTTGAgg  <  1:782118/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                          ACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCGCGTATCTGGCGGATGGCTCTGAATACCGTGGCGATTATCATTGTAGCGTTGAGG  >  NC_000913/276890‑277147

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