Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	337,809	A→G	10.4%	intergenic (‑59/+70)	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,809	0	A	G	10.4%	119.5 / 6.0	48	intergenic (‑59/+70)	PP_0279/PP_0280	hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base G (3/2);  total (25/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 = 2.08e-01

GTTTCCACAGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAC  >  NC_002947/337663‑337948                                                                                                                                                   |                                                                                                                                               gTTTCCACAGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGtata                                                                                                                                              <  5:56404/149‑1 (MQ=255) gTTTCCACAGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGta                                                                                                                                                >  7:182457/1‑147 (MQ=255)      cacaGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGcc                                                                                                                                        <  1:161164/149‑1 (MQ=255)       acaGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCc‑a                                                                                                                                      >  6:343728/1‑149 (MQ=255)         aGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCc‑a‑c                                                                                                                                    >  8:326748/1‑148 (MQ=255)              ctgGTGCTCGATCAGGTACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCaa                                                                                                                              <  4:260567/149‑1 (MQ=255)              ccggTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTc                                                                                                                                >  2:305456/3‑147 (MQ=255)                ggTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAgg                                                                                                                            >  6:90406/1‑149 (MQ=255)                       gATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGgcgc                                                                                                                     <  7:326748/149‑1 (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‑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               <  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)                                                                                                                                                   |                                                                                                                                               GTTTCCACAGCTGCTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATA‑GCGCC‑A‑CATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAC  >  NC_002947/337663‑337948

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