Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,076,092	A→G	10.8%	intergenic (‑40/+34)	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,092	0	A	G	10.8%	113.1 / 5.3	46	intergenic (‑40/+34)	PP_2685/PP_2686	hypothetical protein/transglutaminase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (3/2);  total (24/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 = 8.85e-01

GACGAACAACTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGT  >  NC_002947/3075945‑3076238                                                                                                                                                    |                                                                                                                                                   gACGAACAACTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCAc                                                                                                                                                   >  6:228820/1‑149 (MQ=255)     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)                                                                                                                                                    |                                                                                                                                                   GACGAACAACTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGT  >  NC_002947/3075945‑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: ‑