Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,774,877	Δ1 bp	9.1%	coding (635/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,877	0	C	.	9.1%	129.7 / 7.8	44	coding (635/1353 nt)	PP_5521	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/18);  new base . (2/2);  total (23/21)
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.58e-01

TCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCACTTGCGGGCCGCCCCAGCCATCGGTCACAC  >  NC_002947/2774736‑2775021                                                                                                                                                   |                                                                                                                                                  tccatgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGCGAGTACT‑GCTggtg                                                                                                                                            >  3:301434/1‑147 (MQ=255) tccatgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGCGAGTACT‑GCTggtg                                                                                                                                            >  3:296619/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)                                                                                                                                               tACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCACTTGCGGGCCGCCCCAGCCATCGGTCacac  <  2:187715/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAC‑GAAGATACC‑TTCGGCGCTCACC‑‑AGC‑AGTACTCGCTGGCCATCGGGGAA‑GGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCCTCCAGCCCCGTGGGCCAGGCCAGGCGTTGCACTTGCGGGCCGCCCCAGCCATCGGTCACAC  >  NC_002947/2774736‑2775021

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