Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,360,340	G→A	9.1%	A287V (GCA→GTA)	PP_2075 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,360,340	0	G	A	9.1%	115.5 / 3.0	44	A287V (GCA→GTA)	PP_2075	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (17/23);  new base A (2/2);  total (19/25)
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.30e-01

CGCTCAGGGTTGTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATAT  >  NC_002947/2360194‑2360475                                                                                                                                                   |                                                                                                                                        cGCTCAGGGTTGTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCaa                                                                                                                                      <  2:422288/150‑1 (MQ=38)     cAGGGTTGTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACttt                                                                                                                                  <  2:764207/150‑1 (MQ=40)     cAGGGTTGTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACttt                                                                                                                                  <  1:472019/150‑1 (MQ=255)          ttGTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAAca                                                                                                                             <  1:293415/150‑1 (MQ=255)          ttGTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTACAACTTTGAAca                                                                                                                             >  1:961375/1‑150 (MQ=255)          ttGTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTACAACTTTGAAca                                                                                                                             >  1:961274/1‑150 (MQ=255)            gTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAAcagc                                                                                                                           <  1:997374/150‑1 (MQ=255)               aCGGGTGGGCCGTGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAAcagccag                                                                                                                        >  2:1089552/1‑150 (MQ=255)               aCGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAAca                                                                                                                             <  1:1045456/145‑1 (MQ=255)               aCGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAAca                                                                                                                             >  2:1045456/1‑145 (MQ=255)                     gggACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCtt                                                                                                                  <  2:520798/150‑1 (MQ=40)                      ggACGGGCTTGCCTGGCAGGCCGGCCGCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGGCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCttt                                                                                                                 <  1:1361528/150‑1 (MQ=255)                       gACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGACGCCAGCCCAGCAGCAGGCTGAAACAATGGGCGGAGCCGATGAGGACGAACGCGGGCCAGACAAACAGCAGTGAGCCGGTGGACAGCCAGCTGCtttt                                                                                                                >  1:1142305/1‑150 (MQ=14)                         cTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGc                                                                                                              <  1:926925/150‑1 (MQ=255)                             cTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGgaa                                                                                                          >  1:1388271/1‑150 (MQ=255)                                         gccggccTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAg                                                                                              >  1:893985/1‑150 (MQ=255)                                         gccggccTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTACAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAg                                                                                              <  2:961375/150‑1 (MQ=255)                                         gccggccTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTACAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAg                                                                                              <  2:961274/150‑1 (MQ=255)                                           cggccTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGc                                                                                            <  1:776322/150‑1 (MQ=255)                                                tCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCg                                                                                       >  2:841772/1‑150 (MQ=255)                                                     gATGACCTTCGGAAAGCGCGGCAGCCCGCACAGCAGCAGGTGGAGTCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTg                                                                                  <  2:726305/150‑1 (MQ=14)                                                      aTGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTgg                                                                                 <  1:1089552/150‑1 (MQ=255)                                                        gACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGgc                                                                               <  2:617473/150‑1 (MQ=255)                                                        gACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGgc                                                                               <  2:1272959/150‑1 (MQ=255)                                                           cTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCt                                                                            >  1:1346256/1‑150 (MQ=255)                                                                       cgGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGAt                                                                <  2:1047463/150‑1 (MQ=255)                                                                              cagcacagcaGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATgggg                                                            >  1:321255/1‑146 (MQ=255)                                                                              cagcacagcaGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGt                                                            <  2:321255/147‑1 (MQ=255)                                                                              cagcacagcaGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAg                                                         >  2:925366/1‑150 (MQ=255)                                                                                  acagcaGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTa                                                     >  2:1386797/1‑150 (MQ=255)                                                                                        gcagGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTggg                                               <  2:42963/150‑1 (MQ=255)                                                                                          agGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTg                                             <  1:1331568/150‑1 (MQ=255)                                                                                               gAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCtt                                                                    >  1:364891/1‑122 (MQ=255)                                                                                               gAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCtt                                                                    <  2:364891/122‑1 (MQ=255)                                                                                                  ccAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAggcg                                     <  1:841772/150‑1 (MQ=255)                                                                                                   cAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTgg                                                                                 <  2:1389445/105‑1 (MQ=255)                                                                                                   cAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTgg                                                                                 >  1:1389445/1‑105 (MQ=255)                                                                                                    aGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAggcggg                                   <  2:893985/150‑1 (MQ=255)                                                                                                         gCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAggcgggggcgg                              <  1:872219/150‑1 (MQ=255)                                                                                                             tGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGCGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGcgggggcgggggg                          >  2:378469/1‑149 (MQ=255)                                                                                                                      tgaCGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGtgtg                 >  1:1064805/1‑150 (MQ=255)                                                                                                                      tgaCGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGtgtg                 >  1:1064787/1‑150 (MQ=255)                                                                                                                       gaCGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGATTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTgg                >  1:248884/1‑150 (MQ=255)                                                                                                                                     aGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGGAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCatat  >  1:1150354/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                        CGCTCAGGGTTGTCACGGGTGGGACTGGCTTGCCGGGCAGGCCGGCCTCGCGGATGACCTTCGGCCAGCGCGGTCGCCAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATAT  >  NC_002947/2360194‑2360475

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