Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,655,767	T→C	T198A (ACT→GCT)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,655,767	0	T	C	100.0%	140.3 / NA	40	T198A (ACT→GCT)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/21);  total (19/21)

GCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCT  >  NC_000913/2655632‑2655898                                                                                                                                        |                                                                                                                                    gCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGAt                                                                                                                                  >  1:1183583/1‑139 (MQ=255)  cTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATa                                                                                                                                 <  1:549323/139‑1 (MQ=255)               ttCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCa                                                                                                                    >  1:451317/1‑139 (MQ=255)                 cAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                  <  2:716814/139‑1 (MQ=255)                       gcgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCAcgcg                                                                                                            <  2:930012/139‑1 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGt                                                                                                           >  1:1188143/1‑139 (MQ=255)                             aaTGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACt                                                                                                      <  2:1188143/139‑1 (MQ=255)                                         cAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACATTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTa                                                                                          >  1:939013/1‑139 (MQ=255)                                          aGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTAt                                                                                         <  2:939013/139‑1 (MQ=255)                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGc                                                                                                                 <  1:539925/109‑1 (MQ=255)                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGc                                                                                                                 >  2:539925/1‑109 (MQ=255)                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAg                                                                                   <  2:1183583/139‑1 (MQ=255)                                                 cGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGt                                                                                  >  1:1050461/1‑139 (MQ=255)                                                            gccgcgccgcCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAAt                                                                       >  2:1412699/1‑139 (MQ=255)                                                             ccgcgccgcCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATc                                                                      >  2:843683/1‑139 (MQ=255)                                                                  ccgccCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAgcg                                                                 <  1:421390/139‑1 (MQ=255)                                                                      ccATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAg                                                             <  1:1412699/139‑1 (MQ=255)                                                                                  tCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                  >  2:654276/1‑74 (MQ=255)                                                                                  tCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                  <  1:654276/74‑1 (MQ=255)                                                                                  tCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGgaga                                                 <  1:636198/139‑1 (MQ=255)                                                                                     tCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACg                                              >  1:1440083/1‑139 (MQ=255)                                                                                          tCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGtt                                                                                 >  2:259212/1‑99 (MQ=255)                                                                                          tCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGtt                                                                                 <  1:259212/99‑1 (MQ=255)                                                                                             aTAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAAc                                      <  1:205994/139‑1 (MQ=255)                                                                                              tAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAAcc                                     >  1:1026639/1‑139 (MQ=255)                                                                                                    cAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCa                                                                                                                    <  2:518688/54‑1 (MQ=255)                                                                                                    cAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCa                                                                                                                    >  1:518688/1‑54 (MQ=255)                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTc                                <  1:313055/136‑1 (MQ=255)                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTc                                >  2:313055/1‑136 (MQ=255)                                                                                                       tCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATa                                                                                                                                 >  1:6892/1‑38 (MQ=255)                                                                                                       tCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATa                                                                                                                                 <  2:6892/38‑1 (MQ=255)                                                                                                       tCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGga                                                   <  1:999379/116‑1 (MQ=255)                                                                                                       tCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGga                                                   >  2:999379/1‑116 (MQ=255)                                                                                                        cTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                  <  2:1361935/52‑1 (MQ=255)                                                                                                        cTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                  >  1:1361935/1‑52 (MQ=255)                                                                                                                tGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAg                   <  1:124344/139‑1 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAgcg                                                                 >  2:1384571/1‑78 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAgcg                                                                 <  1:1384571/78‑1 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGTAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTg    >  1:312985/1‑139 (MQ=255)                                                                                                                                 tCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCt  <  1:1307346/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    GCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCT  >  NC_000913/2655632‑2655898

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