Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,707,246	A→G	8.3%	intergenic (+64/+78)	PP_1501 → / ← PP_1502	hypothetical protein/OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,707,246	0	A	G	8.3%	121.9 / 2.4	48	intergenic (+64/+78)	PP_1501/PP_1502	hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base A (24/20);  new base G (2/2);  total (26/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 = 6.85e-01

CTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCAC  >  NC_002947/1707105‑1707392                                                                                                                                               |                                                                                                                                                    cTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcgg                                                                                                                                              >  7:253495/1‑149 (MQ=255) cTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcgg                                                                                                                                              >  7:230206/1‑149 (MQ=255)  tGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcgg                                                                                                                                              >  8:300683/1‑148 (MQ=255)   gACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGca                                                                                                                                            >  8:30255/1‑149 (MQ=255)    aCCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGca                                                                                                                                            >  1:207738/1‑148 (MQ=255)          gTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGca                                                                                                                                            <  1:353311/142‑1 (MQ=255)                gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTCAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTAGATCTAACAGGCGGCACAACATTATCGGc                                                                                                                               >  8:220630/1‑149 (MQ=255)                    aCAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGaa                                                                                                                            >  8:284385/1‑148 (MQ=255)                        ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                                               >  8:312609/1‑141 (MQ=255)                        ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                                               <  7:312609/141‑1 (MQ=255)                        ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGaa                                                                                                                            >  6:191538/1‑144 (MQ=255)                        ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGaa                                                                                                                            <  5:191538/144‑1 (MQ=255)                          gCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGa                                                                                                                      >  6:222320/1‑148 (MQ=255)                            aaGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGGGGCACAACATTATCGGCGAACAAAGaa                                                                                                                     >  5:323482/1‑147 (MQ=255)                              gcgcgcGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGc                                                                                                                 <  7:30255/149‑1 (MQ=255)                                     gTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑t                                                                                                         >  1:241860/1‑149 (MQ=255)                                       gcgcTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑Tc                                                                                                        >  3:174151/1‑148 (MQ=255)                                              aTGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCgg                                                                                                <  5:291931/149‑1 (MQ=255)                                               tGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCggg                                                                                               >  7:48092/1‑149 (MQ=255)                                                 ggCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCgggg                                                                                              <  5:222320/148‑1 (MQ=255)                                                 ggCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGc                                                                                             <  8:230206/149‑1 (MQ=255)                                                  gCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GGAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGACAGGCGGCACAACATTATGGGGGAACAAAGGAAAGCCCCGCG‑TCGATCGCgggggt                                                                                            >  4:276431/1‑147 (MQ=39)                                                    gAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTg                                                                                          >  7:315816/1‑149 (MQ=255)                                                      ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTc                                                                                        <  7:110759/149‑1 (MQ=255)                                                        aaaCACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCt                                                                                      <  2:241860/149‑1 (MQ=255)                                                          acacCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCttt                                                                                    >  7:335882/1‑149 (MQ=255)                                                           cacCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTc                                                                                   <  3:19145/149‑1 (MQ=255)                                                           cacCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTc                                                                                   <  4:174151/149‑1 (MQ=255)                                                            acCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACCCGCGGCCCAACATTATCGGCGGACAAAGAAAAGCCCCGCG‑TCGCACGCGGGGGGTGGGCTTTCa                                                                                  >  2:341899/1‑149 (MQ=38)                                                            acCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATcgcg                                                                                                 <  3:31506/134‑1 (MQ=255)                                                            acCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATcgcg                                                                                                 >  4:31506/1‑134 (MQ=255)                                                            acCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCttt                                                                                    >  8:98475/1‑147 (MQ=255)                                                             ccGCGGCAGGTGCAGTGAGGTTTCGGCGGCAGCATTTTTGCCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCttt                                                                                    >  1:21504/1‑147 (MQ=255)                                                                          aGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTa                                                                    <  1:341899/149‑1 (MQ=255)                                                                            tGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgc                                                                  <  7:284385/149‑1 (MQ=255)                                                                                 tttCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTg                                                              >  1:315736/1‑148 (MQ=255)                                                                                         gcagcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGc                                                      <  8:253495/149‑1 (MQ=255)                                                                                               tttttGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGGACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑GCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCATAGCGCTGAACCTGGCTCACCg                                                >  1:158850/1‑149 (MQ=255)                                                                                                   tGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGc                                             <  8:315816/148‑1 (MQ=255)                                                                                                       gcctgtgagaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                        >  6:225805/7‑149 (MQ=38)                                                                                                       gcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCTGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                        >  4:71174/1‑149 (MQ=255)                                                                                                              agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCg                                 <  5:329277/149‑1 (MQ=255)                                                                                                               gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTagccag                                                                         <  2:327476/108‑1 (MQ=255)                                                                                                               gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTagccag                                                                         >  1:327476/1‑108 (MQ=255)                                                                                                                     gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt                          <  7:300683/149‑1 (MQ=255)                                                                                                                          gcccgcgcggcgCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGt                     <  5:225805/143‑1 (MQ=255)                                                                                                                           cgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt                    <  7:98475/149‑1 (MQ=255)                                                                                                                            gcgggcggCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAATCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTg                   <  8:48092/149‑1 (MQ=255)                                                                                                                                               acaGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCAc  <  2:207738/147‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    CTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGC‑GCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCAC  >  NC_002947/1707105‑1707392

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑