Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,158,876	A→C	9.0%	N296H (AAC→CAC)	PP_3661 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,158,876	0	A	C	9.0%	154.3 / 8.0	56	N296H (AAC→CAC)	PP_3661	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/21);  new base C (2/3);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

GTGGGCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGAT  >  NC_002947/4158729‑4159021                                                                                                                                                    |                                                                                                                                                  gtggGCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGaa                                                                                                                                                  >  2:274734/1‑149 (MQ=255)     gCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGc                                                                                                                                              >  4:253293/1‑149 (MQ=255)     gCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGc                                                                                                                                              <  6:266847/149‑1 (MQ=255)      ccTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGcc                                                                                                                                             >  3:199467/1‑149 (MQ=255)          tgGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTGCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCTGTTCACGGCCTgg                                                                                                                                          <  3:66225/148‑1 (MQ=255)           ggATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGa                                                                                                                                         >  2:140507/1‑148 (MQ=255)            gATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGaca                                                                                                                                       >  2:139363/1‑149 (MQ=255)             aTCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGacac                                                                                                                                      <  1:140507/149‑1 (MQ=255)                aGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCt                                                                                                                                    <  1:274734/148‑1 (MQ=255)                 gCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCttt                                                                                                                                  <  6:233846/149‑1 (MQ=255)                       cgcgCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTTGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCtt                                                                                                                             <  5:214783/148‑1 (MQ=255)                          gcAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGAc                                                                                                                          >  3:259602/1‑148 (MQ=255)                               gcaggcCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGc                                                                                                                    <  3:253293/149‑1 (MQ=255)                                                 aaCGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCAt                                                                                                  >  3:296780/1‑149 (MQ=255)                                                    gggCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGga                                                                                               <  4:296780/149‑1 (MQ=255)                                                        gcATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGAcc                                                                                           >  4:16359/1‑149 (MQ=255)                                                         cATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGAccc                                                                                          >  2:150424/1‑149 (MQ=255)                                                          aTTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCg                                                                                         >  2:181827/1‑149 (MQ=255)                                                          aTTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCg                                                                                         >  2:38574/1‑149 (MQ=255)                                                               cATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGac                                                                                                                                        <  6:235754/97‑1 (MQ=255)                                                               cATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCATGGCAGCGTGACCCAGGCCGTGAACAGCCTGGac                                                                                                                                        >  5:235754/1‑97 (MQ=255)                                                                  gCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACACGCTTCAGc                                                                                  >  8:259199/1‑148 (MQ=255)                                                                   cccTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGcc                                                                                 <  4:259602/148‑1 (MQ=255)                                                                    ccTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGc                                                                                  >  8:183666/1‑146 (MQ=255)                                                                         ttCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCg                                                                          >  4:59993/1‑149 (MQ=255)                                                                            gCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGcc                                                                        <  4:199467/148‑1 (MQ=255)                                                                              tttCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAgg                                                                     >  3:70170/1‑149 (MQ=255)                                                                                    cTTCCTGGCCCTGGTCCTGGTGAATGCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGc                                                               <  8:157900/149‑1 (MQ=255)                                                                                         tGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACgg                                                          >  4:180113/1‑149 (MQ=255)                                                                                             cctggtcctggtGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGgcgc                                                      >  7:12075/1‑149 (MQ=255)                                                                                                     tggtGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCAcc                                              >  1:245834/1‑149 (MQ=255)                                                                                                     ccgcccaTTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCAc                                               <  2:329069/142‑1 (MQ=255)                                                                                                      ggtGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCg                                             <  1:139363/149‑1 (MQ=255)                                                                                                                 tGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTc                                  <  7:183666/149‑1 (MQ=255)                                                                                                                    aaGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGTGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAAc                               >  4:177337/1‑149 (MQ=255)                                                                                                                    aaGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAAc                               >  2:179937/1‑149 (MQ=255)                                                                                                                    aaGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAAc                               >  5:103484/1‑149 (MQ=255)                                                                                                                      gTGCTGCCTGGCAGCGTGACCCTGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCt                             >  5:197452/1‑149 (MQ=255)                                                                                                                         ctgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgtg                          >  4:104144/1‑149 (MQ=255)                                                                                                                         ctgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgt                           <  8:12075/148‑1 (MQ=255)                                                                                                                         ctgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgt                           <  1:38574/148‑1 (MQ=255)                                                                                                                         ctgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCCGGCCACCGGTGCGATCCTCAACCtgt                           <  1:181827/148‑1 (MQ=255)                                                                                                                          tgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTc                                                                           <  5:256808/99‑1 (MQ=255)                                                                                                                          tgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTc                                                                           >  6:256808/1‑99 (MQ=255)                                                                                                                            ccTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCt                       >  8:208010/1‑149 (MQ=255)                                                                                                                                gCAGCGTGACCCTGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGtt                   <  6:197452/149‑1 (MQ=255)                                                                                                                                gCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCAGCGGTGCGATCCTCAACCTGTGGCtcgat                   >  4:215285/1‑145 (MQ=255)                                                                                                                                    aggtgtccAGGCTGTTCACGGCCTGGGCACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGt                                           <  6:236884/115‑1 (MQ=255)                                                                                                                                    aggtgtccAGGCTGTTCACGGCCTGGGCACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGt                                           >  5:236884/7‑121 (MQ=255)                                                                                                                                     gTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCgg              <  3:180113/149‑1 (MQ=255)                                                                                                                                         tccAGGCTGTTCACGGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGGCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgtg                          >  1:279084/2‑133 (MQ=255)                                                                                                                                         tccAGGCTGTTCACGGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgtg                          <  2:279084/132‑1 (MQ=255)                                                                                                                                         cccAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTgg                         <  4:223071/134‑1 (MQ=255)                                                                                                                                         cccAGGCCGTGAACAGCCTGGACACCTTTGCCTAGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTgg                         >  3:223071/1‑134 (MQ=255)                                                                                                                                           cAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCt        >  6:224823/1‑149 (MQ=255)                                                                                                                                            aGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTg       >  3:27840/1‑149 (MQ=255)                                                                                                                                                 gTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGAt  >  8:65032/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GTGGGCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGAT  >  NC_002947/4158729‑4159021

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