Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	188,939	C→T	E190K (GAA→AAA)	map ←	methionine aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	188,939	0	C	T	100.0%	201.4 / NA	62	E190K (GAA→AAA)	map	methionine aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/28);  total (34/28)

TCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAATCGCCGCACCGATTTCGCG  >  NC_000913/188803‑189068                                                                                                                                         |                                                                                                                                  tCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTAc                                                                                                                                 >  1:232660/1‑139 (MQ=255)     atTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGa                                                                                                                             >  2:211304/1‑139 (MQ=255)      ttGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACgg                                                                                                                               >  2:1278361/1‑136 (MQ=255)      ttGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACgg                                                                                                                               <  1:1278361/136‑1 (MQ=255)         tgCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTAc                                                                                                                                 >  2:254817/1‑131 (MQ=255)         tgCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTAc                                                                                                                                 <  1:254817/131‑1 (MQ=255)         tgCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCa                                                                                                                         <  2:155403/139‑1 (MQ=255)          gCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCAt                                                                                                                        <  1:688379/139‑1 (MQ=255)           cAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATa                                                                                                                       >  2:830760/1‑139 (MQ=255)           cAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATa                                                                                                                       >  1:765697/1‑139 (MQ=255)              aCAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTg                                                                                                                    >  1:1305127/1‑139 (MQ=255)               cAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTgc                                                                                                                   >  2:331813/1‑139 (MQ=255)                             tGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGttc                                                                                                     <  1:129529/139‑1 (MQ=255)                                     aCCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGATCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGa                                                                                             >  1:86813/1‑139 (MQ=255)                                       cGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAg                                                                                           <  2:232660/139‑1 (MQ=255)                                         tCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                                         <  2:915662/139‑1 (MQ=255)                                         tCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                                         <  1:1172793/139‑1 (MQ=255)                                         tCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                                         <  1:1151386/139‑1 (MQ=255)                                         tCCAGCCATCCTTCATGGGGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCTAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                                         >  2:251192/1‑139 (MQ=255)                                              ccATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGAc                                                                                    >  1:641095/1‑139 (MQ=255)                                                   tttCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCgg                                                                                                        >  2:588531/1‑114 (MQ=255)                                                   tttCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCgg                                                                                                        <  1:588531/114‑1 (MQ=255)                                                   tttCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                                         <  1:439712/129‑1 (MQ=255)                                                   tttCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                                         >  2:439712/1‑129 (MQ=255)                                                     tCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATAcc                                                                             <  1:331813/139‑1 (MQ=255)                                                        tGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATgg                                                                                              <  1:1365011/119‑1 (MQ=255)                                                        tGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATgg                                                                                              >  2:1365011/1‑119 (MQ=255)                                                          gTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTc                                                                        <  1:923279/139‑1 (MQ=255)                                                           tGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTgcag                                                                                                                 >  1:655072/1‑97 (MQ=255)                                                           tGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTgcag                                                                                                                 <  2:655072/97‑1 (MQ=255)                                                           tGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTcc                                                                       >  2:1390160/1‑139 (MQ=255)                                                            gCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTgcagc                                                                                                                >  2:120459/1‑97 (MQ=255)                                                            gCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTgcagc                                                                                                                <  1:120459/97‑1 (MQ=255)                                                                         ttACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaa                                                         >  1:514047/1‑139 (MQ=255)                                                                          tACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaac                                                        >  1:1355767/1‑139 (MQ=255)                                                                               gcgTTGACCATTGGCTCGATGGTGAACGTCATCCCAGTTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGc                                                                                                          <  2:1071948/84‑1 (MQ=255)                                                                               gcgTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGc                                                                                                          >  1:1071948/1‑84 (MQ=255)                                                                                  ttGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCcgc                                                                                       >  1:939226/1‑100 (MQ=255)                                                                                  ttGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCcgc                                                                                       <  2:939226/100‑1 (MQ=255)                                                                                               cGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTc                                   <  1:211304/139‑1 (MQ=255)                                                                                                   ggTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACa                               <  1:251192/139‑1 (MQ=255)                                                                                                       aaCGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGt                                                                                                       >  1:111806/1‑63 (MQ=255)                                                                                                       aaCGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGt                                                                                                       <  2:111806/63‑1 (MQ=255)                                                                                                       aaCGTCATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTAACGAACGACGGAGAAGCCTTCTGCTTCGACAAAtt                           >  2:833338/1‑139 (MQ=255)                                                                                                            cATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaac                                                        >  2:963103/1‑105 (MQ=255)                                                                                                            cATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaac                                                        <  1:963103/105‑1 (MQ=255)                                                                                                            cATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACgg                                                   <  1:1292331/110‑1 (MQ=255)                                                                                                            cATCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACgg                                                   >  2:1292331/1‑110 (MQ=255)                                                                                                              tCCCAGGTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAAt                    >  2:604840/1‑139 (MQ=255)                                                                                                                   ggTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCg                                                                                                         <  1:834095/49‑1 (MQ=255)                                                                                                                   ggTTTCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCg                                                                                                         >  2:834095/1‑49 (MQ=255)                                                                                                                       tCAGTACGACGTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAATCGCCGCAcc           >  1:523878/1‑139 (MQ=255)                                                                                                                         aGTACGACGTTGGTTTTACGGGAGTCATAGTGTAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAATCGCCGCACCGa         >  1:1181666/1‑139 (MQ=255)                                                                                                                              gacgTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCa                                                                    <  1:1165340/75‑1 (MQ=255)                                                                                                                              gacgTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCa                                                                    >  2:1165340/1‑75 (MQ=255)                                                                                                                               acgTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCcgc                                                                                       >  2:495601/1‑55 (MQ=255)                                                                                                                               acgTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCcgc                                                                                       <  1:495601/55‑1 (MQ=255)                                                                                                                                cgTTGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAATCGCCGCACCGATTTcgcg  >  2:930838/1‑139 (MQ=255)                                                                                                                                  ttGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                                         >  1:962793/1‑50 (MQ=255)                                                                                                                                  ttGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                                         <  2:962793/50‑1 (MQ=255)                                                                                                                                   tGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCtt                                    <  1:240011/102‑1 (MQ=255)                                                                                                                                   tGGTTTTACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCtt                                    >  2:240011/1‑102 (MQ=255)                                                                                                                                         |                                                                                                                                  TCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAATCGCCGCACCGATTTCGCG  >  NC_000913/188803‑189068

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