Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	337,825	G→T	14.6%	intergenic (‑75/+54)	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,825	0	G	T	14.6%	123.8 / 8.9	46	intergenic (‑75/+54)	PP_0279/PP_0280	hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base T (4/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.00e+00

GGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGG  >  NC_002947/337678‑337960                                                                                                                                                       |                                                                                                                                        ggTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAgg                                                                                                                                        >  6:90406/1‑149 (MQ=255)        gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcgc                                                                                                                                 <  1:184546/149‑1 (MQ=255)        gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcgc                                                                                                                                 <  5:90406/149‑1 (MQ=255)        gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcgc                                                                                                                                 <  7:326748/149‑1 (MQ=255)        gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcg                                                                                                                                  <  5:281327/148‑1 (MQ=255)        gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgc                                                                                                                                   <  1:122428/147‑1 (MQ=255)        gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgc                                                                                                                                   <  2:309638/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                           <  6:357041/149‑1 (MQ=255)                                                                                                              gggCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTg                           <  1:305456/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)                                                                                                                                                       |                                                                                                                                        GGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGG  >  NC_002947/337678‑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: ‑