Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,953,927	A→G	E162E (GAA→GAG)	cmoB →	tRNA U34 carboxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,953,927	0	A	G	100.0%	196.6 / NA	54	E162E (GAA→GAG)	cmoB	tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/30);  total (24/30)

CATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGC  >  NC_000913/1953790‑1954063                                                                                                                                          |                                                                                                                                         cATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGc                                                                                                                                        >  1:595548/1‑140 (MQ=255) cATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGc                                                                                                                                        >  1:1048498/1‑140 (MQ=255) cATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGc                                                                                                                                        >  2:651008/1‑140 (MQ=255) cATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGc                                                                                                                                        >  1:466591/1‑140 (MQ=255)  aTCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCg                                                                                                                                       <  2:923382/140‑1 (MQ=255)  aTCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCg                                                                                                                                       <  2:898110/140‑1 (MQ=255)  aTCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCg                                                                                                                                       <  2:595548/140‑1 (MQ=255)      ttCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGc                                                                                                                                   <  2:466591/140‑1 (MQ=255)        cTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGt                                                                                                                                 >  1:211766/1‑140 (MQ=255)         tGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTa                                                                                                                                <  2:130969/140‑1 (MQ=255)                aCCGGACGCACCATTCTGGATGTCGGCTGCGGCATCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGGGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAactact                                                                                                                         <  1:1034786/140‑1 (MQ=255)                aCCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGt                                                                                                                                     <  1:1165928/128‑1 (MQ=255)                aCCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGt                                                                                                                                     >  2:1165928/1‑128 (MQ=255)                aCCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAactact                                                                                                                         <  2:448776/140‑1 (MQ=255)                aCCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAactact                                                                                                                         <  2:185425/140‑1 (MQ=255)                aCCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAactact                                                                                                                         <  2:547314/140‑1 (MQ=255)                aCCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAactact                                                                                                                         <  2:984911/140‑1 (MQ=255)                aCCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAactact                                                                                                                         >  1:564060/1‑140 (MQ=255)                          ccATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGAt                                                                                                               <  2:335562/140‑1 (MQ=255)                                  gATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAgcgcgc                                                                                                       <  2:211766/140‑1 (MQ=255)                                   aTGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAgcgcgcg                                                                                                      >  1:797662/1‑140 (MQ=255)                                     gTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTaa                                                                                                                   <  2:525659/125‑1 (MQ=255)                                     gTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTaa                                                                                                                   >  1:525659/1‑125 (MQ=255)                                      tCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCAt                                                                                                   <  2:564060/140‑1 (MQ=255)                                       cGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCAtt                                                                                                  >  2:1013527/1‑140 (MQ=255)                                          cTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGt                                                                                               <  1:60062/140‑1 (MQ=255)                                          cTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGt                                                                                               <  2:1048498/140‑1 (MQ=255)                                            gCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTa                                                                                             >  1:795140/1‑140 (MQ=255)                                                    ggTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAgg                                                                                     <  2:1171032/140‑1 (MQ=255)                                                             aTGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGaaca                                                                            <  2:41341/140‑1 (MQ=255)                                                                tgGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACt                                                                         >  2:1061679/1‑140 (MQ=255)                                                                tgGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACt                                                                         >  1:276966/1‑140 (MQ=255)                                                                  gcgcATGATTGGCGCAGGCGCGCTCCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTc                                                                       <  2:817241/140‑1 (MQ=255)                                                                  gcgcATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTc                                                                       >  1:817241/1‑140 (MQ=255)                                                                  gcgcATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTc                                                                       >  1:902137/1‑140 (MQ=255)                                                                    gcATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCg                                                                     >  1:947590/1‑140 (MQ=255)                                                                     cATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCgg                                                                    <  2:947590/140‑1 (MQ=255)                                                                             gcgcAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGaaa                                                            <  2:902137/140‑1 (MQ=255)                                                                                        cACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATAc                                                 <  1:469148/140‑1 (MQ=255)                                                                                             cgcgGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCt                                            >  1:1212247/1‑140 (MQ=255)                                                                                                   ggggATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGa                                      <  1:1145599/140‑1 (MQ=255)                                                                                                       aTCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATggg                                  <  1:1061679/140‑1 (MQ=255)                                                                                                        tCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGc                                 <  1:1013527/140‑1 (MQ=255)                                                                                                              cTACGCAACTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCtt                           >  2:839428/1‑140 (MQ=255)                                                                                                                      cTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCg                   >  2:110880/1‑140 (MQ=255)                                                                                                                      cTGTTTCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCg                   >  1:292672/1‑140 (MQ=255)                                                                                                                         tttCTGTGCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCg                >  2:1178522/1‑140 (MQ=255)                                                                                                                             tgtgCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCa            <  1:110880/140‑1 (MQ=255)                                                                                                                                gCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCg         <  1:1178522/140‑1 (MQ=255)                                                                                                                                gCCAGTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCg         >  1:31248/1‑140 (MQ=255)                                                                                                                                    gTTTGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTgg     >  1:990875/1‑140 (MQ=255)                                                                                                                                     tttGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGa    <  2:795140/140‑1 (MQ=255)                                                                                                                                     tttGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGa    <  2:797662/140‑1 (MQ=255)                                                                                                                                       tGAGGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGc  >  2:1115839/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                         CATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGC  >  NC_000913/1953790‑1954063

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