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%	204.0 / NA	43	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 . (30/13);  total (30/13)

GGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAA  >  NC_000913/4485664‑4485944                                                                                                                                            |                                                                                                                                              ggTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCa                                                                                                                                   <  1:131303/151‑1 (MQ=255)    tttCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAggg                                                                                                                                >  1:293227/1‑151 (MQ=255)           aTAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGAt                                                                                                                         >  1:174072/1‑151 (MQ=255)             aaCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGt                                                                                                                       <  2:293513/151‑1 (MQ=255)             aaCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGt                                                                                                                       >  2:442296/1‑151 (MQ=255)             aaCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGt                                                                                                                       >  2:561609/1‑151 (MQ=255)                cTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTAcc                                                                                                                    >  2:711796/1‑151 (MQ=255)                        aCTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCt                                                                                                             <  2:244033/150‑1 (MQ=255)                            aCGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGAttt                                                                                                        >  2:16458/1‑151 (MQ=255)                                   tCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCg                                                                                                   >  2:723041/1‑149 (MQ=255)                                         tCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTc                                                                                           >  2:767794/1‑151 (MQ=255)                                          cACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGtt                                                                                            <  1:754981/149‑1 (MQ=255)                                          cACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCt                                                                                          <  2:590831/151‑1 (MQ=255)                                             gTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGc                                                                                        >  2:216605/1‑150 (MQ=255)                                             gTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGc                                                                                        >  1:264819/1‑150 (MQ=255)                                                cTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATc                                                                                    >  1:419242/1‑151 (MQ=255)                                                     ccGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                  >  1:655619/1‑148 (MQ=255)                                                           ccAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGc                                                                          <  1:183518/150‑1 (MQ=255)                                                                aaGGGGAATTCGCTCGGACAGTAAATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTGCG‑CCCGACGTAGCAGTGCGCTGATGTACCCCTCGCAGATATAATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGt                                                                     >  1:200453/1‑150 (MQ=255)                                                                      aaTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCg                                                               >  2:626173/1‑150 (MQ=255)                                                                      aaTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCg                                                               >  1:222822/1‑150 (MQ=255)                                                                         tCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAg                                                            <  1:546194/150‑1 (MQ=255)                                                                          cGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGAc                                                          >  2:757922/1‑151 (MQ=255)                                                                                      aCATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGAt                                              >  2:449832/1‑151 (MQ=255)                                                                                       cATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATg                                             <  2:305289/151‑1 (MQ=255)                                                                                          tCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCa                                           <  1:458069/150‑1 (MQ=255)                                                                                           cGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAgg                                         >  1:144395/1‑151 (MQ=255)                                                                                            ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                        >  2:638844/1‑151 (MQ=255)                                                                                            ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                        >  1:151522/1‑151 (MQ=255)                                                                                              aaCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCAc                                      <  1:745424/151‑1 (MQ=255)                                                                                               aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCGCACGACGATGCAGGCACt                                     <  1:710744/151‑1 (MQ=255)                                                                                                 aTGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTcc                                   >  2:664612/1‑151 (MQ=255)                                                                                                 aTGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTcc                                   >  1:231058/1‑151 (MQ=255)                                                                                                   ggTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGc                                 >  1:428471/1‑151 (MQ=255)                                                                                                      gcagcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTg                               <  2:419293/150‑1 (MQ=255)                                                                                                              aaTGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCgg                      >  1:558890/1‑151 (MQ=255)                                                                                                                tGTCTGCCCCGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGc                    >  1:44885/1‑151 (MQ=255)                                                                                                                        ccGGTTTTCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCtt            >  1:165910/1‑151 (MQ=255)                                                                                                                            ttttCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttttt         >  2:372325/1‑150 (MQ=255)                                                                                                                            ttttCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttttt         <  1:626116/150‑1 (MQ=255)                                                                                                                            ttttCCTTCCAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCtttt          >  1:471117/1‑149 (MQ=255)                                                                                                                                ccttccAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt    >  1:717535/1‑151 (MQ=255)                                                                                                                                   tccAGTTCG‑CCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa  >  1:509689/1‑150 (MQ=255)                                                                                                                                            |                                                                                                                                              GGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAA  >  NC_000913/4485664‑4485944

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