Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,194,361	C→G	10.0%	R228G (CGG→GGG)	PP_4574 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,194,361	0	C	G	10.0%	107.2 / 4.9	40	R228G (CGG→GGG)	PP_4574	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base G (2/2);  total (18/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 = 8.19e-01

GCTTGGCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTACTACTTCCACAATTGCGTTTAT  >  NC_002947/5194216‑5194504                                                                                                                                                   |                                                                                                                                                gCTTGGCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAgcggc                                                                                                                                              <  1:359265/149‑1 (MQ=255)     ggCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGgcgc                                                                                                                                            <  5:27742/147‑1 (MQ=255)     ggCACCCGCAACATCAAGCTGGCCCTGCGCAGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCaa                                                                                                                                          >  3:78083/1‑149 (MQ=255)      gCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAac                                                                                                                                         <  1:54691/149‑1 (MQ=255)           cGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGt                                                                                                                                    <  4:118398/149‑1 (MQ=255)                 aTCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCa                                                                                                                                           >  1:69282/1‑137 (MQ=255)                 aTCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCa                                                                                                                                           <  2:69282/136‑1 (MQ=255)                 aTCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAgct                                                                                                                              >  7:341643/1‑149 (MQ=255)                    aaGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCGCTCGGGGCGCatctgga                                                                                                                                     <  5:188150/139‑7 (MQ=255)                    aaGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCGCTCGGGGCGCatctgga                                                                                                                                     >  6:188150/1‑133 (MQ=255)                          gCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCtgctgttgc                                                                                                                      <  2:10506/148‑1 (MQ=255)                             cTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTgctgttgctgtt                                                                                                                  <  3:283202/149‑1 (MQ=255)                                cgcCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCg                                                                                                                <  1:222773/148‑1 (MQ=255)                                   cGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACAt                                                                                                            >  3:328114/1‑149 (MQ=255)                                     gCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTg                                                                                                          <  1:144498/149‑1 (MQ=255)                                          gCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACa                                                                                                             <  8:47985/141‑1 (MQ=255)                                          gCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACa                                                                                                             >  7:47985/1‑141 (MQ=255)                                          gCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGt                                                                                                     >  5:178852/1‑149 (MQ=255)                                           cAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGtt                                                                                                    <  8:109355/149‑1 (MQ=255)                                               ttCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCGCTCGGGGCGCatctgg                                                                                                                                      <  7:38062/111‑6 (MQ=255)                                               ttCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCGCTCGGGGCGCatctgg                                                                                                                                      >  8:38062/1‑106 (MQ=255)                                                            ggCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGt                                                                                    <  4:78083/149‑1 (MQ=255)                                                                   ccgAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTc                                                                             <  8:187203/149‑1 (MQ=255)                                                                              gACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCAc                                                                                      >  8:217172/1‑129 (MQ=255)                                                                              gACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCAc                                                                                      <  7:217172/129‑1 (MQ=255)                                                                                cATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGtt                                                                                                    <  1:190115/113‑1 (MQ=255)                                                                                cATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGtt                                                                                                    >  2:190115/1‑113 (MQ=255)                                                                                    gATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTtctc                                                            >  4:350653/1‑149 (MQ=255)                                                                                        gCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGc                                                         <  6:178852/148‑1 (MQ=255)                                                                                             aTCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCaaa                                                                                <  4:156359/120‑1 (MQ=255)                                                                                                gACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGAc                                                >  5:17759/1‑149 (MQ=255)                                                                                                 aCCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGAcc                                               >  5:146285/1‑149 (MQ=255)                                                                                                          cgcgTGATTCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTTTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAg                                      <  5:149597/149‑1 (MQ=255)                                                                                                                    cGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAg                             >  1:28020/1‑148 (MQ=255)                                                                                                                          gttgAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGc                                                     <  2:241502/118‑1 (MQ=255)                                                                                                                          gttgAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGc                                                     >  1:241502/1‑118 (MQ=255)                                                                                                                            tgAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTtctc                                                            <  5:72551/109‑1 (MQ=255)                                                                                                                            tgAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTtctc                                                            >  6:72551/1‑109 (MQ=255)                                                                                                                             gAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTactact                   >  3:46747/1‑149 (MQ=255)                                                                                                                               aCATCCAGATGCGCCCCGAGCGGCGCAACTCCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGTCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTACTACtt                  <  8:341643/148‑1 (MQ=255)                                                                                                                                              ccGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTACTACTTCCACAATTGCGTTTAt  >  1:50636/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                GCTTGGCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAA‑GGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTACTACTTCCACAATTGCGTTTAT  >  NC_002947/5194216‑5194504

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