Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,500,570	G→A	12.5%	intergenic (+17/+8)	PP_2193 → / ← PP_2194	ferric siderophore receptor/LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,500,570	0	G	A	12.5%	78.2 / 3.3	32	intergenic (+17/+8)	PP_2193/PP_2194	ferric siderophore receptor/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (16/12);  new base A (2/2);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

AAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGCCCATAACAGCCCC  >  NC_002947/2500425‑2500707                                                                                                                                                  |                                                                                                                                          aaGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGc                                                                                                                                        >  8:91507/1‑149 (MQ=255) aaGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACAGGc                                                                                                                                        >  3:274481/1‑149 (MQ=255)      gTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGtt                                                                                                                                   <  7:243596/149‑1 (MQ=255)         cAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCa                                                                                                                                 <  7:152239/148‑1 (MQ=255)           gCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACAGGCCCGTTCAGGc                                                                                                                              <  4:274481/149‑1 (MQ=255)                gTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGTGCCGGAGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAAcg                                                                                                                         >  6:292255/1‑149 (MQ=255)                  aTGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACAGGCCCGTTCAGGccaccac                                                                                                                        >  6:144201/1‑142 (MQ=255)                  aTGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACAGGCCCGTTCAGGccaccac                                                                                                                        <  5:144201/148‑7 (MQ=255)                    ggggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGccac                                                                                                                           >  3:267924/1‑140 (MQ=255)                    ggggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGccac                                                                                                                           <  4:267924/143‑4 (MQ=255)                      ggACAGCGTGCTACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCa                                                                                                                   <  4:166656/149‑1 (MQ=255)                                             ttCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATgg                                                                                            >  1:43630/1‑149 (MQ=255)                                                gCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCt                                                                                         >  7:23770/1‑149 (MQ=255)                                                       tACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATgg                                                                                            >  5:189546/1‑139 (MQ=255)                                                       tACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATgg                                                                                            <  6:189546/139‑1 (MQ=255)                                                                      cTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGccc                                                                   >  3:211851/1‑149 (MQ=255)                                                                       tGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCg                                                                  >  4:199925/1‑149 (MQ=255)                                                                           tGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTg                                                               <  3:199925/148‑1 (MQ=255)                                                                              cTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCgg                                                           <  4:211851/149‑1 (MQ=255)                                                                               tACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGccctt                                                          >  1:110111/1‑146 (MQ=255)                                                                                    ggggCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTG‑GCCCGGCACCGGGcacagg                                                    <  8:23770/149‑6 (MQ=255)                                                                                         gTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCCCCAACCCCCTGTCCCCGGAGCCGGGCCCggcagag                                                >  6:237063/1‑144 (MQ=255)                                                                                         gTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGGGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAgg                                                     <  5:237063/144‑1 (MQ=255)                                                                                            ggCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACgg                                             <  2:43630/149‑1 (MQ=255)                                                                                                              gCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCATCGCCAAACCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCt                            >  7:43193/1‑148 (MQ=255)                                                                                                                           ttACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAAc                                               <  7:12555/116‑1 (MQ=255)                                                                                                                           ttACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAAc                                               >  8:12555/1‑116 (MQ=255)                                                                                                                           ttACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGcc              >  7:51605/1‑149 (MQ=255)                                                                                                                                     tggtggCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGCCCATAACAGcc    >  8:114512/1‑149 (MQ=255)                                                                                                                                       gtggCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGCCCATAACAGcccc  >  1:203454/1‑149 (MQ=255)                                                                                                                                               aaCGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCATATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAgcgc               <  4:78132/128‑1 (MQ=255)                                                                                                                                               aaCGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAGATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAgcgc               >  3:78132/1‑128 (MQ=255)                                                                                                                                                  |                                                                                                                                          AAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGCCCATAACAGCCCC  >  NC_002947/2500425‑2500707

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