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

GCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCTCA  >  NC_000913/2655632‑2655900                                                                                                                                        |                                                                                                                                      gCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGAt                                                                                                                                    >  1:1121501/1‑139 (MQ=255)  cTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCACGTCCGACTTCATCGAGTCCATAAACGCAGTATGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATa                                                                                                                                   <  1:1050705/139‑1 (MQ=255)      aCGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCtt                                                                                                                               <  2:142653/139‑1 (MQ=255)        gcgcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                             <  1:1695928/139‑1 (MQ=255)          gcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtacc                                                                                                                           <  2:302341/139‑1 (MQ=255)            ttGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtaccta                                                                                                                         <  2:671220/139‑1 (MQ=255)             tGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtacctac                                                                                                                        >  1:640878/1‑139 (MQ=255)                 cAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    <  2:511569/139‑1 (MQ=255)                 cAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    <  2:352490/139‑1 (MQ=255)                 cAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    <  1:928214/139‑1 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                             >  1:1448480/1‑123 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                             <  2:1448480/123‑1 (MQ=255)                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGt                                                                                                             >  1:682433/1‑139 (MQ=255)                                   aaaTGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGct                                                                                                  <  2:218524/139‑1 (MQ=255)                                       gCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGAt                                                                                                                                    >  1:1737161/1‑101 (MQ=255)                                       gCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGAt                                                                                                                                    <  2:1737161/101‑1 (MQ=255)                                          aGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTAt                                                                                           >  2:223503/1‑139 (MQ=255)                                           gcgcCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTacc                                                                                          >  1:1319250/1‑137 (MQ=255)                                                 cGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGt                                                                                    >  1:629596/1‑139 (MQ=255)                                                   gTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATTGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTg                                                                                  <  2:1211478/139‑1 (MQ=255)                                                   gTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTACCGCCAGTTg                                                                                  <  2:1319250/139‑1 (MQ=255)                                                             ccgcgccgcCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATc                                                                        >  1:1473896/1‑139 (MQ=255)                                                                  ccgccCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAg                                                                                     >  2:946082/1‑121 (MQ=255)                                                                  ccgccCATGTCCGACTTCACCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAg                                                                                     <  1:946082/121‑1 (MQ=255)                                                                     cccATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCa                                                                <  2:563654/139‑1 (MQ=255)                                                                         tGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCaa                                                            <  2:640878/139‑1 (MQ=255)                                                                              gACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCg                                                       <  2:1750483/139‑1 (MQ=255)                                                                               aCTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCgg                                                      >  1:757606/1‑139 (MQ=255)                                                                               aCTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCgg                                                      >  1:365621/1‑139 (MQ=255)                                                                                cTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGga                                                     >  1:1358876/1‑139 (MQ=255)                                                                                  tCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGgaga                                                   <  1:657151/139‑1 (MQ=255)                                                                                           ccATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGtt                                                                                   >  2:457850/1‑98 (MQ=255)                                                                                           ccATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGtt                                                                                   <  1:457850/98‑1 (MQ=255)                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGa                               <  1:223503/139‑1 (MQ=255)                                                                                                       tCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGAc                              >  2:1530734/1‑139 (MQ=255)                                                                                                        cTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTaca                                                                                                          <  2:1236256/62‑1 (MQ=255)                                                                                                        cTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTaca                                                                                                          >  1:1236256/1‑62 (MQ=255)                                                                                                        cTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACt                             >  1:1131242/1‑139 (MQ=255)                                                                                                             tGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    <  1:806844/47‑1 (MQ=255)                                                                                                             tGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    >  2:806844/1‑47 (MQ=255)                                                                                                             tGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACtgtgtg                        <  2:629596/139‑1 (MQ=255)                                                                                                                tGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAg                     <  2:682433/139‑1 (MQ=255)                                                                                                                 gCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGc                    >  2:1027873/1‑139 (MQ=255)                                                                                                                  cTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGcc                   <  2:1358876/139‑1 (MQ=255)                                                                                                                   tGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGgaga                                                   >  2:1445891/1‑106 (MQ=255)                                                                                                                   tGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGgaga                                                   <  1:1445891/106‑1 (MQ=255)                                                                                                                   tGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGccc                  >  2:1077741/1‑139 (MQ=255)                                                                                                                   tGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGccc                  >  2:271808/1‑139 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTa                                                                                            >  1:1247014/1‑63 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTa                                                                                            <  2:1247014/63‑1 (MQ=255)                                                                                                                     tAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCa                >  2:664960/1‑139 (MQ=255)                                                                                                                          gccgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACtgtgtg                        >  2:515901/1‑126 (MQ=255)                                                                                                                          gccgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACtgtgtg                        <  1:515901/126‑1 (MQ=255)                                                                                                                            cgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAgcg                                                                   <  1:764045/81‑1 (MQ=255)                                                                                                                            cgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAgcg                                                                   >  2:764045/1‑81 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGcc                   >  2:1018870/1‑126 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGcc                   <  1:1018870/126‑1 (MQ=255)                                                                                                                               aGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTg      >  1:465709/1‑139 (MQ=255)                                                                                                                                   aaaaGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCTCa  >  1:1589977/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                      GCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCTCA  >  NC_000913/2655632‑2655900

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