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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,405,492	0	A	T	7.3%	166.1 / 3.2	55	F23L (TTT→TTA) ‡	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 = 9.89e-01

GATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGCTTGTCAAT  >  NC_002947/4405345‑4405624                                                                                                                                                    |                                                                                                                                      gaTGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCa                                                                                                                                       >  3:218413/1‑148 (MQ=255)     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                                                                                                                     <  4:205944/120‑1 (MQ=255)                                               ccACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGa                                                                                                                     >  3:205944/1‑120 (MQ=255)                                                    cTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGt                                                                                  <  3:206660/149‑1 (MQ=255)                                                       aCTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGc                                                                                                      <  1:48951/126‑1 (MQ=255)                                                       aCTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGc                                                                                                      >  2:48951/1‑126 (MQ=255)                                                       aCTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTc                                                                                 <  2:181574/147‑1 (MQ=255)                                                       aCTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTc                                                                                 <  5:77334/147‑1 (MQ=255)                                                        cTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAg                                                                              <  7:14001/149‑1 (MQ=255)                                                         tGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGc                                                                             >  4:32490/1‑149 (MQ=255)                                                           ttACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCTTTGGCGaacagat                                                                                                                          >  2:266308/1‑98 (MQ=255)                                                           ttACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCTTTGGCGaacagat                                                                                                                          <  1:266308/103‑6 (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                                          <  4:76354/148‑1 (MQ=255)                                                                                             ctcCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGtt                                          >  3:76354/1‑148 (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)                                                                                                                                                    |                                                                                                                                      GATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCTTTTTTCGATCTGTTCGCCAAAGGCGAAATTATCACGATGGCGA‑AATAATGCCAAAACACTCGCCAATTCGTCAAGCGGCTTTTTCACGCTCTCTCATTCGGCCTACAAGTTGGGCAACAGGCCGCAGCGCTTGCTTGTCCAGCTTGTCAAT  >  NC_002947/4405345‑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: ‑