Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,958,537	A→G	G116G (GGT→GGC)	ptrA ←	protease 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,958,537	0	A	G	100.0%	125.1 / NA	38	G116G (GGT→GGC)	ptrA	protease 3
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/18);  total (20/18)

TCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCA  >  NC_000913/2958406‑2958670                                                                                                                                    |                                                                                                                                      tCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTg                                                                                                                               <  1:138744/140‑1 (MQ=255) tCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTg                                                                                                                               <  1:525175/140‑1 (MQ=255)   ggCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCa                                                                                                                             >  1:2072/1‑140 (MQ=255)         tttCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTga                                                                                                                       >  2:473110/1‑140 (MQ=255)               gTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATAtt                                                                                                                 >  2:344525/1‑140 (MQ=255)               gTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTCCTGGCATAGTGGCTGCCGCCGCGCATTTTGAGATAtt                                                                                                                 >  1:220335/1‑140 (MQ=255)                  gAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCgg                                                                                                              >  1:562219/1‑140 (MQ=255)                                cAATAGCATCGGCCAAGCGGTCTACAGCATCAGGCAACGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGcc                                                                                                <  1:170417/140‑1 (MQ=255)                                           gCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACtt                                                                                     >  1:82874/1‑140 (MQ=255)                                                  ggTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGAc                                                                              <  1:495557/140‑1 (MQ=255)                                                   gTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGAcc                                                                             <  1:562106/140‑1 (MQ=255)                                                    tCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGAccc                                                                            >  2:115078/1‑140 (MQ=255)                                                       aCCGCACCAGGCAAGGCTTCGTTCTCAACTTTCAGCTAGTAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCAt                                                                         <  2:183481/140‑1 (MQ=255)                                                        ccGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTACATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATc                                                                        <  2:365597/140‑1 (MQ=255)                                                           cACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGa                                                                     <  2:107608/140‑1 (MQ=255)                                                              cAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTc                                                                  >  2:551585/1‑140 (MQ=255)                                                              cAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTc                                                                  >  1:57917/1‑140 (MQ=255)                                                                ggCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCat                                                                <  1:344525/140‑1 (MQ=255)                                                                    aGGCGTCGTTCTCAACTTCCAGCTAGAAAGCTGTGCGATTCGGCGCAGGGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGt                                                            <  2:220335/140‑1 (MQ=255)                                                                     ggCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTCCTTATTAGACCCCATCAGACTCATATGtt                                                           >  1:214365/1‑140 (MQ=255)                                                                            ttCTCAACTTCCAGATAAAAAGGCGTGCGATGCGGTGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTa                                                    <  2:244097/140‑1 (MQ=255)                                                                               tcAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGTCTGCGGGTACTTGTTCGACCCCATCAGACTGATATGTTCAAGGTAAtg                                                 >  2:195775/1‑140 (MQ=255)                                                                                    ttCCAGATAGAAAGCTGTGCGATACCGCGCAGTGCTGCCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCa                                            <  1:195775/140‑1 (MQ=255)                                                                                      ccAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCCGACTCATATGTTCAAGGTAATGTGCCAGc                                          >  1:244097/1‑140 (MQ=255)                                                                                       cAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGcc                                         >  2:593129/1‑140 (MQ=255)                                                                                              aaaGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATGTGTGCAGGGTAATGTGCCAGCCCCTGGTa                                  >  1:187363/1‑140 (MQ=255)                                                                                              aaaGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTa                                  >  1:555959/1‑140 (MQ=255)                                                                                              aaaGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTa                                  <  2:53480/140‑1 (MQ=255)                                                                                              aaaGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTa                                  >  1:53480/1‑140 (MQ=255)                                                                                              aaaGCCGTGCGATACGGCGCAGCGCTGGCATCGTGCCTGCCGCCGTGCCTTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTa                                  <  2:187363/140‑1 (MQ=255)                                                                                                  ccGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGcc                              <  1:593129/140‑1 (MQ=255)                                                                                                   cGTGCGATACGGCGCAGTGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCt                             <  2:645869/140‑1 (MQ=255)                                                                                                                    tGCTGGTATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGAc            >  1:489398/1‑140 (MQ=255)                                                                                                                    tGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGAc            >  2:460391/1‑140 (MQ=255)                                                                                                                    tGCTGGCATTGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGAc            <  2:344126/140‑1 (MQ=255)                                                                                                                            ttGTGACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACggg    <  2:489398/140‑1 (MQ=255)                                                                                                                              gtgACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGca  >  1:468551/1‑140 (MQ=255)                                                                                                                              gtgACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGca  >  2:555861/1‑140 (MQ=255)                                                                                                                              gtgACTGCCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGca  >  1:26095/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                      TCGGCATATTTCTTGTCGAGCAAAGGTTCAGCAATAGCATCGGCCAGGCGGTCTACCGCACCAGGCAAGGCGTCGTTCTCAACTTCCAGATAGAAAGCCGTGCGATACGGCGCAGTGCTGGCATTGTGACTACCGCCGTGCATTTTGAGATATTCGGCCAGACTGTCAGCCTGCGGGTACTTTTTCGACCCCATCAGACTCATATGTTCAAGGTAATGTGCCAGCCCCTGGTACGCCTCGGGATCTTCCAGCGACCCAACGGGCA  >  NC_000913/2958406‑2958670

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