Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,957,284	A→G	V534A (GTG→GCG)	ptrA ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,957,284	0	A	G	100.0%	159.1 / NA	46	V534A (GTG→GCG)	ptrA	protease III
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/28);  total (18/28)

AACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAG  >  NC_000913/2957137‑2957431                                                                                                                                                    |                                                                                                                                                    aaCTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCgcgc                                                                                                                                                  <  1:581844/151‑1 (MQ=255)    tGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCgcgcgc                                                                                                                                                >  2:328064/1‑150 (MQ=255)                aGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACtcat                                                                                                                                   <  1:239913/151‑1 (MQ=255)                 gCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACtcatc                                                                                                                                  <  1:89823/151‑1 (MQ=255)                      gCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATcaacaa                                                                                                                             <  2:279685/151‑1 (MQ=255)                              aTCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGctct                                                                                                                     >  1:547932/1‑151 (MQ=255)                                cATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTgg                                                                                                                   >  2:27208/1‑151 (MQ=255)                                cATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTgg                                                                                                                   <  2:189175/151‑1 (MQ=255)                                      gcgcAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTc                                                                                                             <  2:219142/151‑1 (MQ=255)                                        gcAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGt                                                                                                           >  2:196432/1‑151 (MQ=255)                                         cAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGt                                                                                                           >  2:386405/1‑150 (MQ=255)                                              aTCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATttct                                                                                                     >  2:32742/1‑151 (MQ=255)                                              aTCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATttct                                                                                                     <  1:612542/151‑1 (MQ=255)                                                     gATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTg                                                                                               <  1:208418/150‑1 (MQ=255)                                                           gggCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAAt                                                                                        <  2:683232/151‑1 (MQ=255)                                                              cgcTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAg                                                                                     <  2:463009/151‑1 (MQ=255)                                                               gcTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGc                                                                                    <  2:224704/151‑1 (MQ=255)                                                                    ccATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCgaga                                                                                >  2:414650/1‑150 (MQ=255)                                                                         gCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAatcat                                                                          <  2:636285/151‑1 (MQ=255)                                                                            ttCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAAtcatca                                                                        <  2:301449/150‑1 (MQ=255)                                                                              cGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGa                                                                     <  2:417764/151‑1 (MQ=255)                                                                               ggATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGaa                                                                    <  1:546687/151‑1 (MQ=255)                                                                                      gCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAgg                                                             >  1:407319/1‑151 (MQ=255)                                                                                        aaaaTCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGGATATAAggg                                                            >  1:370289/1‑150 (MQ=255)                                                                                             cAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGc                                                      >  1:609824/1‑151 (MQ=255)                                                                                              aGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGct                                                     <  2:644863/151‑1 (MQ=255)                                                                                                 cTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTg                                                  >  1:547472/1‑151 (MQ=255)                                                                                                  tGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTgg                                                 >  2:358423/1‑151 (MQ=255)                                                                                                   gACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATAAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTg                                                  <  1:8021/149‑1 (MQ=255)                                                                                                     cATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCaa                                              <  2:675734/151‑1 (MQ=255)                                                                                                      aTCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCaaa                                             <  1:137383/151‑1 (MQ=255)                                                                                                             ttGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAagagag                                        <  2:385829/149‑1 (MQ=255)                                                                                                               gggCTCGCTGGCAAAATAACGGCTTGGCGCATACCCCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAgcgc                                     <  1:386392/150‑1 (MQ=255)                                                                                                               gggCTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCa                                    <  2:218987/151‑1 (MQ=255)                                                                                                                  cTCGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAAt                                  >  1:121661/1‑150 (MQ=255)                                                                                                                    cGCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATg                                 <  1:470559/149‑1 (MQ=255)                                                                                                                     gCTGGCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTcg                              >  2:530743/1‑151 (MQ=255)                                                                                                                        ggCAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTcggc                            <  1:358412/150‑1 (MQ=255)                                                                                                                          cAAAATAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCggcggc                         >  1:637895/1‑151 (MQ=255)                                                                                                                           aaaaTAACGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCt                        <  1:51599/151‑1 (MQ=255)                                                                                                                            aaaTAACGGCTTGGAGCATACACCGCGAGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCtt                       >  2:584488/1‑151 (MQ=255)                                                                                                                                  cGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCt                 <  2:547958/151‑1 (MQ=255)                                                                                                                                  cGGCTTGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCt                 >  2:628696/1‑151 (MQ=255)                                                                                                                                      ttGGCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGcc              >  1:576457/1‑150 (MQ=255)                                                                                                                                        ggCGCATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCa             <  1:533499/149‑1 (MQ=255)                                                                                                                                           gcATACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCgg        <  2:513104/151‑1 (MQ=255)                                                                                                                                             aTACACCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCg      <  1:480098/151‑1 (MQ=255)                                                                                                                                               acacCGCGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCg      >  2:346694/1‑149 (MQ=255)                                                                                                                                                cacctcGCGCAGATTCGACTCATCAACAAGCAGCTCTGGATGGTCGTAGTTCTTCGCTGACATAATCAGCGAGAAAGCTTCAGGAATTTAAGGGTTCAGCTCTGGCAAAGAGAGCTCAATGTCGGCTGCTTTTTTCTGCCAGTCGGCGaaa   <  1:98830/151‑1 (MQ=255)                                                                                                                                                 accgcGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAg  <  1:315203/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTACGCAAAATCAGGCTGACATCAGCTTTGGGCTCGCTGGCAAAATAACGGCTTGGCGCATACACCACGCGCAGATTCGACTCATCAACAATCAGCTCTGGATGGTCGTATTTCTTCTCTGACTTAATCAGCGAGAAATCATCAGGAATATAAGGGTTAAGCTCTGGCAAAGAGAGCGCAATGTCGGCGGCTTTTTTCTGCCAGTCGGCGAAAG  >  NC_000913/2957137‑2957431

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