Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,423,158	G→A	8.9%	A83V (GCG→GTG)	PP_4763 ←	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,423,158	0	G	A	8.9%	124.4 / 2.8	45	A83V (GCG→GTG)	PP_4763	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (24/17);  new base A (2/2);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAGGCAGGTTGCATTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGCCGATCAACGCAGCCGCCTCAAC  >  NC_002947/5423013‑5423306                                                                                                                                                  |                                                                                                                                                     cGAAGGCAGGTTGCATTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCAccgccg                                                                                                                                                   <  5:277570/149‑1 (MQ=255) cGAAGGCAGGTTGCATTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCAccgccg                                                                                                                                                   <  7:223170/149‑1 (MQ=255)  gAAGGCAGGTTGCATTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCgg                                                                                                                                                  <  8:26335/149‑1 (MQ=255)    aGGCAGGTTGCATTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCgggg                                                                                                                                                <  5:19727/149‑1 (MQ=255)        aGGTTGCATTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCaaa                                                                                                                                            <  7:305408/149‑1 (MQ=255)           ttGCATTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAgg                                                                                                                                         <  6:307971/149‑1 (MQ=255)               aTTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGcgcg                                                                                                                                     >  8:289884/1‑149 (MQ=255)                    cTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCg                                                                                                                                  >  4:144581/1‑147 (MQ=255)                      ggTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCAccgccg                                                                                                                                                   <  6:193722/128‑1 (MQ=255)                      ggTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCAccgccg                                                                                                                                                   >  5:193722/1‑128 (MQ=255)                       gTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGtt                                                                                                                             >  3:159319/1‑149 (MQ=255)                       gTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGt                                                                                                                              <  1:282147/148‑1 (MQ=255)                             gacacgacCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGCTCTCCCCCGCCGGGGCAAAAGGCGCGGCGAAGtggcgagt                                                                                                                       >  2:299497/1‑142 (MQ=255)                                     ctcctcCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAAcc                                                                                                               >  2:280590/1‑149 (MQ=255)                                       cctccAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCg                                                                                                             >  3:101011/1‑149 (MQ=255)                                                gcAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGgc                                                                                                                                        <  8:141164/113‑1 (MQ=255)                                                gcAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGgc                                                                                                                                        >  7:141164/1‑113 (MQ=255)                                                 cAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGa                                                                                                   <  4:159319/149‑1 (MQ=255)                                                                acaACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGt                                                                                    >  3:164594/1‑149 (MQ=255)                                                                 caACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTc                                                                                   >  4:57243/1‑149 (MQ=255)                                                                 caACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTc                                                                                   <  3:278733/149‑1 (MQ=255)                                                                   aCGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCt                                                                                 >  3:123865/1‑149 (MQ=255)                                                                        tgcaggTCCACGCCACTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTgcgc                                                                            <  2:171326/148‑1 (MQ=255)                                                                         gcaggTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCt                                                                           >  8:29780/1‑149 (MQ=255)                                                                         ccgctctcACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCACCGGGGCAAAAGGCGCGGCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTgcgc                                                                            <  5:45572/141‑1 (MQ=255)                                                                             gTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAg                                                                       >  6:97127/1‑149 (MQ=255)                                                                                 aCGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGcc                                                                   >  5:40938/1‑149 (MQ=255)                                                                                 aCGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCACCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCCTGCCGATGAACGCGCCACTGTCCTTgcgc                                                                            >  6:50329/1‑140 (MQ=255)                                                                                 aCGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCACCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTgcgc                                                                            <  5:50329/140‑1 (MQ=255)                                                                                 aCGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGAAGGCGTGCCAGACGCCAGCCGATCTCCACCACCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTgcgc                                                                            >  6:45572/1‑140 (MQ=255)                                                                                  cgccgcTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGccc                                                                  <  7:29780/149‑1 (MQ=255)                                                                                       cTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACaaa                                                             <  3:144581/149‑1 (MQ=255)                                                                                       cTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACaaa                                                             >  1:178928/1‑149 (MQ=255)                                                                                           gCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCaa                                                         <  3:57243/149‑1 (MQ=255)                                                                                               gCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAAc                                                                                              >  4:269318/1‑108 (MQ=255)                                                                                               gCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAAc                                                                                              <  3:269318/108‑1 (MQ=255)                                                                                                     cccAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACt                                                >  1:265758/1‑148 (MQ=255)                                                                                                     cccAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTc                                               <  4:101011/149‑1 (MQ=255)                                                                                                            cccAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGaaa                                        >  5:283765/1‑149 (MQ=255)                                                                                                             ccAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAAt                                       <  6:40938/149‑1 (MQ=255)                                                                                                               aGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCAt                                             <  4:151585/141‑1 (MQ=255)                                                                                                               aGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCAt                                             >  3:151585/1‑141 (MQ=255)                                                                                                                          gCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCAcgcg                          >  8:96622/1‑149 (MQ=255)                                                                                                                            cAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGAAGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGc                        >  2:120613/1‑149 (MQ=255)                                                                                                                                    aGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCGCAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGCCGATCAAc                >  8:193986/1‑149 (MQ=255)                                                                                                                                      ccGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGCCGATCAACg               >  4:45991/1‑148 (MQ=255)                                                                                                                                         aTCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGCCGATCAACGCAgc           >  3:259108/1‑149 (MQ=255)                                                                                                                                           ctcCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGCCGATCAACGCAgccg         >  4:156379/1‑149 (MQ=255)                                                                                                                                              cACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGCCGATCAACGCAGCCGCCt      <  4:123865/149‑1 (MQ=255)                                                                                                                                                  gccgGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGCCGATCAACGCAGCCGCCTCAAc  >  7:151884/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CGAAGGCAGGTTGCATTCGCTGGTGAACGACACGACCTCCTCCAGGCGCAATTGGGCGAAAGCACAACGCAGGCAGGTCCACGCCGCTTCGCTGGCAAAGCCCAGGCCCCAGTGACGGCGTGCCAGACGCCAGCCGATCTCCACCGCCGGGGCAAAAGGCGCGTCGAAGTTGACATTGAGCAACCCGGTCATGCCGATGAACGCGCCACTGTCCTTGCGCTCCAGGGCCCACAAACCAAAACCGTACTCATTGAAATGCCCACGCACGCGGCCGATCAACGCAGCCGCCTCAAC  >  NC_002947/5423013‑5423306

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