Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,675,334	A→G	13.3%	intergenic (+33/+71)	PP_4137 → / ← chrR	outer membrane siderophore receptor/chromate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,675,334	0	A	G	13.3%	69.6 / 3.8	30	intergenic (+33/+71)	PP_4137/chrR	outer membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand):  ref base A (9/17);  new base G (2/2);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCTTGTCA  >  NC_002947/4675193‑4675480                                                                                                                                              |                                                                                                                                                    tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGaa                                                                                                                                                <  8:164562/147‑1 (MQ=255)  gCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAttt                                                                                                                                             >  8:287036/1‑149 (MQ=255)             gACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCAtt                                                                                                                                  >  3:17157/1‑149 (MQ=255)              aCGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCAtt                                                                                                                                  <  7:287036/148‑1 (MQ=255)                gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTcgc                                                                                                                               <  4:17157/149‑1 (MQ=255)                 ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTcgcg                                                                                                                              <  7:22750/149‑1 (MQ=37)                  cGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCgg                                                                                                                             <  8:325827/149‑1 (MQ=255)                  cGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCgg                                                                                                                             <  5:135528/149‑1 (MQ=255)                            cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGt                                                                                                                           <  8:66503/141‑1 (MQ=37)                            cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGt                                                                                                                           >  7:66503/1‑141 (MQ=37)                              acgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTc                                                                                                                 <  8:261885/149‑1 (MQ=255)                                             aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTg                                                                                                                                                  <  2:272346/101‑1 (MQ=255)                                             aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTg                                                                                                                                                  >  1:272346/1‑101 (MQ=255)                                                                 gAAGGATTTCGCGTGCATGCCGGGGTTTCCTACAGCGACCGCTGACGGGGGGGCGCGGGTTTACCCGCGCATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCGCCTACACAGGGGCGTTGCACAGCAGGATTAgggggt                                                                              <  8:174875/149‑3 (MQ=38)                                                                           gCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGTGTTTACCCGTGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                     <  1:109932/148‑1 (MQ=255)                                                                                tATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGa                                                                          <  5:319707/138‑1 (MQ=255)                                                                                tATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                               <  3:157935/149‑1 (MQ=255)                                                                                   gCCGGGGTTTCCTACCGCTTCCGCGGACGGGGGGGCGCGGGTTTACCCGCTAATGCGTCATTGAATTTGCCGCCGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                             <  3:60119/148‑1 (MQ=38)                                                                                                       ccGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                          <  7:61290/147‑1 (MQ=255)                                                                                                        cGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGAATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAACAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                        <  5:18954/148‑1 (MQ=255)                                                                                                           tGACGGGGGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCTACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCg                                    <  7:214370/149‑1 (MQ=255)                                                                                                                           gTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAg                    >  8:301108/1‑148 (MQ=255)                                                                                                                               aCCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                                      >  4:99085/1‑127 (MQ=255)                                                                                                                               aCCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                                      <  3:99085/127‑1 (MQ=255)                                                                                                                                 ccGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTc             >  6:31295/1‑149 (MQ=37)                                                                                                                                 ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTc             >  6:280945/1‑149 (MQ=255)                                                                                                                                 ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTc             >  3:316372/1‑149 (MQ=255)                                                                                                                                    cgAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTcgcg          >  2:209617/1‑149 (MQ=255)                                                                                                                                          cGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGCACAAACGCCTGCAAGAACGGTCGCGTCtt      <  1:209617/148‑1 (MQ=255)                                                                                                                                            tCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa  >  6:330613/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                    TGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCTTGTCA  >  NC_002947/4675193‑4675480

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