Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,033,553	C→G	9.8%	R32G (CGC→GGC)	PP_2646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,033,553	0	C	G	9.8%	116.2 / 2.7	39	R32G (CGC→GGC)	PP_2646	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/17);  new base G (2/2);  total (20/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

CCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACG  >  NC_002947/3033407‑3033692                                                                                                                                                   |                                                                                                                                                    cccAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                                   >  1:74273/1‑149 (MQ=255) cccAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                                   >  8:158782/1‑149 (MQ=255)           gAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                                   <  1:108541/139‑1 (MQ=255)           gAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                                   >  2:108541/1‑139 (MQ=255)             aGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑c                                                                                                                                <  2:185751/149‑1 (MQ=255)              gACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑cc                                                                                                                               >  8:20937/1‑149 (MQ=255)                cccTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCgg                                                                                                                             <  3:243892/149‑1 (MQ=255)                  ctctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGg‑t                                                                                                                          >  1:141320/1‑149 (MQ=255)                              cgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACgg                                                                                                              <  2:141320/149‑1 (MQ=255)                                    ggCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAg                                                                                                        >  7:284690/1‑149 (MQ=255)                                     gCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGc                                                                                                       <  7:158782/149‑1 (MQ=255)                                       cAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCa                                                                                                     <  2:18799/149‑1 (MQ=255)                                                          aaGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                                  <  2:74273/149‑1 (MQ=255)                                                            ggAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGa                                                                                <  2:84252/149‑1 (MQ=255)                                                             gagaGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGa                                                                                >  4:85944/1‑148 (MQ=255)                                                             gagaGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCATGGTTGTGAAATCGa                                                                                >  6:185472/1‑148 (MQ=255)                                                             gagaGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCAGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAA‑CCCGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGaaaaaga                                                                                >  2:302298/1‑144 (MQ=255)                                                                  cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCCCGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTg                                                                          >  7:106042/1‑149 (MQ=255)                                                                         tGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGgcgcgc                                                                   <  1:302298/149‑1 (MQ=255)                                                                                aCTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGTCGCGCGAGCATg                                                            >  4:201179/1‑149 (MQ=255)                                                                                     tatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGaa                                                                               <  4:242398/125‑1 (MQ=255)                                                                                     tatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGaa                                                                               >  3:242398/1‑125 (MQ=255)                                                                                      atGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAg                                                                >  1:123951/1‑139 (MQ=255)                                                                                      atGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAg                                                                <  2:123951/139‑1 (MQ=255)                                                                                            aGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCAGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑ACGGG‑TGTTGATCCACGGTGGCAGCCAGTTTATTGGTTGGGAAATAGAAAGCTGGGCGCGCGAGCATGGAGCCGATCttg                                                >  6:205956/1‑147 (MQ=255)                                                                                             tggCTCAACTACTCCTACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCATCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGt                                               <  8:284690/148‑1 (MQ=255)                                                                                                 gcAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGt                                           >  7:150951/1‑149 (MQ=255)                                                                                                   aaCTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTa                                         <  5:190482/149‑1 (MQ=255)                                                                                                    aCTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCGAGCGGC‑AGCGGCTCACCCTTTCCCGGGCATTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCTCGAGCATGGCGCCGATCTGGt                                               <  1:100222/149‑1 (MQ=255)                                                                                                         gCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATAGAAAGCTGGGCGCGCGAGCATAGCGCCGATCTGGTGCGTTACCCGc                                    >  8:118134/1‑148 (MQ=255)                                                                                                             acacCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGGGTTACCGGCCAAAc                               >  7:300601/1‑149 (MQ=255)                                                                                                             acacCCGGATCGGCCG‑CTGGCGCTGGACGTCGTCCAGGGCCTGGCGG‑‑‑‑‑CAGCG‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGt                                               <  4:255747/133‑1 (MQ=17)                                                                                                             acacCCGGATCGGCCG‑CTGGCGCTGGACGTCGTCCAGGGCCTGGCGG‑‑‑‑‑CAGCG‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGt                                               >  3:255747/1‑133 (MQ=17)                                                                                                                 ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTAcc                                       >  8:302020/1‑137 (MQ=255)                                                                                                                 ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTAcc                                       <  7:302020/137‑1 (MQ=255)                                                                                                                        tGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGGCAGTTTCTTGGTGGGGACATCGAAAGCAGGGCGCGCGAGGATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGCCAGCTGCa                    >  5:309374/1‑149 (MQ=255)                                                                                                                         gccgccAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCa                          >  3:185456/1‑142 (MQ=255)                                                                                                                         gccgccAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCa                          <  4:185456/142‑1 (MQ=255)                                                                                                                                   ccTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAgg         <  8:150951/149‑1 (MQ=255)                                                                                                                                    cTGGACGACGTCTATGGCC‑AGCTG‑‑‑‑‑CCGAT‑CGGTG‑TGGTGATCCACTGTGTCAGCCAGTTTCTTTGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGc        <  5:185472/149‑1 (MQ=255)                                                                                                                                    cTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCCTGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGc        >  5:37276/1‑149 (MQ=255)                                                                                                                                         cgacgTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAAc   <  7:20937/149‑1 (MQ=255)                                                                                                                                          gtcgTCCAGGGCC‑TGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTAGATCACCGTTGTCACCCAATAAATAAAAAAAAAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACg  >  2:246103/3‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGG‑‑‑‑‑CCGAT‑CCGGG‑TGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACG  >  NC_002947/3033407‑3033692

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