Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,405,493	A→T	7.3%	F23Y (TTT→TAT) ‡	PP_3886 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,405,493	0	A	T	7.3%	166.7 / 3.0	55	F23Y (TTT→TAT) ‡	PP_3886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base T (2/2);  total (28/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 = 8.33e-01

GCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGCTTGTCAAT  >  NC_002947/4405349‑4405624                                                                                                                                                 |                                                                                                                                     gCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGc                                                                                                                                  <  4:250043/149‑1 (MQ=255)      ccTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCTTTGGCGaa                                                                                                                               >  6:33362/1‑147 (MQ=255)          tCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCTTTGGCGaacagat                                                                                                                          <  5:33362/148‑6 (MQ=255)            aGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCAc                                                                                                                       >  6:77334/1‑149 (MQ=255)              gTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGa                                                                                                                     <  8:85332/149‑1 (MQ=255)               tCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGaaa                                                                                                                              >  4:202293/1‑139 (MQ=255)               tCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGaaa                                                                                                                              <  3:202293/139‑1 (MQ=255)                cGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATg                                                                                                                   <  4:218413/149‑1 (MQ=255)                          gcgcTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑aataat                                                                                                        <  8:249809/149‑1 (MQ=255)                                   tCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAaca                                                                                               >  1:181574/1‑149 (MQ=255)                                      aGTCCCCACGCTGACTGTTACCCCTGCACGTCGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTc                                                                                            <  7:15374/149‑1 (MQ=255)                                         ccccACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGc                                                                                          >  4:206660/1‑148 (MQ=255)                                          cccACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGcc                                                                                         >  2:178862/1‑148 (MQ=255)                                           ccACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGa                                                                                                                     >  3:205944/1‑120 (MQ=255)                                           ccACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGa                                                                                                                     <  4:205944/120‑1 (MQ=255)                                                cTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGt                                                                                  <  3:206660/149‑1 (MQ=255)                                                   aCTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGc                                                                                                      >  2:48951/1‑126 (MQ=255)                                                   aCTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGc                                                                                                      <  1:48951/126‑1 (MQ=255)                                                   aCTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTc                                                                                 <  5:77334/147‑1 (MQ=255)                                                   aCTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTc                                                                                 <  2:181574/147‑1 (MQ=255)                                                    cTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAg                                                                              <  7:14001/149‑1 (MQ=255)                                                     tGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGc                                                                             >  4:32490/1‑149 (MQ=255)                                                       ttACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCTTTGGCGaacagat                                                                                                                          <  1:266308/103‑6 (MQ=255)                                                       ttACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCTTTGGCGaacagat                                                                                                                          >  2:266308/1‑98 (MQ=255)                                                           cccTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCttt                                                                       >  7:147380/1‑149 (MQ=255)                                                           cccTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCttt                                                                       >  3:113562/1‑149 (MQ=255)                                                           cccTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACGCGCCAATTCGTCAAGCGGCttt                                                                       >  3:242752/1‑149 (MQ=255)                                                                       ccgccgGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGctctc                                                            <  4:113562/148‑1 (MQ=255)                                                                       ccgccgGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGctctc                                                            <  3:32490/148‑1 (MQ=255)                                                                                         ctcCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGtt                                          >  3:76354/1‑148 (MQ=255)                                                                                         ctcCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGtt                                          <  4:76354/148‑1 (MQ=255)                                                                                          tcCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTgg                                        >  6:310004/1‑149 (MQ=255)                                                                                            cAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTggg                                       >  2:18406/1‑148 (MQ=255)                                                                                              aaaaCGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGtt                                          >  4:196950/1‑143 (MQ=255)                                                                                              aaaaCGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGtt                                          <  3:196950/143‑1 (MQ=255)                                                                                                     cAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCTTTCGGCCTACAAGTTGGGCAACAGGCCg                             <  1:18406/149‑1 (MQ=255)                                                                                                      aaaGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCg                             >  1:118829/1‑148 (MQ=255)                                                                                                      aaaGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGc                            >  7:92683/1‑149 (MQ=255)                                                                                                        aGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAg                          >  3:58165/1‑149 (MQ=255)                                                                                                        aGCCCAGCTCGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAg                          >  3:84107/1‑149 (MQ=255)                                                                                                            cAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGc                                                   <  2:194383/120‑1 (MQ=255)                                                                                                            cAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATAACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGc                                                   >  1:194383/1‑120 (MQ=255)                                                                                                                      gCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTc              >  5:123505/1‑147 (MQ=255)                                                                                                                        cGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGATAATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTc              >  1:232865/1‑146 (MQ=255)                                                                                                                        cGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGc          <  2:106644‑M2/149‑2 (MQ=255)                                                                                                                        cGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGc          >  4:84107‑M2/1‑148 (MQ=255)                                                                                                                           gCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCgcttgc                   >  6:67105/1‑137 (MQ=255)                                                                                                                           gCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCgcttgc                   <  5:67105/137‑1 (MQ=255)                                                                                                                           gCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGcttg       >  5:310004‑M2/1‑145 (MQ=255)                                                                                                                            cTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGcttgt      <  6:90315‑M2/149‑6 (MQ=255)                                                                                                                             ttttttCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGctggacaagcaagcgctgcggcctgttgcccaacttg       >  6:203529‑M1/1‑110 (MQ=255)                                                                                                                             ttttttCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGcttgt      <  5:268391‑M2/148‑6 (MQ=255)                                                                                                                             ttttttCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACAATCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCAGCAGCGCTTGCTTGTCCAGcttgt      <  6:212052‑M2/148‑6 (MQ=255)                                                                                                                               ttttCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGctggacaagcaagcgctgcggcctgttgcccaacttgtag    <  5:203529‑M1/148‑41 (MQ=255)                                                                                                                                 ttCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGCTTGTCAat  <  6:123505/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                     GCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGCTTGTCAAT  >  NC_002947/4405349‑4405624

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