Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,311,315	T→A	7.6%	F12Y (TTT→TAT) ‡	PP_3784 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,311,315	0	T	A	7.6%	147.3 / 3.5	53	F12Y (TTT→TAT) ‡	PP_3784	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base A (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

ACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGG  >  NC_002947/4311169‑4311455                                                                                                                                                   |                                                                                                                                             acaACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTg                                                                                                                                            >  6:221993/1‑149 (MQ=255)   aaCCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCg                                                                                                                                          >  6:246245/1‑149 (MQ=255)    aCCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGt                                                                                                                                         >  8:91444/1‑149 (MQ=255)      cgcgTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTg                                                                                                                                       <  8:279691/149‑1 (MQ=255)        cgTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAAATGCGTTGcc                                                                                                                                     <  3:167674/149‑1 (MQ=255)             gcTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCg                                                                                                                                <  6:116716/149‑1 (MQ=255)             gcTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGCGCAACGCATTTGCGTTGCCCTTCg                                                                                                                                >  3:267705/1‑149 (MQ=255)                   aaaCCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTACCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGct                                                                                                                          >  6:277109/1‑149 (MQ=255)                     aCCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGctct                                                                                                                        >  3:246504/1‑149 (MQ=255)                      ccACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCtctc                                                                                                                       <  5:221993/149‑1 (MQ=255)                            gTGGCAAGGCGTTTCACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtca                                                                                                                   <  6:80580/147‑1 (MQ=255)                             tGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTACGCCGGCTCTCGtcatc                                                                                                                 >  5:299918/1‑148 (MQ=255)                             tGGCAAGGAGTTTACCCGACTTCCGTAACGCGCTGGCCCCCGAACCGATTCCGCGCCCCCCCGTCACCGACAGGAACATCATCACAGGCTAACAGCACTCCCCCCGGGGGCAACGCATTTGCGTTGCCCGGCGCCGGCTCTCGtcatca                                                                                                                <  2:308371/149‑1 (MQ=255)                                cAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAAtt                                                                                                             <  5:376927/149‑1 (MQ=255)                                 aaGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAAtt                                                                                                             >  2:944/1‑148 (MQ=255)                                   ggCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTTAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcat                                                                                                          <  7:91444/149‑1 (MQ=255)                                    gCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTGCGCCGGCTCTCGTCATCAATtcatc                                                                                                         <  3:54129/149‑1 (MQ=255)                                     cGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcatca                                                                                                        <  1:944/149‑1 (MQ=255)                                           aCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCa                                                                                                  >  7:206984/1‑149 (MQ=255)                                            ccGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGcc                                                                                                                                     <  6:29255/113‑1 (MQ=255)                                            ccGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGcc                                                                                                                                     >  5:29255/1‑113 (MQ=255)                                             cGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACg                                                                                                >  2:205440/1‑149 (MQ=255)                                             cGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACg                                                                                                >  6:383555/1‑149 (MQ=255)                                              gAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACgt                                                                                               <  4:267705/149‑1 (MQ=255)                                                     gTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGcc                                                                                        >  5:114999/1‑149 (MQ=255)                                                        aCGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGa                                                                                     <  7:31897/149‑1 (MQ=255)                                                        aCGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGa                                                                                     <  6:166826/149‑1 (MQ=255)                                                                  ccccGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGc                                                                           >  7:18030/1‑149 (MQ=255)                                                                                 gcgcCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCaa                                                            >  5:280334/1‑149 (MQ=255)                                                                                         ccGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTg                                                      <  4:246504/147‑1 (MQ=255)                                                                                          cGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTgccg                                                   <  5:383555/149‑1 (MQ=255)                                                                                                 gAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATc                                              >  7:143046/1‑147 (MQ=255)                                                                                                      gAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCCTTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGGGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCCAc                                       >  4:191852/1‑149 (MQ=255)                                                                                                      gAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAAc                                       <  8:18030/149‑1 (MQ=255)                                                                                                      gAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAAc                                       <  5:246245/149‑1 (MQ=255)                                                                                                              tcaATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCACTCATTGCCGCCATCCAGCAACGACCAGGc                               >  5:71783/1‑149 (MQ=255)                                                                                                                     cTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACtt                        >  7:114045/1‑149 (MQ=255)                                                                                                                        aCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTATGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGc                     >  4:330303/1‑149 (MQ=255)                                                                                                                        aCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGc                     >  3:241071/1‑149 (MQ=255)                                                                                                                        aCAGCATTCCCTCCGGGGGCAACGCATTTGCGATGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGc                     >  1:41575/1‑149 (MQ=255)                                                                                                                          aGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGt                   <  6:114999/149‑1 (MQ=255)                                                                                                                                cccTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtgctgc              <  2:190794/148‑1 (MQ=255)                                                                                                                                 ccTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCGGCAACGAACAGGCACCACTTGGCGTGCTGCGt            <  6:280334/149‑1 (MQ=255)                                                                                                                                 ccTCCGGGGGCAACGCATTTGCGTTGCCCTGCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGt            <  6:299918/149‑1 (MQ=255)                                                                                                                                  cTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTc           <  2:41575/149‑1 (MQ=255)                                                                                                                                        gggCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCg                    <  5:235814/134‑1 (MQ=255)                                                                                                                                        gggCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCg                    >  6:235814/1‑134 (MQ=255)                                                                                                                                          gCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCtacgagcg   >  7:117369/1‑142 (MQ=255)                                                                                                                                           cAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg  <  8:143046/149‑1 (MQ=255)                                                                                                                                           cAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg  >  6:122871/1‑149 (MQ=255)                                                                                                                                           cAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg  <  5:277109/149‑1 (MQ=255)                                                                                                                                            aaCGCATTTGCGTTGCCCTTCGCGGGCTATCGTCTTGAATTCATGAGTGCCACGTGTACTCCTGACACATTCATCGACTGCCTTTCTGATTCATTGCCGCCATCCAGCAACGAACTGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg  <  4:27595/148‑1 (MQ=255)                                                                                                                                            aaCGCATTTGCGTTGCCCTTCGCAGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGAGACTGCCTGTCGAAATCATTGCCGCCATCCAGCAACTAACAGGCACAACATGGCGTGCTGCTTCTTGTCTTgg  >  3:27595/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                             ACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGG  >  NC_002947/4311169‑4311455

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