Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,511,070	T→A	10.9%	Q166L (CAG→CTG)	PP_1325 ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,511,070	0	T	A	10.9%	118.3 / 3.7	44	Q166L (CAG→CTG)	PP_1325	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (22/17);  new base A (3/2);  total (25/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

CTCCAGGGTGCCGGCGCTTTTCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTGCGCGCTTGCAGGTCTGCCGGCAGTTCGGCGAGGCGCTGCAGGCTGGGGT  >  NC_002947/1510928‑1511214                                                                                                                                                |                                                                                                                                                  ctCCAGGGTGCCGGCGCTTTTCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGgccgc                                                                                                                                             >  4:221475/1‑149 (MQ=255) ctCCAGGGTGCCGGCGCTTTTCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGgccgc                                                                                                                                             >  3:260592/1‑149 (MQ=255)    cAGGGTGCCGGCGCTTTTCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCaa                                                                                                                                          <  3:302307/149‑1 (MQ=255)          gCCGGCGCTTTTCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCg                                                                                                                                    <  5:317690/149‑1 (MQ=255)                gcTTTTCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCgg                                                                                                                                   >  8:206810/1‑144 (MQ=255)                gcTTTTCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCgg                                                                                                                                   <  7:206810/144‑1 (MQ=255)                  ttttCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCcagc                                                                                                                            >  4:21088/1‑149 (MQ=255)                   tttCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCcag                                                                                                                             >  5:47594/1‑147 (MQ=255)                      cACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATTCGTTGGTGCTGGCGGCCAGGCCGCCag                                                                                                                                          <  1:177912/132‑2 (MQ=255)                      cACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCAGGCCGCCag                                                                                                                                          >  2:177912/1‑130 (MQ=255)                         ggCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATgc                                                                                                                     <  6:47594/149‑1 (MQ=255)                               ggCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCgctgct                                                                                                               <  4:260592/149‑1 (MQ=255)                                         gctggaccagcctgGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCa                                                                                                       <  3:236193/141‑1 (MQ=255)                                         gctggaccagccagGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCa                                                                                                       <  5:13927/141‑1 (MQ=255)                                           ccagccagccagGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCgctgct                                                                                                               >  2:119452/1‑137 (MQ=255)                                            cagccagccagGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAg                                                                                                  >  4:22591/1‑149 (MQ=255)                                              gccagccagGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAggg                                                                                                <  3:221475/149‑1 (MQ=255)                                                    cagGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGc                                                                                          >  5:63278/1‑149 (MQ=255)                                                             aTCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCtt                                                                                >  6:205127/1‑149 (MQ=255)                                                                          tgctgcAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCgcggcggc                                                                   >  5:58206/1‑149 (MQ=255)                                                                            ctgcAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAg                                                                            <  8:121737/138‑1 (MQ=255)                                                                            ctgcAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAg                                                                            >  7:121737/1‑138 (MQ=255)                                                                             tgcAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCatg                                                                                                 <  6:240967/117‑3 (MQ=255)                                                                             tgcAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCatg                                                                                                 >  5:240967/1‑115 (MQ=255)                                                                                aGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGca                                                              >  2:15509/1‑148 (MQ=255)                                                                                    gCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACt                                                          <  6:58206/148‑1 (MQ=255)                                                                                     ctctGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCGCGCTTCGAGCGCGGCGGCGTGCACACTg                                                         <  1:15509/149‑1 (MQ=255)                                                                                      tctGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCAt                                                                                       >  5:83454/1‑118 (MQ=255)                                                                                      tctGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCAt                                                                                       <  6:83454/118‑1 (MQ=255)                                                                                          caggcagCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTg                                                   >  8:160066/1‑149 (MQ=255)                                                                                          caggcagCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTg                                                   >  3:82012/1‑149 (MQ=255)                                                                                                  gAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTCTAGGTGCGCGCTTg                                           >  6:292503/1‑149 (MQ=255)                                                                                                          cAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAg                                                      <  4:101506/130‑1 (MQ=255)                                                                                                          cAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAg                                                      >  3:101506/1‑130 (MQ=255)                                                                                                              gCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTGCGCGCTTGCAGGTCTGCCgg                               <  6:63278/149‑1 (MQ=255)                                                                                                                          tCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTGCGCGCTTGCAGGTCTGCCGGCAGTTCGGCGAg                   <  3:22591/149‑1 (MQ=255)                                                                                                                                gttgGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTGCGCGCTTGCAGGTCTGCCGGCAGTTCGGCGAGGCGCt              >  6:228235/1‑149 (MQ=255)                                                                                                                                 ttgGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTGCGCGCTTGCAGGTCTGCCGGCAGTTCGGCGAGGCGCTg             >  8:245473/1‑149 (MQ=255)                                                                                                                                    tgcttGGCGGCCAGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTGCGCGCTTGCAGGTCTGCCGGCAGTTCGGCGAGGCGCTGCAg          >  5:134740/5‑149 (MQ=255)                                                                                                                                      cttGGCGGCCAGGCCGCC‑AGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCgcggcgg                                                                    >  2:155642/3‑88 (MQ=255)                                                                                                                                      cttGGCGGCCAGGCCGCC‑AGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCgcggcgg                                                                    <  1:155642/86‑1 (MQ=255)                                                                                                                                        tGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAg                                                      >  6:36250/1‑101 (MQ=255)                                                                                                                                        tGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAg                                                      <  5:36250/101‑1 (MQ=255)                                                                                                                                            ggcCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTGCGCGCTTGCAGGTCTGCCGGCAGTTCGGCGAGGCGCTGCAGGCTGGGGt  >  7:121253/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                  CTCCAGGGTGCCGGCGCTTTTCACGGCGAAGGCCAGGCTGGTCCAGCCAGCCAGGGTCTGATCGTTGGTTGGCTGCTGCAGTTGCTCTGCAGGCAGCGAGGCGACCAAGGCCCAGATGGCGTCATT‑GTTGGTGCTGGCGGCCTGGCCGCCAAGCAGCGGGGCCAGCAGGATGCGCTGCTGTGCGGCAGCCAGGGCCTGGCCATCG‑GCTTCGAGCGCGGCGGCGTGCACACTGTAGGTGCGCGCTTGCAGGTCTGCCGGCAGTTCGGCGAGGCGCTGCAGGCTGGGGT  >  NC_002947/1510928‑1511214

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: ‑