Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,196,645	G→T	9.5%	R231S (CGC→AGC)	PP_4577 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,196,645	0	G	T	9.5%	115.9 / 3.7	42	R231S (CGC→AGC)	PP_4577	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base T (2/2);  total (20/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 = 2.77e-01

GGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGCTGGCAGTGTATG  >  NC_002947/5196502‑5196791                                                                                                                                                |                                                                                                                                                   ggTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGggcggc                                                                                                                                              >  2:698307/1‑150 (MQ=255)  gTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCa                                                                                                                                             >  2:216964/1‑150 (MQ=255)       aCAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTg                                                                                                                                        <  1:216964/150‑1 (MQ=255)         aGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCa                                                                                                                                      <  1:251599/150‑1 (MQ=255)              tACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGgcg                                                                                                                                 >  1:706490/1‑150 (MQ=255)                       gCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCACCATGCACGCAGAGGGTGCGggcggc                                                                                                                                              <  2:191816/128‑1 (MQ=255)                       gCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCACCATGCACGCAGAGGGTGCGggcggc                                                                                                                                              >  1:191816/1‑128 (MQ=255)                         tGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTc                                                                                                                      <  1:210617/150‑1 (MQ=255)                            cGATGGCGCCCACTTCGATTCGCGGCTTCATCACCCCTCCAGCCTTTCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTAGCCATGCACGCAGAGGGGGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGc                                                                                                                   <  2:632257/150‑1 (MQ=255)                             gATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGcc                                                                                                                  >  1:741575/1‑150 (MQ=255)                                             aTTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCTATCGTCAGCCAACAGCGTTTCGCCTc                                                                                                  <  2:10092/150‑1 (MQ=255)                                              ttCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCTGCCAACAGCGTTTCGCCTCg                                                                                                 <  2:404862/150‑1 (MQ=255)                                                 gccggttCATCACACCTTCCGCCCATGGAAGGCCTTTCGAATATCCCGCACCGCTCCCACCGAACTGTCGTTATCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCTCTGCCATCGTCATCCACCAGCGTTTCGCCTCGGGc                                                                                              >  2:61081/6‑150 (MQ=255)                                                    gCTTCATCGCACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAg                                                                                           >  2:218897/1‑150 (MQ=255)                                                        catcaCACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGc                                                                                       <  2:465833/150‑1 (MQ=255)                                                          tcaCACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGCCCTGGCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCAc                                                                                     >  1:656028/1‑150 (MQ=255)                                                                 tCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCAcc                                                                                    >  1:212451/1‑144 (MQ=255)                                                                 tCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCAcc                                                                                    <  2:212451/144‑1 (MQ=255)                                                                    cgCCCATCGAGGGCCTGGAGAATCTTCAGCCCCGCCCCCGACGAAATGTCGGTGTGGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACCGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTcc                                                                           <  2:163589/149‑1 (MQ=255)                                                                     gCCCACCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGGGCGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTcca                                                                          <  1:424959/150‑1 (MQ=255)                                                                            gAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCcg                                                                   >  1:306216/1‑150 (MQ=255)                                                                             aGGGCCTGGCAAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCcgc                                                                  >  2:474029/1‑150 (MQ=255)                                                                                  cTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCgg                                                             <  1:119075/150‑1 (MQ=255)                                                                                       gAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGt                                                        <  1:207942/150‑1 (MQ=255)                                                                                         aTCGGCCGAACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGa                                                      <  2:274626/150‑1 (MQ=255)                                                                                          tCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGAt                                                     >  1:148794/1‑150 (MQ=255)                                                                                               cGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCAcc                                                <  2:706490/150‑1 (MQ=255)                                                                                                     gccgccACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGccc                                          <  2:656028/150‑1 (MQ=255)                                                                                                               aaCTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCa                                                                                                                           >  2:208352/1‑59 (MQ=255)                                                                                                               aaCTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCa                                                                                                                           <  1:208352/59‑1 (MQ=255)                                                                                                                    tCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGAc                           >  1:531582/1‑150 (MQ=255)                                                                                                                       ttGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGAcagc                        >  1:525864/1‑150 (MQ=255)                                                                                                                           cGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGAt                    <  1:474029/150‑1 (MQ=255)                                                                                                                               aTGCACGCAGAGGGTGCTGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGcc                >  1:162568/1‑150 (MQ=255)                                                                                                                               aTGCACGCAGAGGGTGCTGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGcc                >  1:162489/1‑150 (MQ=255)                                                                                                                                 gCACGCAGAGGGTGCTGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGc              <  2:162568/150‑1 (MQ=255)                                                                                                                                 gCACGCAGAGGGTGCTGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGc              <  2:162489/150‑1 (MQ=255)                                                                                                                                 gCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACTGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGc              <  1:518976/150‑1 (MQ=255)                                                                                                                                   aggCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCAACCTCGGTGCTGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGCTg            <  2:255936/148‑1 (MQ=255)                                                                                                                                   aCGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGCTg            >  1:624500/1‑150 (MQ=255)                                                                                                                                     gCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGACAGCCGACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGCTGGc          >  1:717957/1‑150 (MQ=255)                                                                                                                                             tGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGAGATTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGCTGGCAGTGTATg  >  1:736618/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                   GGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCCGCCACCGAACTGTCGTTGTCGCCATGCACGCAGAGGGTGCGGGCGGCAAGCTGCAAGGCGCTGCCATCGTCAGCCAACAGCGTTTCGCCTCGGGCCAGACGCACCGCCTGCTCCACCACCCGCGCCGGGTCGTGATGCACCGCGCCCGGCAGGCGCCGCGACAGCAGATGGCCGCTGGCAGTGTATG  >  NC_002947/5196502‑5196791

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