Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	337,827	Δ1 bp	13.0%	intergenic (‑77/+52)	PP_0279 ← / ← PP_0280	hypothetical protein/amino acid ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	337,827	0	T	.	13.0%	127.8 / 15.7	46	intergenic (‑77/+52)	PP_0279/PP_0280	hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base . (3/3);  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 = 1.39e-01

GATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGG  >  NC_002947/337685‑337960                                                                                                                                                  |                                                                                                                                      gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcgc                                                                                                                                 <  1:184546/149‑1 (MQ=255) gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcgc                                                                                                                                 <  7:326748/149‑1 (MQ=255) gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcgc                                                                                                                                 <  5:90406/149‑1 (MQ=255) gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcg                                                                                                                                  <  5:281327/148‑1 (MQ=255) gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgc                                                                                                                                   <  2:309638/147‑1 (MQ=255) gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgc                                                                                                                                   <  1:122428/147‑1 (MQ=255)     aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTaca                                                                                                                             >  1:306476/1‑149 (MQ=255)      ggAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTacac                                                                                                                            >  2:217841/1‑149 (MQ=255)            ggTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                      >  3:224098/1‑149 (MQ=255)             gTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTA‑GCGCCAA‑CCTTGAATG‑TGGCGCTATACCTGTgg                                                                                                                      <  3:154881/148‑1 (MQ=25)             gTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTgg                                                                                                                      >  3:243436/1‑148 (MQ=255)              tCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGGCGCA‑A‑CCGTGAATGTTGG‑GCGATACCTGTgggg                                                                                                                    >  7:348460/1‑148 (MQ=37)                 aGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCg                                                                                                                 <  4:224098/149‑1 (MQ=255)                  gggCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCgg                                                                                                                >  5:147382/1‑149 (MQ=255)                         gAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGcc                                                                                                         >  5:357041/1‑149 (MQ=255)                             ggggTGCTTCGGGGGGGGCGGGGGTGGGGGTTGGGCTTGGCATGGATGGGGGGTGGTTTTGCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTTATA‑CCGCA‑ACCCTTTAATG‑TGGCGCTATACCTGGGGGAGGGGGCGTGCCcgcg                                                                                                     <  8:348460/148‑1 (MQ=9)                               ggTGCTTTGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                   <  6:147382/149‑1 (MQ=255)                                   cTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCTCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcag                                                                                               >  7:235461/1‑149 (MQ=255)                                    ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCt                                                                                                                          <  3:46008/121‑1 (MQ=255)                                    ttAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCt                                                                                                                          >  4:46008/1‑121 (MQ=255)                                          gggggCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGa                                                                                        <  5:343728/149‑1 (MQ=255)                                                             gggCATGGCAAGGAGGTGTGGTGGTTTTGTGGGCCTCTTCGCGGGCACGCCGGCTCCCACAGGTATA‑GCGCC‑A‑CATTCACGGTGGGCGCTACACGTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggca                                                                     <  3:233103/149‑1 (MQ=38)                                                              ggCAGGGCACGTATGTGTGGTGGTTTTGAGGGCGTATTCCCGGGCACGCCCGCTCCCCCTGTTATC‑GCGCC‑A‑CATTCATGGTTGGCGCTACACGTGTGGGAGGGGGGGGGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggcag                                                                    <  4:243436/149‑1 (MQ=255)                                                                   ggCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATc                                                               >  5:141382/1‑149 (MQ=255)                                                                                            ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCagcagc                                      <  1:217841/149‑1 (MQ=255)                                                                                                tctTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑C‑CATTCAAGGTTGGCGCTACCCCTGTGGGAGCGGGCGTGCCCGCGAGGCAGTCAGTGCGGGCTCAGCGGCCTTGGCAGGCAGCCACGCGCAACCAGCGCTCCAGCAGCTTg                                   >  5:65638/1‑148 (MQ=38)                                                                                                  ttCGCGGGCACGCCCGCTCCCACAGGTGTA‑GCGCCAA‑CCTTGAATG‑TGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAg                                >  7:75095/1‑149 (MQ=21)                                                                                                      cgggaccGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCACGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCtt                            <  8:235461/143‑1 (MQ=255)                                                                                                       gggCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                           <  1:305456/149‑1 (MQ=255)                                                                                                       gggCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                           <  6:357041/149‑1 (MQ=255)                                                                                                         gCACGCCCGCTCCCACACGTATA‑GCGCC‑A‑CATTCAAGGGTGGGGCTACACCAGTGGGAGCGGGCCTGCCCGCGAAGGAGTGATTGAGGTCTAAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCa                         >  2:340710/1‑149 (MQ=255)                                                                                                             gCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCCCCAg                     >  6:365418/1‑149 (MQ=255)                                                                                                               ccGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCa                   >  8:83459/1‑149 (MQ=255)                                                                                                               ccGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCa                   >  2:119944/1‑149 (MQ=255)                                                                                                                     ccACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCagcagc                                      <  7:288876/124‑1 (MQ=255)                                                                                                                     ccACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCagcagc                                      >  8:288876/1‑124 (MQ=255)                                                                                                                     ccACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc              <  5:365418/148‑1 (MQ=255)                                                                                                                      cacaGGTGTA‑GCGCCAA‑CCTTGAATG‑TGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggc                                                                      >  6:317484/1‑91 (MQ=14)                                                                                                                      cacaGGTGTA‑GCGCCAA‑CCTTGAATG‑TGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTggc                                                                      <  5:317484/91‑1 (MQ=14)                                                                                                                        caGGTGTA‑GCGCCAA‑CCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc              >  8:15897/1‑145 (MQ=25)                                                                                                                         aGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc              >  3:26221/1‑144 (MQ=255)                                                                                                                         aGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc              >  2:223809/1‑144 (MQ=255)                                                                                                                                  cgcC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg  >  8:157368/1‑148 (MQ=255)                                                                                                                                    cc‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg  >  4:102464/1‑146 (MQ=255)                                                                                                                                     cga‑caTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCACGCGAAGCAGTCAGAGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTAGCACCAGCATGACCGAGATCACCAgg  >  5:257167/3‑146 (MQ=255)                                                                                                                                       a‑caTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCt                             <  2:297492/117‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                      GATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGG  >  NC_002947/337685‑337960

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