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

GCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACC  >  NC_000913/2655632‑2655893                                                                                                                                        |                                                                                                                               gCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGAt                                                                                                                             >  1:1358267/1‑139 (MQ=255)    tCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATAcc                                                                                                                          >  2:961797/1‑139 (MQ=255)       cgcgCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCttt                                                                                                                       >  1:967917/1‑139 (MQ=255)        gcgcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGACAAAAGCGATACCTTta                                                                                                                      <  1:961797/139‑1 (MQ=255)         cgcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtac                                                                                                                     <  2:981190/139‑1 (MQ=255)            ttGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtaccta                                                                                                                  >  1:781157/1‑139 (MQ=255)                 cAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                             <  1:710236/139‑1 (MQ=255)                               tGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCGATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTgg                                                                                               <  2:967917/139‑1 (MQ=255)                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                      >  1:1371742/1‑99 (MQ=255)                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                      <  2:1371742/99‑1 (MQ=255)                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAg                                                                              <  1:96912/139‑1 (MQ=255)                                                                               aCTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCgg                                               >  2:104820/1‑139 (MQ=255)                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATAc                                                                                                                           <  1:12441/58‑1 (MQ=255)                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATAc                                                                                                                           >  2:12441/1‑58 (MQ=255)                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAgc                                                             <  1:1225546/120‑1 (MQ=255)                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAgc                                                             >  2:1225546/1‑120 (MQ=255)                                                                                      cGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATc                                                                                   <  1:101213/96‑1 (MQ=255)                                                                                      cGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATc                                                                                   >  2:101213/1‑96 (MQ=255)                                                                                            cATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGaa                                  <  1:470753/139‑1 (MQ=255)                                                                                            cATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGaa                                  <  2:1358267/139‑1 (MQ=255)                                                                                              tAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAAcc                                <  1:104820/139‑1 (MQ=255)                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                             >  2:1495273/1‑54 (MQ=255)                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                             <  1:1495273/54‑1 (MQ=255)                                                                                                             tGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAAcc                                >  1:704303/1‑124 (MQ=255)                                                                                                             tGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAAcc                                <  2:704303/124‑1 (MQ=255)                                                                                                             tGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACtgtgtg                 <  2:1436407/139‑1 (MQ=255)                                                                                                             tGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGACACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACtgtgtg                 <  2:194783/139‑1 (MQ=255)                                                                                                                tGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                             >  1:881441/1‑44 (MQ=255)                                                                                                                tGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                             <  2:881441/44‑1 (MQ=255)                                                                                                                 gCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTaca                                                                                                   >  2:1113832/1‑53 (MQ=255)                                                                                                                 gCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTaca                                                                                                   <  1:1113832/53‑1 (MQ=255)                                                                                                                   tGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGccc           >  1:172641/1‑139 (MQ=255)                                                                                                                   tGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGccc           >  2:1018783/1‑139 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGcc            <  2:1112886/136‑1 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGcc            >  1:1112886/1‑136 (MQ=255)                                                                                                                            cgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTac                                                                                                    <  2:954195/41‑1 (MQ=255)                                                                                                                            cgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTac                                                                                                    >  1:954195/1‑41 (MQ=255)                                                                                                                            cgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAAcc  >  2:535397/1‑139 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGct                                                                                           <  1:1382595/47‑1 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGct                                                                                           >  2:1382595/1‑47 (MQ=255)                                                                                                                                   aaaaGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTg                                                                           >  2:327881/1‑59 (MQ=255)                                                                                                                                   aaaaGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTg                                                                           <  1:327881/59‑1 (MQ=255)                                                                                                                                        |                                                                                                                               GCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACC  >  NC_000913/2655632‑2655893

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