Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,125,301	A→T	16.7%	intergenic (+30/+51)	PP_3628 → / ← PP_3629	choline/carnitine/betaine transporter family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,125,301	0	A	T	16.7%	80.0 / 9.8	36	intergenic (+30/+51)	PP_3628/PP_3629	choline/carnitine/betaine transporter family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base T (4/2);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

ATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGG  >  NC_002947/4125155‑4125444                                                                                                                                                   |                                                                                                                                                aTCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAag                                                                                                                                               >  4:132193/1‑149 (MQ=255) aTCGAGAAGGCCGGACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTa                                                                                                                                                 <  5:253189/147‑1 (MQ=255)    gagaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagca                                                                                                                                             <  4:305212/148‑1 (MQ=255)     agaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                           <  3:171177/149‑1 (MQ=255)     agaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                           <  7:93850/149‑1 (MQ=255)         ggCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGt                                                                                                                                       <  3:72976/149‑1 (MQ=255)          gCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGtg                                                                                                                                      <  6:274090/149‑2 (MQ=255)                cTGTAGTGCGCGGTGATGCATGGGGAGGAGTTCTACTTTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAATCCACGGGCTGCTAAGCAGCCCGTTCGCGa                                                                                                                                 <  6:15085/148‑1 (MQ=255)                     tGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCga                                                                                                                            <  4:105445/148‑1 (MQ=255)                                       ggaggaGTTCTACATTGTGGCGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGTAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACa                                                                                                         >  6:1307/1‑149 (MQ=255)                                        gaggagTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACAt                                                                                                        <  5:322495/149‑1 (MQ=255)                                                   cATTGTGGTGCAGGAGAAGAAGACATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCtt                                                                                              <  8:37145/148‑1 (MQ=255)                                                    aTTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGaaa                                                                                                           >  6:171267/1‑134 (MQ=255)                                                                gagaAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAAc                                                                                                          >  6:175223/1‑123 (MQ=255)                                                                gagaAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAAc                                                                                                          <  5:175223/123‑1 (MQ=255)                                                                               gCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGc                                                                 <  3:132193/149‑1 (MQ=255)                                                                                       aGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGAAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCg                                                          >  8:29801/1‑148 (MQ=255)                                                                                         cGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATc                                                       <  8:164465/149‑1 (MQ=255)                                                                                            cAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTg                                                    >  5:324079/1‑149 (MQ=255)                                                                                                         aaCCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGcatca                                       >  8:15136/1‑149 (MQ=255)                                                                                                           cccGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATg                                     <  7:29801/149‑1 (MQ=255)                                                                                                                  cTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCg                               >  5:158931/1‑148 (MQ=255)                                                                                                                     aCGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAac                            <  5:1307/148‑1 (MQ=255)                                                                                                                                ccGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAacgcacgc                      >  5:286757/1‑143 (MQ=255)                                                                                                                                ccGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAacgcacgc                      <  6:286757/143‑1 (MQ=255)                                                                                                                                ccGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgcc                >  7:83698/1‑149 (MQ=255)                                                                                                                                 cGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgccg               >  8:229750/1‑149 (MQ=255)                                                                                                                                 cGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgccg               >  1:87252/1‑149 (MQ=255)                                                                                                                                     aaCCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCcgccgc           <  6:324079/149‑1 (MQ=255)                                                                                                                                       gaaCGGGCTGCTTAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCACCCCCCCCCCTATCTTGGCCCCATCGATATTCCCAATCTGAAAATTGCCCTCGAACGCCCCATTGGcccccgag          >  2:339380/3‑146 (MQ=255)                                                                                                                                        cACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACg        >  6:136205/1‑149 (MQ=255)                                                                                                                                        cACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACg        >  4:217015/1‑149 (MQ=255)                                                                                                                                        aaCGGGCTGCTTAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACg        >  5:153292/2‑149 (MQ=255)                                                                                                                                         aCGGGCTGCTTAGCAGCCCGTTCGCGATCCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACgg       >  6:49622/1‑149 (MQ=255)                                                                                                                                          cGGGCTGCTTAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGt      >  8:263332/1‑149 (MQ=255)                                                                                                                                                 ctAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCggg  >  6:122851/1‑146 (MQ=255)                                                                                                                                                 ctAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCggg  <  5:122851/146‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                ATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGG  >  NC_002947/4125155‑4125444

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