Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,928,699	A→G	G47G (GGT→GGC)	ptrB ←	oligopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,928,699	0	A	G	95.1%	127.8 / ‑2.0	41	G47G (GGT→GGC)	ptrB	oligopeptidase B
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (17/22);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTT  >  NC_000913/1928563‑1928829                                                                                                                                         |                                                                                                                                   aGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCg                                                                                                                                 <  2:794834/140‑1 (MQ=255)  gCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGt                                                                                                                                >  2:688584/1‑140 (MQ=255)  gCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGt                                                                                                                                >  2:661736/1‑140 (MQ=255)  gCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGt                                                                                                                                <  1:810718/140‑1 (MQ=255)  gCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGt                                                                                                                                >  2:163109/1‑140 (MQ=255)    cTGGTTCATAAATATGCCGATACCGGTAGCCTTTTCTGATGTAGGGCGTAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTaa                                                                                                                              <  2:299206/140‑1 (MQ=255)      ggTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACt                                                                                                                            >  2:23642/1‑140 (MQ=255)       gTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTa                                                                                                                           <  1:213741/140‑1 (MQ=255)             aaaTATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCt                                                                                                                     >  2:485803/1‑140 (MQ=255)                 atGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCttgtt                                                                                                                 <  1:688584/140‑1 (MQ=255)                 atGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCATTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCttgtt                                                                                                                 <  2:124834/140‑1 (MQ=255)                     cGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCttgtt                                                                                                                 <  2:757653/136‑1 (MQ=255)                     cGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCttgtt                                                                                                                 >  1:757653/1‑136 (MQ=255)                     cGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAg                                                                                                             <  2:508508/140‑1 (MQ=255)                      gATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAgg                                                                                                            <  1:23642/140‑1 (MQ=255)                      gATAGCGGTAGCCATTTGTGATGTAGGGCGCCGAAGCTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAgg                                                                                                            <  2:320147/140‑1 (MQ=255)                           cGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTc                                                                                                       >  2:607316/1‑140 (MQ=255)                             gTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCa                                                                                                     >  2:168524/1‑140 (MQ=255)                             gTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCa                                                                                                     <  1:163109/140‑1 (MQ=255)                                     ttttGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCt                                                                                             >  1:647311/1‑140 (MQ=255)                                               gggCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAAc                                                                                   >  2:11584/1‑140 (MQ=255)                                                            tctcGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCg                                                                      <  2:215322/140‑1 (MQ=255)                                                             ctcGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGc                                                                     >  2:659488/1‑140 (MQ=255)                                                              tcGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCa                                                                    <  1:670191/140‑1 (MQ=255)                                                                gTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGc                                                                  <  2:164907/140‑1 (MQ=255)                                                                               cGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCg                                                   <  1:168524/140‑1 (MQ=255)                                                                                        ccTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTa                                          <  1:485803/140‑1 (MQ=255)                                                                                                aTGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATg                                  <  2:647311/140‑1 (MQ=255)                                                                                                  gCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGaa                                <  1:11584/140‑1 (MQ=255)                                                                                                  gCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGaa                                <  1:659488/140‑1 (MQ=255)                                                                                                   cGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAg                               >  2:507620/1‑140 (MQ=255)                                                                                                          gCAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATg                        >  2:445428/1‑140 (MQ=255)                                                                                                           cAAGGCTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATgg                       <  1:445428/140‑1 (MQ=255)                                                                                                                cTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTACTTAACGATGCGCGTATAGCCATGAAGCGTCATGGCGTgg                  >  2:387333/1‑140 (MQ=255)                                                                                                                cTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGAATCGCCATGTAGCGTCATGGCGTgg                  >  1:397000/1‑140 (MQ=255)                                                                                                                cTTGTTGTGAGGCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATGATCGATGCGCGTATGGCCATGAAGCGTCATGGCGTgg                  >  1:329483/1‑140 (MQ=255)                                                                                                                          ggCCATCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcggg        <  1:661736/140‑1 (MQ=255)                                                                                                                              aTCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGggcggc    >  2:677053/1‑140 (MQ=255)                                                                                                                               tCACCCGATGGCCGTTACTATTGTCTTGTCGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCt   <  1:387333/140‑1 (MQ=255)                                                                                                                               tCACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCt   <  1:607316/140‑1 (MQ=255)                                                                                                                                cACCCGATGGCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGGCTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtt  >  2:573613/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                   AGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTT  >  NC_000913/1928563‑1928829

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