Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,972,279	C→T	8.3%	A22A (GCC→GCT)	PP_5238 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,972,279	0	C	T	8.3%	134.1 / 3.2	49	A22A (GCC→GCT)	PP_5238	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/25);  new base T (2/2);  total (22/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.36e-01

AACAAGCCTTAAACGACGCTGTCCGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCTGGACAGTGCCAG  >  NC_002947/5972139‑5972416                                                                                                                                             |                                                                                                                                          aaCAAGCCTTAAACGACGCTGTCCGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGAt                                                                                                                                   <  4:122399/149‑1 (MQ=255)     aGCCTTAAACGACGCTGTCCGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGaa                                                                                                                               <  4:266434/149‑1 (MQ=255)            aaCGACGCTGTCCGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAAtcgtcg                                                                                                                        <  5:125409/149‑1 (MQ=255)            aaCGACGCTGTCCGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAAtcgtc                                                                                                                         <  2:162094/148‑1 (MQ=255)              cgacgCTGTCCGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGt                                                                                                                      <  6:271707/149‑1 (MQ=255)                      tCCGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGTGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAAtt                                                                                                              <  1:287199/149‑1 (MQ=255)                        cGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTc                                                                                                             <  7:103746/148‑1 (MQ=255)                              aCGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCa                                                                                                      <  2:274663/149‑1 (MQ=255)                                gtcattgtcattTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAg                                                                                                     >  4:169073/1‑148 (MQ=255)                                gtcattgtcattTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGc                                                                                                    <  1:126288/149‑1 (MQ=255)                                       tcattTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCa                                                                                             <  1:320003/149‑1 (MQ=255)                                          ttttCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATaa                                                                                                                 <  4:6627/126‑1 (MQ=255)                                          ttttCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATaa                                                                                                                 >  3:6627/1‑126 (MQ=255)                                             tCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCa                                                                                        <  7:90909/148‑1 (MQ=255)                                             tCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCa                                                                                        <  5:220225/148‑1 (MQ=255)                                             tCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAg                                                                                       >  6:309732/1‑149 (MQ=255)                                                gACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCTAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCa                                                                                    >  3:186947/1‑149 (MQ=255)                                                gACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCTAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCa                                                                                    >  3:190956/1‑149 (MQ=255)                                                gACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCa                                                                                    >  2:241422/1‑149 (MQ=255)                                                      ttGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGc                                                                                >  7:163157/1‑147 (MQ=255)                                                        gTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCaa                                                                              <  2:60957/147‑1 (MQ=255)                                                        gTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCaa                                                                              >  1:60957/1‑147 (MQ=255)                                                         ttACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATc                                                                                                                                  >  7:216883/1‑94 (MQ=255)                                                         ttACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATc                                                                                                                                  <  8:216883/94‑1 (MQ=255)                                                           aCCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCt                                                                                          >  7:150219/1‑132 (MQ=255)                                                           aCCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCt                                                                                          <  8:150219/132‑1 (MQ=255)                                                               gggAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCgg                                                                     <  5:309732/149‑1 (MQ=255)                                                                ggAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCggg                                                                    >  6:57117/1‑149 (MQ=255)                                                                    tGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCTAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGgc                                                                <  4:190956/149‑1 (MQ=255)                                                                     ggAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCTAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGgc                                                                <  4:186947/148‑1 (MQ=255)                                                                     ggAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGgc                                                                <  7:125852/148‑1 (MQ=255)                                                                      gAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGAGCTGGACCAGCAATACGAAAACTGGATCGAAATCGTCGGTTATAAATTCGGTTTCATCCAAAGCACATCAGCCACCGCAAGCTCTGCGGGTGgcg                                                               >  5:325170/1‑148 (MQ=255)                                                                            aTGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCa                                                                                                      >  8:34359/1‑103 (MQ=255)                                                                            aTGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCa                                                                                                      <  7:34359/103‑1 (MQ=255)                                                                            aTGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCtt                                                         >  5:45297/1‑148 (MQ=255)                                                                                cTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAggg                                                    <  6:180182/149‑1 (MQ=255)                                                                                 tttCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGgc                                                   >  3:144659/1‑149 (MQ=255)                                                                                              atCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCCACTCAGCCACCGCAAGCTCTGCCGGTGGCGAATCTTCAGGGCGCACCACCCacac                                      >  3:63915/1‑149 (MQ=25)                                                                                                   aaTCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCaccacc                                           <  3:14704/139‑1 (MQ=255)                                                                                                   aaTCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCaccacc                                           >  4:14704/1‑139 (MQ=255)                                                                                                                    gCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCTCCTTCACCAAAt                  <  2:24481/147‑1 (MQ=255)                                                                                                                    gCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAAt                  <  4:127923/147‑1 (MQ=255)                                                                                                                    gCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAAt                  <  4:132662/147‑1 (MQ=255)                                                                                                                      gAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCt              >  4:258896/1‑149 (MQ=255)                                                                                                                       aaGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCTg             >  2:70598/1‑149 (MQ=255)                                                                                                                         gcgcTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCTGGa           <  2:3314/149‑1 (MQ=255)                                                                                                                         gcgcTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCTGGa           >  1:3314/1‑149 (MQ=255)                                                                                                                                 ccAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCTGGACAGTGc     >  1:322800/1‑147 (MQ=255)                                                                                                                                  cAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCTGGACAGTGCCAg  >  2:197294/1‑149 (MQ=255)                                                                                                                                             cAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCTg             <  8:223927/127‑1 (MQ=255)                                                                                                                                             |                                                                                                                                          AACAAGCCTTAAACGACGCTGTCCGAGCAACGTCATTGTCATTTTCAGACCAATTGTTACCAGGGAGTGGAACACATGGCTTTCGATGCCTATATCCAAATCGCGGAAATTACTGGCGAAGCGCTGGACCAGCAATACGCCAACTGGATCGAAATCGTCGGTTATAAATTCGGTGTCAGCCAAAGCACCTCAGCCACCGCAAGCTCTGCGGGTGGCGCATCTTCAGGGCGCACCACCCTCACCAATTTCACCTTCACCAAATACCTGGACAGTGCCAG  >  NC_002947/5972139‑5972416

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