Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,745,547	C→A	9.3%	E86* (GAA→TAA) ‡	PP_1543 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,745,547	0	C	A	9.3%	126.3 / 3.0	43	E86* (GAA→TAA) ‡	PP_1543	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/21);  new base A (2/2);  total (20/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 = 7.92e-01

CGGTGGGCGCCGCCGCACTCATCGAGCGAGCTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGCCCCAGTAGGCGTGCCAGGTTTTGTCGTAGCAG  >  NC_002947/1745402‑1745686                                                                                                                                                  |                                                                                                                                            cgGTGGGCGCCGCCGCACTCATCGAGCGAGCTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAgcg                                                                                                                                         <  2:542162/150‑1 (MQ=255)         gccgccgcACTCATCGAGCGAGCTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTGAAGCGCTggccga                                                                                                                                 <  2:295643/150‑6 (MQ=255)                   tCATCGAGCGAGCTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGAGCCCGGTCACGCCGCATTTGCACGATCACGTTCCGCGCTTTGTCGGCCAGCGCTGCAGCGCTGAACCGGCGAGAACTc                                                                                                                       >  2:228238/1‑150 (MQ=255)                        gagcgagcTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCt                                                                                                                                       <  2:481664/129‑1 (MQ=255)                        gagcgagcTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCt                                                                                                                                       >  1:481664/1‑129 (MQ=255)                         agcgagcTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGAc                                                                                                                 <  2:217994/150‑1 (MQ=255)                           cgagcTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCg                                                                                                               <  2:416078/150‑1 (MQ=255)                                  gtgtgtCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGaa                                                                                                         >  2:753682/1‑149 (MQ=255)                                  gtgtgtCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGaa                                                                                                         <  1:753682/149‑1 (MQ=255)                                         gAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCga                                                                                                 >  1:792086/1‑150 (MQ=255)                                                       tcgtcgAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAgcgc                                                                                                                                        <  2:74576/97‑1 (MQ=255)                                                       tcgtcgAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAgcgc                                                                                                                                        >  1:74576/1‑97 (MQ=255)                                                          tcgAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGc                                                                                <  2:40213/150‑1 (MQ=255)                                                              gCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCt                                                                            <  2:435844/150‑1 (MQ=255)                                                                atgCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTg                                                                          >  1:381731/3‑150 (MQ=255)                                                                       gggCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCagaaga                                                                   >  2:186528/1‑150 (MQ=255)                                                                        ggCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAgaagaa                                                                  >  1:486599/1‑150 (MQ=255)                                                                                  gCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATg                                                        <  1:57230/150‑1 (MQ=255)                                                                                   cGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATgg                                                       >  2:345521/1‑150 (MQ=255)                                                                                    gTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTg                                                               >  2:763453/1‑141 (MQ=255)                                                                                    gTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGt                                                      >  2:213575/1‑150 (MQ=255)                                                                                    gTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGt                                                      >  2:64557/1‑150 (MQ=255)                                                                                    gTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGGAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTg                                                               <  1:763453/141‑1 (MQ=255)                                                                                       cAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAg                                                   >  1:713281/1‑150 (MQ=255)                                                                                              cGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATcc                                            <  1:64557/150‑1 (MQ=255)                                                                                                gTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCca                                          <  1:186186/150‑1 (MQ=255)                                                                                                  cACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCcaca                                        >  1:531203/1‑150 (MQ=255)                                                                                                   aCGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACAt                                       >  1:788997/1‑150 (MQ=255)                                                                                                        gcATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATg                                                        >  2:165252/1‑128 (MQ=255)                                                                                                        gcATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATg                                                        <  1:165252/128‑1 (MQ=255)                                                                                                         cATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGcc                                 <  2:454413/150‑1 (MQ=255)                                                                                                          aTTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGccc                                <  1:345521/150‑1 (MQ=255)                                                                                                              gCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGCCCCAGt                            <  1:213575/150‑1 (MQ=255)                                                                                                                     cACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATc                                             >  2:332863/1‑126 (MQ=255)                                                                                                                     cACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATc                                             <  1:332863/126‑1 (MQ=255)                                                                                                                             gcgcTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGCCCCAGt                            >  1:582215/1‑135 (MQ=255)                                                                                                                             gcgcTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGCCCCAGt                            <  2:582215/135‑1 (MQ=255)                                                                                                                                cTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGCCCCAGTAGGCGTGCCAGGTTTTGt          <  2:792086/150‑1 (MQ=255)                                                                                                                                    gTCGGCCAGCGCTCCCCCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATACCACATACCGCCCCAGTAGGCGTGCCAGGTTTTGTCGTa      <  1:228238/150‑1 (MQ=255)                                                                                                                                       gttcAGCGCTGAAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGCCCCAGTAGGCGTGCCAGGTTTTGTCGTagca   >  2:490897/4‑150 (MQ=255)                                                                                                                                        ttcAGCGCTGAAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAg                                                   >  2:176169/3‑101 (MQ=255)                                                                                                                                        ttcAGCGCTGAAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAg                                                   <  1:176169/99‑1 (MQ=255)                                                                                                                                        gCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGCCCCAGTAGGCGTGCCAGGTTTTGTCGTagcag  <  2:713281/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                            CGGTGGGCGCCGCCGCACTCATCGAGCGAGCTGGTGTGTCGAACATCCTCAGCCTCGTCGAGCAGGCTTCGGGCATCTTCTGCGTCCAGATCCCGGTCACGCCGCATTTGCACGATCACCTTCCGCGCTTTGTCGGCCAGCGCTTCAGCGCTGAACCGGCGAGAACTCAGCGACCGGTCGAAGTAGCCGATGATGTAGGCGTCGTGCAGCTTGCAGAAGAACTGGCCGATGGTCAGGCCATCCCACATACCGCCCCAGTAGGCGTGCCAGGTTTTGTCGTAGCAG  >  NC_002947/1745402‑1745686

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