Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,327,252	A→T	17.5%	N74Y (AAC→TAC) ‡	PP_3797 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,327,252	0	A	T	17.5%	87.1 / 12.8	40	N74Y (AAC→TAC) ‡	PP_3797	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (19/14);  new base T (2/5);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

TGGTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGAC  >  NC_002947/4327105‑4327397                                                                                                                                                    |                                                                                                                                                   tGGTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGa                                                                                                                                                    >  8:288444/1‑148 (MQ=255)   gTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAAc                                                                                                                                                  >  7:24781/1‑148 (MQ=255)      cTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAg                                                                                                                                               <  3:249054/148‑1 (MQ=255)         tGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGcc                                                                                                                                            <  8:24781/148‑1 (MQ=255)           tCCTGATCGTGTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                          <  1:109259/148‑3 (MQ=255)           tCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                          <  6:161983/148‑3 (MQ=255)                aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                          <  8:101568/143‑3 (MQ=255)                aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                          >  7:101568/1‑141 (MQ=255)                aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCATCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                          <  3:82588/143‑3 (MQ=255)                aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCATCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                          >  4:82588/1‑141 (MQ=255)                        tttCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑tgg                                                                                                                              <  1:200330/146‑1 (MQ=255)                          tCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑tggtggtg                                                                                                                         >  5:219453/1‑149 (MQ=255)                                     cTCGCCGTGGGCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGa                                                                                                               <  7:288444/148‑1 (MQ=255)                                      tCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGagc                                                                                                             >  5:167726/1‑149 (MQ=255)                                                    cagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGtt                                                                                                <  6:219453/148‑1 (MQ=255)                                                    cagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTc                                                                                               >  3:62674/1‑149 (MQ=255)                                                                  gCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTa                                                                                  >  7:145074/1‑148 (MQ=255)                                                                  gCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTa                                                                                  >  3:64297/1‑148 (MQ=255)                                                                       ccGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATAAGCTGTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGc                                                                             >  8:90872/1‑148 (MQ=255)                                                                        cGATCCGGCCGCACGCGCCCACAAGCCAGACGCGGACACCGAGAGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGCGCCGGCCA‑TGGTGGTGGAGGGCATGAGCCAGCCGGATGGTTCGCCCTTCGCGGGaactgcc                                                                            >  2:295204/1‑143 (MQ=255)                                                                        cGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCa                                                                           <  4:62674/149‑1 (MQ=255)                                                                                          caacaaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAAt                                                                         >  8:235076/1‑133 (MQ=255)                                                                                          caacaaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGAct                                                         >  1:112682/1‑148 (MQ=255)                                                                                              gaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAc                                                     <  5:68448/148‑1 (MQ=255)                                                                                                gCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAc                                                     >  6:102547/1‑147 (MQ=255)                                                                                                 ccAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAg                                                  <  8:145074/149‑1 (MQ=255)                                                                                                         ggACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGAt                                                                                                    >  7:187254/1‑91 (MQ=255)                                                                                                         ggACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGAt                                                                                                    <  8:187254/91‑1 (MQ=255)                                                                                                         ggACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGc                                          <  8:6270/149‑1 (MQ=255)                                                                                                         ggACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGc                                          <  6:167726/149‑1 (MQ=255)                                                                                                             acCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCtgta                                      <  2:112682/149‑1 (MQ=255)                                                                                                                      aTGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAAacc                             >  7:68273/1‑149 (MQ=255)                                                                                                                            tttCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTacac                       >  2:289850/1‑149 (MQ=255)                                                                                                                              tCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGAACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCg                     >  8:296123/1‑149 (MQ=255)                                                                                                                                agcagGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGa                   <  7:90872/149‑1 (MQ=255)                                                                                                                                          ggAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAgg         >  7:125988/1‑149 (MQ=255)                                                                                                                                           gAGGCCTTGAACAAGGCCGAGCCGGACAGTGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAgg         >  1:179362/1‑149 (MQ=255)                                                                                                                                             ggCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCa      <  2:179362/149‑1 (MQ=255)                                                                                                                                                 ttGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGAc  >  2:186990/1‑149 (MQ=255)                                                                                                                                                   gAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCa           >  6:212946/1‑138 (MQ=255)                                                                                                                                                   gAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCa           <  5:212946/138‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TGGTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGAC  >  NC_002947/4327105‑4327397

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