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%	161.0 / NA	45	T198A (ACT→GCT)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/22);  total (23/22)

GCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGC  >  NC_000913/2655632‑2655897                                                                                                                                        |                                                                                                                                   gCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCTCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCTGCCGGAGTCAAAAGCGAt                                                                                                                                 >  2:1083056/1‑139 (MQ=255)       cgcgCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCttt                                                                                                                           >  1:1180358/1‑139 (MQ=255)        gcgcTTGTTCAGTCCGCGCGTAATGTCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                          <  2:1042617/139‑1 (MQ=255)               ttCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCa                                                                                                                   >  2:1155806/1‑139 (MQ=255)                 cAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGTCTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 <  1:1155806/139‑1 (MQ=255)                    tCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCAc                                                                                                              >  2:1476654/1‑139 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCttt                                                                                                                           >  2:392504/1‑122 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCttt                                                                                                                           <  1:392504/122‑1 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                          <  1:806616/123‑1 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGt                                                                                                          >  2:516783/1‑139 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGCCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                          >  2:806616/1‑123 (MQ=255)                           gTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTaca                                                                                                       <  2:1180358/139‑1 (MQ=255)                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtacctac                                                                                                                     <  2:456024/104‑1 (MQ=255)                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtacctac                                                                                                                     >  1:456024/1‑104 (MQ=255)                                                  ggTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGtt                                                                                <  1:1083056/139‑1 (MQ=255)                                                                             cGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATAcc                                                     >  1:326255/1‑139 (MQ=255)                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGc                                                                                     <  2:1261724/100‑1 (MQ=255)                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGc                                                                                     >  1:1261724/1‑100 (MQ=255)                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTg                                                                          <  1:1250937/111‑1 (MQ=255)                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTg                                                                          >  2:1250937/1‑111 (MQ=255)                                                                                           ccATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGa                                       >  2:860590/1‑139 (MQ=255)                                                                                             aTAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGcc                                                                                    >  1:1541379/1‑92 (MQ=255)                                                                                             aTAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGcc                                                                                    <  2:1541379/92‑1 (MQ=255)                                                                                                    cAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 >  2:1131525/1‑56 (MQ=255)                                                                                                    cAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 <  1:1131525/56‑1 (MQ=255)                                                                                                    cAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTcc                              >  2:226229/1‑139 (MQ=255)                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 >  1:37261/1‑54 (MQ=255)                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 <  2:37261/54‑1 (MQ=255)                                                                                                       tCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGgaga                                                <  1:700826/118‑1 (MQ=255)                                                                                                       tCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGgaga                                                >  2:700826/1‑118 (MQ=255)                                                                                                                tGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAg                  <  1:1476654/139‑1 (MQ=255)                                                                                                                tGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAg                  <  2:1455099/139‑1 (MQ=255)                                                                                                                 gCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGc                 <  2:326255/139‑1 (MQ=255)                                                                                                                 gCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACC‑ACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGcc                >  1:1306663/1‑139 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 >  2:5862/1‑39 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 <  1:5862/39‑1 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 <  1:768666/39‑1 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                 >  2:768666/1‑39 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCTAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACtgtgtg                     >  1:1128357/1‑131 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACtgtgtg                     <  2:1128357/131‑1 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGc                 >  2:604181/1‑135 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGc                 <  1:604181/135‑1 (MQ=255)                                                                                                                       gccgccGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGgcg                                                                                                 >  2:1328838/1‑53 (MQ=255)                                                                                                                       gccgccGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGgcg                                                                                                 <  1:1328838/53‑1 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACC‑ACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGc  <  2:1306663/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   GCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGC  >  NC_000913/2655632‑2655897

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