Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,485,806	Δ1 bp	coding (146/1512 nt)	pepA ←	multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,485,803	0	C	.	100.0%	234.7 / NA	50	coding (149/1512 nt)	pepA	multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (36/14);  total (36/14)

GCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTC  >  NC_000913/4485654‑4485948                                                                                                                                                      |                                                                                                                                                  gCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcc                                                                                                                                                  <  2:707866/151‑1 (MQ=255)     aTTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcccgac                                                                                                                                              >  2:37035/1‑148 (MQ=255)        aaTGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGt                                                                                                                                           >  2:816748/1‑150 (MQ=255)         aTGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGt                                                                                                                                           >  2:193611/1‑149 (MQ=255)          tGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGc                                                                                                                                        >  1:415095/1‑151 (MQ=255)               ttCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCa                                                                                                                                       >  2:438973/1‑147 (MQ=255)                tCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGgc                                                                                                                                   >  2:507714/1‑150 (MQ=255)                  tGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCt                                                                                                                                >  1:392635/1‑151 (MQ=255)                      tAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCt                                                                                                                                >  1:849386/1‑147 (MQ=255)                        aCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGt                                                                                                                           >  2:98261/1‑150 (MQ=255)                        aCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTa                                                                                                                          <  2:749437/151‑1 (MQ=255)                         ccTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTAc                                                                                                                         >  1:689548/1‑151 (MQ=255)                            gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTAccc                                                                                                                       >  2:54692/1‑150 (MQ=255)                            gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTAccc                                                                                                                       >  2:634452/1‑150 (MQ=255)                                gTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCg                                                                                                                   <  1:322593/150‑1 (MQ=255)                                  aCTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTg                                                                                                                <  1:609057/151‑1 (MQ=255)                                    tGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGAt                                                                                                              >  1:14363/1‑151 (MQ=255)                                      aCGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGAttt                                                                                                            <  1:8935/151‑1 (MQ=255)                                         atCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTAt                                                                                                         >  1:316240/2‑151 (MQ=255)                                           cATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTAt                                                                                                         <  1:459784/149‑1 (MQ=255)                                            aTCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCg                                                                                                       >  1:187528/1‑150 (MQ=255)                                             tCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAACAGGGCGCTGATGTACCCATCGCTGATTTTATCGa                                                                                                      >  1:680287/1‑150 (MQ=255)                                               cAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCt                                                                                                   <  1:728683/151‑1 (MQ=255)                                               cAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCt                                                                                                   <  2:532992/151‑1 (MQ=255)                                                           tttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCg                                                                                       >  2:81166/1‑151 (MQ=255)                                                           tttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCg                                                                                       >  2:73655/1‑151 (MQ=255)                                                               ccGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGgaga                                                                                   <  1:157570/151‑1 (MQ=255)                                                                 gCAACAAATAGGGAGAATTCGCTCGGAAAGTACATTCGGAAAATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCg                                                                                       <  1:738307/145‑1 (MQ=255)                                                                     ccAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                             >  2:682070/1‑151 (MQ=255)                                                                       aaTAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGAc                                                                           >  2:515883/1‑151 (MQ=255)                                                                          aaGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTg                                                                        >  2:468128/1‑151 (MQ=255)                                                                          aaGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTg                                                                        >  2:629434/1‑151 (MQ=255)                                                                              agaATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTc                                                                    <  1:187021/151‑1 (MQ=255)                                                                              agaATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTc                                                                    >  2:725659/1‑151 (MQ=255)                                                                                aaTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCg                                                                   >  2:716198/1‑150 (MQ=255)                                                                                      cTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTAAG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGAGACGTGGTTCGAAGACg                                                             >  1:548650/1‑150 (MQ=255)                                                                                         ggAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGccc                                                          >  1:232141/1‑150 (MQ=255)                                                                                               tACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgacg                                                    >  2:251931/1‑150 (MQ=255)                                                                                                   ttCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGc                                                >  1:302707/1‑150 (MQ=255)                                                                                                      ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                            >  2:795818/1‑151 (MQ=255)                                                                                                      ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                            >  2:233215/1‑151 (MQ=255)                                                                                                             ggTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGc                                     >  2:332636/1‑151 (MQ=255)                                                                                                                     aacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTAACCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtctc                             >  2:183673/1‑151 (MQ=255)                                                                                                                     aacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCGGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGATGTTtctc                             >  1:786275/1‑151 (MQ=255)                                                                                                                         aTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCggg                         >  1:289607/1‑151 (MQ=255)                                                                                                                                 cccGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCt                 >  2:809019/1‑151 (MQ=255)                                                                                                                                    ggTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttt               <  1:114458/150‑1 (MQ=255)                                                                                                                                     gTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttttt             >  2:124160/1‑151 (MQ=255)                                                                                                                                             tccAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa       <  2:625034/149‑1 (MQ=255)                                                                                                                                              ccAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa     <  1:514266/150‑1 (MQ=255)                                                                                                                                                aGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc  >  1:174954/1‑151 (MQ=255)                                                                                                                                                aGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTc  <  1:679809/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTC  >  NC_000913/4485654‑4485948

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