Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,138,588	G→C	9.8%	S29S (TCG→TCC) ‡	PP_3643 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,138,588	0	G	C	9.8%	145.7 / 6.0	49	S29S (TCG→TCC) ‡	PP_3643	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base C (2/3);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGG  >  NC_002947/4138441‑4138725                                                                                                                                                    |                                                                                                                                          aCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCgg                                                                                                                                          <  4:404583/149‑1 (MQ=255)  ccATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCgg                                                                                                                                          <  4:130224/148‑1 (MQ=255)   cATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCAACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtc                                                                                                                                        >  4:32618/1‑149 (MQ=255)        gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCgg                                                                                                                                   >  1:58825/1‑149 (MQ=255)           acaTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTCCt                                                                                                                                <  8:22415/149‑1 (MQ=255)             aTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGa                                                                                                                              >  4:238147/1‑149 (MQ=255)                aaCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGtcgtcgtc                                                                                                                                  >  1:313695/1‑140 (MQ=255)                aaCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGtcgtcgtc                                                                                                                                  <  2:313695/142‑3 (MQ=255)                    aTAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAg                                                                                                                       >  7:318417/1‑149 (MQ=255)                     tAAAACCCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAgc                                                                                                                      <  3:279112/149‑1 (MQ=255)                     tAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCTATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAgc                                                                                                                      <  3:273825/149‑1 (MQ=255)                     tAAAAACCGCGCACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCAACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAgc                                                                                                                      <  3:32618/149‑1 (MQ=255)                           ccGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtc                                                                                                                <  2:135913/149‑1 (MQ=255)                           ccGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtc                                                                                                                <  3:238147/149‑1 (MQ=255)                                           cccTTCAGAGGTACTCCCGATGATCGACGCAAGTTTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATgg                                                                                                <  7:67520/149‑1 (MQ=255)                                             cTTCAGAGGTACTCCCGATGATCGACGAAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAc                                                                                              <  8:318417/149‑1 (MQ=255)                                                 agagGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACccgc                                                                                          <  6:308938/149‑1 (MQ=255)                                                     gTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCt                                                                                      >  8:377741/1‑149 (MQ=255)                                                      tACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCGCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTg                                                                                     <  6:300611/149‑1 (MQ=255)                                                        cTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                  <  3:394618/102‑1 (MQ=255)                                                        cTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                  >  4:394618/1‑102 (MQ=255)                                                           ccGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGcc                                                                                                                   <  4:79472/114‑1 (MQ=255)                                                           ccGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGcc                                                                                                                   >  3:79472/1‑114 (MQ=255)                                                                     aCGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGcccc                                                                      <  7:377741/149‑1 (MQ=255)                                                                          aCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGAt                                                                                                  <  5:121810/116‑1 (MQ=255)                                                                          aCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGAt                                                                                                  >  6:121810/1‑116 (MQ=255)                                                                          aCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTa                                                                 <  3:350896/149‑1 (MQ=255)                                                                              tCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACAcc                                                             <  2:58825/149‑1 (MQ=255)                                                                                   aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCCCCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGa                                                        >  1:359840/1‑149 (MQ=255)                                                                                   aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAccc                                                                                            >  7:213644/1‑113 (MQ=255)                                                                                   aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAccc                                                                                            <  8:213644/113‑1 (MQ=255)                                                                                     gAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACt                                                      >  7:15048/1‑149 (MQ=255)                                                                                     gAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACt                                                      >  2:145012/1‑149 (MQ=255)                                                                                          ttCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTcc                                                          >  3:272580/1‑140 (MQ=255)                                                                                          ttCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTcc                                                          <  4:272580/140‑1 (MQ=255)                                                                                               ggCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCgg                                             <  8:15048/148‑1 (MQ=255)                                                                                                    cTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTAcc                                                >  2:211924/1‑140 (MQ=255)                                                                                                    cTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTAcc                                                <  1:211924/140‑1 (MQ=255)                                                                                                        ccGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCaa                                   <  2:359840/149‑1 (MQ=255)                                                                                                                 gTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAgcg                          >  3:108894/1‑149 (MQ=255)                                                                                                                     ccGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCACCAACTACACCTcc                                                          >  8:381404/1‑113 (MQ=255)                                                                                                                        tcatcaCCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGc                    >  5:247721/1‑148 (MQ=255)                                                                                                                        tcatcaCCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGc                    >  4:45084/1‑148 (MQ=255)                                                                                                                             accaccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCaa                                   <  7:2840/128‑1 (MQ=255)                                                                                                                             accaccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCaa                                   >  8:2840/1‑128 (MQ=255)                                                                                                                               caccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCAcc            >  6:307232/1‑149 (MQ=255)                                                                                                                                accCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACct           <  3:45084/149‑1 (MQ=255)                                                                                                                                 cccAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACctg          <  1:145012/149‑1 (MQ=255)                                                                                                                                  tcAGGCCGACGACCGAGCCGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                           >  4:312203/2‑115 (MQ=25)                                                                                                                                  tcAGGCCGACGACCGAGCCGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                           <  3:312203/114‑1 (MQ=25)                                                                                                                                         gacgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTTTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCgg  >  6:225932/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                          ACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGG  >  NC_002947/4138441‑4138725

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