Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	904,251	G→A	9.1%	L169L (CTG→CTA)	PP_0788 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	904,251	0	G	A	9.1%	120.6 / 3.2	44	L169L (CTG→CTA)	PP_0788	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/24);  new base A (2/2);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

ATGGAACGGGCCGCCGCCATCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTATGTGCTGCGCGATGTCGATGTGATCTATGGCGAAGCCAACCAGTTGCTGGCCGCTG  >  NC_002947/904105‑904390                                                                                                                                                   |                                                                                                                                            aTGGAACGGGCCGCCGCCATCGAGGGCATGGTGGCGCGGGGGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTgtgt                                                                                                                                          >  1:449426/1‑150 (MQ=255)  tGGAACGGGCCGCCGCCATCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTc                                                                                                                                         <  2:342833/150‑1 (MQ=255)  tGGAACGGGCCGCCGCCATCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTc                                                                                                                                         >  1:342833/1‑150 (MQ=255)  tGGAACGGGCCGCCGCCATCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTc                                                                                                                                         <  2:45133/150‑1 (MQ=255)  tGGAACGGGCCGCCGCCATCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTc                                                                                                                                         >  1:450522/1‑150 (MQ=255)           ccgccgccATCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTcc                                                                                                                                <  1:250922/150‑1 (MQ=255)            cgccgccATCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCt                                                                                                                               <  2:593692/150‑1 (MQ=255)                   aTCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGcc                                                                                                                        >  1:543587/1‑150 (MQ=255)                     cGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCgg                                                                                                                      <  1:107714/150‑1 (MQ=255)                     cGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCgg                                                                                                                      <  2:331494/150‑1 (MQ=255)                            aTGGTGGCGCGGGTGCTGTCACGCGATGGCGTCTCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAAt                                                                                                               <  2:344355/150‑1 (MQ=255)                                      gggTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTg                                                                                                     <  2:51949/150‑1 (MQ=255)                                      gggTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTATCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTg                                                                                                     <  1:119483/150‑1 (MQ=255)                                           cTGTCACTCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTATGTCGTTACTTCCTCGACGCCGGGCGTAATCTggg                                                                                                          <  2:417787/140‑1 (MQ=255)                                           cTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTATGTCGTTACTTCCTCGACGCCGGGCGTAATCTggg                                                                                                          >  1:417787/1‑140 (MQ=255)                                           cTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTATGTCGTTACTTCCTCGACGCCGGGCGTAATCTggg                                                                                                          <  2:41776/140‑1 (MQ=255)                                           cTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTATGTCGTTACTTCCTCGACGCCGGGCGTAATCTggg                                                                                                          >  1:41776/1‑140 (MQ=255)                                                aCGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGt                                                                                           >  1:554204/1‑150 (MQ=255)                                                 cgcgATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTc                                                                                          <  1:233085/150‑1 (MQ=255)                                                     aTGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGATGCTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGc                                                                                      >  1:235894/1‑150 (MQ=255)                                                       ggCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTg                                                                                    <  2:154944/150‑1 (MQ=255)                                                              cccTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGtgctgc                                                                             >  2:599270/1‑150 (MQ=255)                                                              cccTGCAGCAACTGGGTCTGCGCCTAAAGGCATTGCTCGACCGACCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCGTCGAGTGCGTTCGTAGTCTGGGCGTTGGCCTGCGGGTCAAGCTGGtgctgc                                                                             <  1:693932/150‑1 (MQ=255)                                                               ccTGCAGCAACTGGGTCTGCGCCTGAAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTgctgct                                                                            <  1:543687/150‑1 (MQ=255)                                                                cTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCg                                                                                              <  1:184368/131‑1 (MQ=255)                                                                cTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCg                                                                                              >  2:184368/1‑131 (MQ=255)                                                                   cagcaACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAg                                                                        <  2:554204/150‑1 (MQ=255)                                                                     gcaACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCt                                                                      >  1:251690/1‑150 (MQ=255)                                                                        aCTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGtgctgc                                                                             >  2:550598/1‑140 (MQ=255)                                                                        aCTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGtgctgc                                                                             <  1:550598/140‑1 (MQ=255)                                                                             gTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGc                                                              <  1:668983/150‑1 (MQ=255)                                                                             gTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGc                                                              <  2:58716/150‑1 (MQ=255)                                                                             gTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGc                                                              <  2:543587/150‑1 (MQ=255)                                                                                     cTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTAtgtg                                                      <  1:37284/150‑1 (MQ=255)                                                                                      tGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCg                                                                                                                       <  1:451915/84‑1 (MQ=255)                                                                                      tGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCg                                                                                                                       >  2:451915/1‑84 (MQ=255)                                                                                                  tCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTATGTGCTGCGcgatgtcg                                         >  1:514074/1‑150 (MQ=255)                                                                                                     aCCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTAtgt                                                       >  2:362522/1‑133 (MQ=255)                                                                                                     aCCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTAtgt                                                       <  1:362522/133‑1 (MQ=255)                                                                                                         cccccTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTATGTGCTGCGCGATGTCGATGTGAt                                  <  2:450522/150‑1 (MQ=255)                                                                                                                   gAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGc                                                                                                                    <  2:584810/58‑1 (MQ=255)                                                                                                                   gAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGc                                                                                                                    >  1:584810/1‑58 (MQ=255)                                                                                                                                 aTTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTATGTGCTGCGCGATGTCGATGTGATCTATGGCGAAGCCAACCAGTTGCt          >  1:59619/1‑150 (MQ=255)                                                                                                                                    ggggCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTATGTGCTGCGCGATGTCGATGTGATCTATGGCGAAGCCAACCAGTTGCTggg       >  2:529369/1‑149 (MQ=255)                                                                                                                                         cAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTATGTGCTGCGCGATGTCGATGTGATCTATGGCGAAGCCAACCAGTTGCTGGCCGCTg  >  1:68132/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                            ATGGAACGGGCCGCCGCCATCGAGGGCATGGTGGCGCGGGTGCTGTCACGCGATGGCGTCGCCCTGCAGCAACTGGGTCTGCGCCTGCAGGCATTGCTCGACCGCCCCCTGTACGAGCAGCACAACCCATTGGGGCCAGCCGCGCTGTGTCGTTACTTCCTCGACGCCGGGCGTAATCTGGGCGTTGGCCTGCGGGTCAAGCTGGTGCTGCTGAAGCTGTTCGAGCGTTATGTGCTGCGCGATGTCGATGTGATCTATGGCGAAGCCAACCAGTTGCTGGCCGCTG  >  NC_002947/904105‑904390

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