Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,774,870:1	+G	9.3%	coding (642/1353 nt)	PP_5521 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,774,870	1	.	G	9.3%	155.8 / 12.2	43	coding (642/1353 nt)	PP_5521	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (21/18);  new base G (2/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCACTTGCGGGCCGCCCCAGCCATC  >  NC_002947/2774724‑2775013                                                                                                                                                        |                                                                                                                                                 acCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACc‑‑agc‑ag                                                                                                                                                <  2:181087/149‑1 (MQ=255)             tccatgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGCGAGTACT‑GCTggtg                                                                                                                                    >  3:296619/1‑147 (MQ=255)             tccatgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGCGAGTACT‑GCTggtg                                                                                                                                    >  3:301434/1‑147 (MQ=255)               catgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGCGAGTACT‑GCTggtgag                                                                                                                                  <  4:301434/149‑5 (MQ=255)               catgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGCGAGTACT‑GCTggtgag                                                                                                                                  <  4:296619/149‑5 (MQ=255)                   ctcAGCCCCAGGCTGATGTCTTCCGGGTCTGTGCGCTCGGCCAGCTCTTGGTGGACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCAGCTTGAGTGGCCAGTAG‑GATGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTGGCTGGCCATCgg                                                                                                                               <  8:186196/148‑1 (MQ=255)                    tcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCggg                                                                                                                              >  3:68016/1‑148 (MQ=255)                           cAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTgg                                                                                                                                      <  8:227371/133‑1 (MQ=255)                           cAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTgg                                                                                                                                      >  7:227371/1‑133 (MQ=255)                              gCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCa                                                                                                                  <  3:262745/149‑1 (MQ=255)                                tGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAgc                                                                                                                >  7:6394/1‑149 (MQ=255)                                  aTGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGACTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAgc                                                                                                                <  3:240460/147‑1 (MQ=255)                                     tcgtcCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGcggt                                                                                                           >  7:223817/1‑149 (MQ=255)                                                 gTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGAGCGGTCGGTGTGggcgg                                                                                               >  7:209987/1‑149 (MQ=255)                                                       tCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACCTTTTTTAAATGACCATAGC‑AGAACACGCGGGCCAGCGGGGAA‑GGGGGTCAGCGCGGGCGGTGTGggcggcgg                                                                                            >  1:101609/1‑147 (MQ=255)                                                         ggCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGGCGGTGTGGGCGGCAGGTc                                                                                         <  6:115185/147‑1 (MQ=255)                                                            cAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACg                                                                                    <  8:6394/149‑1 (MQ=255)                                                               ctctTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCgg                                                                                                                               >  8:12952/1‑104 (MQ=255)                                                               ctctTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCgg                                                                                                                               <  7:12952/104‑1 (MQ=255)                                                                  ttCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTc                                                                                <  8:209987/147‑1 (MQ=255)                                                                  ttCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCg                                                                               >  1:35249/1‑148 (MQ=255)                                                                      aGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACgg                                                                          <  2:35249/149‑1 (MQ=255)                                                                        cACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGt                                                                        <  2:233373/149‑1 (MQ=255)                                                                           gcgTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑gg                                                                                                                        >  6:79486/1‑98 (MQ=255)                                                                           gcgTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAAc                                                                     >  1:208706/1‑149 (MQ=255)                                                                             gTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCt                                                                   >  3:89379/1‑149 (MQ=255)                                                                              cgCCTGCTGCCACAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGGAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGggc                                                                                                 >  1:69829/2‑119 (MQ=255)                                                                              cgCCTGCTGCCACAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGggc                                                                                                 <  2:69829/117‑1 (MQ=255)                                                                                                  gAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGcccc                                              <  2:208706/149‑1 (MQ=255)                                                                                                       gCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTggg                                         >  6:266199/1‑149 (MQ=255)                                                                                                           cTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGggccag                                     >  7:211235/1‑149 (MQ=255)                                                                                                           cTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCAGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTggc                                                                <  7:4524/122‑1 (MQ=255)                                                                                                           cTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCAGAGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTggc                                                                >  8:4524/1‑122 (MQ=255)                                                                                                              gAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGggccaggc                                   <  2:101609/148‑1 (MQ=255)                                                                                                                    ccAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGt                            >  8:32808/1‑149 (MQ=255)                                                                                                                     cAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGtt                           <  4:89379/149‑1 (MQ=255)                                                                                                                        tAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCa                        >  3:286415/1‑149 (MQ=255)                                                                                                                               gATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCACTTGCg                  <  4:68016/149‑1 (MQ=255)                                                                                                                                    c‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCACTTGCGGgccg             >  1:65696/1‑149 (MQ=255)                                                                                                                                    c‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCACGTGGGCCAGGCCAGGCGTTGCACTTGCgg                 >  6:206579/1‑145 (MQ=255)                                                                                                                                    c‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCACGTGGGCCAGGCCAGGCGTTGCACTTGCgg                 <  5:206579/145‑1 (MQ=255)                                                                                                                                                aCC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAAGGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGggccaggcca                                 >  1:38008/1‑119 (MQ=255)                                                                                                                                                aCC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGggccaggcca                                 <  2:38008/118‑1 (MQ=255)                                                                                                                                                aCC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCACTTGCGGGCCGCCCCAGCCATc  >  3:7661/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                 ACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCACTTGCGGGCCGCCCCAGCCATC  >  NC_002947/2774724‑2775013

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑