Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,076,095	C→A	11.4%	intergenic (‑43/+31)	PP_2685 ← / ← PP_2686	hypothetical protein/transglutaminase domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,076,095	0	C	A	11.4%	116.4 / 7.2	44	intergenic (‑43/+31)	PP_2685/PP_2686	hypothetical protein/transglutaminase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base A (3/2);  total (22/22)
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.86e-01

AACAACTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGT  >  NC_002947/3075949‑3076238                                                                                                                                                   |                                                                                                                                                aacaacTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTcgc                                                                                                                                               >  8:281131/1‑149 (MQ=255)               cGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGt                                                                                                                                 <  6:321869/149‑1 (MQ=255)               cGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTACGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGt                                                                                                                                 <  8:26643/149‑1 (MQ=255)                gTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTc                                                                                                                                <  2:184066/149‑1 (MQ=255)                    agaTCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGaa                                                                                                                            <  6:184628/149‑1 (MQ=255)                     gaTCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCCGTCAGAAg                                                                                                                           >  6:111803/1‑149 (MQ=255)                         tGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAg                                                                                                                       >  7:7367/1‑149 (MQ=255)                            tCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCt                                                                                                                    >  6:311767/1‑149 (MQ=255)                              gATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGt                                                                                                                  <  8:7367/149‑1 (MQ=255)                                    ggCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTgg                                                                                                            <  7:281131/149‑1 (MQ=255)                                       aTTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGa                                                                                                         >  6:289693/1‑149 (MQ=255)                                       aTTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCAGGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTTGAAGGGTCTTGCCGCAGTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGa                                                                                                         <  5:228820/149‑1 (MQ=255)                                           gTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGGCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGc                                                                                                      <  5:111803/148‑1 (MQ=255)                                             gctttaGTCGCTGATTAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTATGTCATGCGGTTTCCCTTGTTTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGAAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCAc                                                                                                    <  4:307586/143‑1 (MQ=255)                                                   gTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACg                                                                                             <  8:5911/149‑1 (MQ=255)                                                              aGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCAt                                                                                  <  2:181773/149‑1 (MQ=255)                                                                aCAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCGCGAGTGCGGCAAGa                                                                                                                                      >  3:92594/1‑95 (MQ=39)                                                                aCAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCGCGAGTGCGGCAAGa                                                                                                                                      <  4:92594/95‑1 (MQ=39)                                                                  aaGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGAt                                                                              >  3:183350/1‑149 (MQ=255)                                                                     cccGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCATTg                                                                                                             >  7:34367/1‑115 (MQ=255)                                                                     cccGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCATTg                                                                                                             <  8:34367/115‑1 (MQ=255)                                                                       cggcctccAGGGGCAGGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGGGGAAGGGTCTGGCCGCACGCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGgc                                                                         <  1:153422/146‑1 (MQ=255)                                                                       cGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGgc                                                                         <  5:289693/149‑1 (MQ=255)                                                                       cGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGGCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGgc                                                                         <  5:311767/149‑1 (MQ=255)                                                                           ctccAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGACGCACTCGCGCGGAAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGAACGTGAACGCTGGCATGCATCGATTCGGcgccg                                                                     >  4:279218/1‑149 (MQ=255)                                                                                     aTGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCgcggc                                                           <  4:183350/149‑1 (MQ=255)                                                                                               aGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCt                                                 >  5:230258/1‑149 (MQ=255)                                                                                                         gCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGggcag                                       >  3:306165/1‑149 (MQ=255)                                                                                                          cGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGggcagg                                      <  4:306165/149‑1 (MQ=255)                                                                                                                   ttGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTAACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAgg                              >  3:314919/1‑148 (MQ=255)                                                                                                                          gCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAg                      >  6:303510/1‑149 (MQ=255)                                                                                                                           cGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAg                      <  6:230258/148‑1 (MQ=255)                                                                                                                            gtgtGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGAt                                                                              <  6:196521/91‑1 (MQ=255)                                                                                                                            gtgtGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGAt                                                                              >  5:196521/1‑91 (MQ=255)                                                                                                                                    gggTCTTGCCGCACTCGCGCGGCAAGCCCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGcc            >  7:5398/1‑149 (MQ=255)                                                                                                                                    gggTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACAGGGCAGGCATCAAGGTAGTCCAGGCCCACCGc             >  7:739/1‑148 (MQ=255)                                                                                                                                     ggTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGcc                   >  5:326925/1‑141 (MQ=255)                                                                                                                                     ggTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGcc                   <  6:326925/141‑1 (MQ=255)                                                                                                                                      gTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGTCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAg          >  6:135171/1‑149 (MQ=255)                                                                                                                                       tCTTGCCGCGCGAGTGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAg                               >  2:225453/1‑127 (MQ=255)                                                                                                                                       tCTTGCCGCGCGAGTGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAg                               <  1:225453/127‑1 (MQ=255)                                                                                                                                         ttGCCGCGCGAGTGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCtt       >  1:157048/1‑149 (MQ=255)                                                                                                                                         ttGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCtt       >  1:192230/1‑149 (MQ=255)                                                                                                                                          tGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCtt       >  2:12032/1‑148 (MQ=255)                                                                                                                                               cACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGCCCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCGTCAGGt  >  5:230873/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                AACAACTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGT  >  NC_002947/3075949‑3076238

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