Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,053,565	(T)5→4	coding (1276/1326 nt)	pepP ←	proline aminopeptidase P II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,053,561	0	T	.	100.0%	198.9 / NA	47	coding (1280/1326 nt)	pepP	proline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (22/25);  total (22/25)

TCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTAC  >  NC_000913/3053425‑3053694                                                                                                                                         |                                                                                                                                      tcaaATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCtttt                                                                                                                                    >  2:338780/1‑140 (MQ=255)    aaTGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCttttcac                                                                                                                                 >  2:851127/1‑137 (MQ=255)         aCCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGc                                                                                                                           <  1:814761/140‑1 (MQ=255)          ccGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCt                                                                                                                          >  2:998064/1‑140 (MQ=255)          ccGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCt                                                                                                                          >  1:173240/1‑140 (MQ=255)          ccGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCt                                                                                                                          >  2:511456/1‑140 (MQ=255)           cGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTg                                                                                                                         >  1:550506/1‑140 (MQ=255)               aGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCgg                                                                                                                     <  1:998064/140‑1 (MQ=255)               aGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCgg                                                                                                                     <  2:133206/140‑1 (MQ=255)                 cgcCCCGTGACTTAACGGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTg                                                                                                                   <  1:438003/140‑1 (MQ=255)                   ccccGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAg                                                                                                                 >  2:1002132/1‑140 (MQ=255)                       gTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGttt                                                                                                             <  2:707070/140‑1 (MQ=255)                        tGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGtttt                                                                                                            >  2:793409/1‑140 (MQ=255)                                   ggAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTa                                                                                                       >  2:400386/1‑134 (MQ=255)                                   ggAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTa                                                                                                       <  1:400386/134‑1 (MQ=255)                                   ggAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGtt                                                                                                 >  2:68554/1‑140 (MQ=255)                                          gccagcgccagcGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTaa                                                                                          >  1:606103/1‑140 (MQ=255)                                           ccaccgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAAt                                                                                         <  2:316367/140‑1 (MQ=255)                                              gcgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATcac                                                                                      >  2:410095/1‑140 (MQ=255)                                               cgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATcaca                                                                                     <  2:606103/140‑1 (MQ=255)                                               cgccagcgTCGCCCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATcaca                                                                                     <  2:157289/140‑1 (MQ=255)                                                    gcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATgtc                                                                                >  2:975027/1‑140 (MQ=255)                                                      gTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATgtcgt                                                                              <  1:457735/140‑1 (MQ=255)                                                      gTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATgtcgt                                                                              <  1:97709/140‑1 (MQ=255)                                                          cgcccgccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTc                                                                          <  1:851127/140‑1 (MQ=255)                                                                ccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACg                                                                    <  2:621799/140‑1 (MQ=255)                                                                 cATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGa                                                                   >  2:352974/1‑140 (MQ=255)                                                                             gacgaTGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACcg                                                       <  1:975027/140‑1 (MQ=255)                                                                                gatgatTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGa                                                    <  1:338780/140‑1 (MQ=255)                                                                                             aTTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCa                                       <  1:246818/140‑1 (MQ=255)                                                                                             aTTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCa                                       >  1:404851/1‑140 (MQ=255)                                                                                             aTTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCa                                       >  2:246818/1‑140 (MQ=255)                                                                                               tGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACt                                     >  2:8093/1‑140 (MQ=255)                                                                                                     ctcGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCa                               >  1:52734/1‑140 (MQ=255)                                                                                                       cGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATc                             <  1:848833/140‑1 (MQ=255)                                                                                                       cGCAGCAACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATc                             >  1:138431/1‑140 (MQ=255)                                                                                                         cagcaACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCg                           <  2:292977/140‑1 (MQ=255)                                                                                                         cagcaACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCg                           <  1:793409/140‑1 (MQ=255)                                                                                                          agcaACCATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCgg                          <  1:8093/140‑1 (MQ=255)                                                                                                               ccATCAACGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCaa                     >  2:453496/1‑140 (MQ=255)                                                                                                                    aaCGCTTCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATAc                <  1:68554/140‑1 (MQ=255)                                                                                                                         ttCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGccc           >  1:109425/1‑140 (MQ=255)                                                                                                                          tCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATAc                <  2:111742/134‑1 (MQ=255)                                                                                                                          tCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATAc                >  1:111742/1‑134 (MQ=255)                                                                                                                          tCGATTTCTTCCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCt          >  1:983327/1‑140 (MQ=255)                                                                                                                                tcttcCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGctct    <  2:52734/140‑1 (MQ=255)                                                                                                                                tcttcCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGctct    <  1:453496/140‑1 (MQ=255)                                                                                                                                tcttcCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGctct    <  1:1002132/140‑1 (MQ=255)                                                                                                                                  ttcCGGC‑TTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTAc  >  1:601799/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      TCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTAC  >  NC_000913/3053425‑3053694

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: ‑