Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,675,343	G→A	13.3%	intergenic (+42/+62)	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,343	0	G	A	13.3%	70.7 / 3.2	30	intergenic (+42/+62)	PP_4137/chrR	outer membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand):  ref base G (10/16);  new base A (2/2);  total (12/18)
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.48e-01

GACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCT  >  NC_002947/4675205‑4675487                                                                                                                                           |                                                                                                                                                  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)                                                     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)                                                                                                                                   ttGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCtt             <  3:81774/142‑1 (MQ=255)                                                                                                                                   ttGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCtt             >  4:81774/1‑142 (MQ=255)                                                                                                                                   ttGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt       >  6:120436/1‑148 (MQ=255)                                                                                                                                       aTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt  >  4:175797/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                  GACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCG‑ACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCT  >  NC_002947/4675205‑4675487

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