Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,732,030	A→G	7.0%	L404S (TTG→TCG)	PP_3299 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,732,030	0	A	G	7.0%	160.4 / 2.3	57	L404S (TTG→TCG)	PP_3299	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base G (2/2);  total (31/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 = 9.42e-01

ATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCAAAGAGTCTCCACCTCAAACCGAG  >  NC_002947/3731884‑3732175                                                                                                                                                   |                                                                                                                                                  aTGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAAt                                                                                                                                                 >  2:55337/1‑149 (MQ=255)    tCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGc                                                                                                                                               <  1:55337/148‑1 (MQ=255)       tGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGa                                                                                                                                            >  4:181888/1‑148 (MQ=255)       tGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGa                                                                                                                                            >  5:353701/1‑148 (MQ=255)          gaAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACaaa                                                                                                                                        <  6:353701/149‑1 (MQ=255)                    aTCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCCCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCa                                                                                                                               >  5:214946/1‑148 (MQ=255)                       tGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACa                                                                                                                            >  3:274853/1‑148 (MQ=255)                       tGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACAt                                                                                                                           >  3:349052/1‑149 (MQ=255)                        gCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTAGGAGGGACGAACGTCGTGCTCGCTCCAGGGATACATCACCGCGTTTGAGCGACTGTGCCTGCGCTTGACGATTGATAAATGAAGACAAAGCGGTTTCAACATc                                                                                                                          >  8:309701/1‑149 (MQ=255)                        gCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCGGATTGATCAAt                                                                                                                                                 <  2:326795/126‑1 (MQ=255)                        gCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCGGATTGATCAAt                                                                                                                                                 >  1:326795/1‑126 (MQ=255)                        gCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATc                                                                                                                          >  4:335080/1‑149 (MQ=255)                        gCGTTCAGAACGTCTAGCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATc                                                                                                                          >  2:178371/1‑149 (MQ=255)                         cGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCaa                                                                                                                              <  3:181888/144‑1 (MQ=255)                          gTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCtt                                                                                                                        <  3:335080/149‑1 (MQ=255)                                           aaCTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGca                                                                                                         >  3:347399/1‑147 (MQ=255)                                           aaCTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTCTCAACATCTTCCGaa                                                                                                                   >  7:89060/1‑137 (MQ=255)                                           aaCTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGGCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTCTCAACATCTTCCGaa                                                                                                                   <  8:89060/137‑1 (MQ=255)                                            aCTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGca                                                                                                         >  5:85118/1‑146 (MQ=255)                                             cTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGa                                                                                                     <  4:274853/149‑1 (MQ=255)                                                         ttCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGAt                                                                                         <  2:129999/149‑1 (MQ=255)                                                           cTGCAATGCTCCGATGGACGAACGTGGTGCTGGCTCCAGGGATACCTCACCGCGTTTGCGCGAGTGTGCGTGCGCTTGCCGATTTACCAATGCAGACAACGGTGTTGCAACATCTTCCGAAGCCCCCAGCACAGACTGCCTATAGATg                                                                                        <  6:312581/148‑1 (MQ=255)                                                            tGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGc                                                                                      >  4:205790/1‑149 (MQ=255)                                                            tGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGc                                                                                      >  2:313748/1‑149 (MQ=255)                                                            tGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGAGTTTGAGCGACTGTGCCTCCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGc                                                                                      >  6:333708/1‑149 (MQ=255)                                                              gAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCCCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCa                                                                                    <  6:214946/149‑1 (MQ=255)                                                               aaTGCTAGGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTTCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCa                                                                                    <  7:164673/148‑1 (MQ=255)                                                               aaTGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAg                                                                                   >  5:52126/1‑149 (MQ=255)                                                                          ggACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGggcttcggctt                                                                        <  4:347399/149‑1 (MQ=255)                                                                            aCGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCa                                                                      <  6:52126/149‑1 (MQ=255)                                                                            aCGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCa                                                                      >  7:69242/1‑149 (MQ=255)                                                                                     gtGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCttt                                                              >  4:70532/1‑148 (MQ=255)                                                                                     gtGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCttt                                                              <  3:70532/148‑1 (MQ=255)                                                                                         tcgctcCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTTCCGATAGATTGACAGGGCTTCGGCTTCACTTCCTTTggaat                                                         >  3:66688/1‑146 (MQ=255)                                                                                           gctcCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTgg                                                       <  3:205790/149‑1 (MQ=255)                                                                                           gctcCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTgg                                                       >  1:42999/1‑149 (MQ=255)                                                                                                        cTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGggcttcggct                                                                         >  5:194224/1‑118 (MQ=255)                                                                                                        cTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGggcttcggct                                                                         <  6:194224/118‑1 (MQ=255)                                                                                                             cgcgTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCGATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATg                                     >  1:328002/1‑149 (MQ=255)                                                                                                             cgcgTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATg                                     <  2:42999/149‑1 (MQ=255)                                                                                                                 tttGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGc                                 <  6:85118/149‑1 (MQ=255)                                                                                                                 tttGAGCGACTGCGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTg                                                        <  4:131749/126‑1 (MQ=255)                                                                                                                 tttGAGCGACTGAGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTg                                                        >  3:131749/1‑126 (MQ=255)                                                                                                                             tgCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCaaa                       <  4:349052/147‑1 (MQ=255)                                                                                                                              gCCTGCGCTTGCCGATTGATCGATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTg                                                        <  3:245870/113‑1 (MQ=255)                                                                                                                              gCCTGCGCTTGCCGATTGATCGATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTg                                                        >  4:245870/1‑113 (MQ=255)                                                                                                                                cTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGGTGATGCGGCTAACATCAAAGAGtc                  >  8:105586/1‑149 (MQ=255)                                                                                                                                cTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGACTGGTCTATCAGCGCGTTGATGCGGCTAATATCAAAGAGtc                  >  6:351591/1‑149 (MQ=255)                                                                                                                                    gcTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCaaa                       >  3:34956/1‑140 (MQ=255)                                                                                                                                    gcTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCaaa                       <  4:34956/140‑1 (MQ=255)                                                                                                                                      ttGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCAAAGAGtctc                >  3:130896/1‑145 (MQ=255)                                                                                                                                      ttGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCAAAGAGtctc                <  4:130896/145‑1 (MQ=255)                                                                                                                                        gacGATTTATCAATTCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTACAGCACAGACTGCCTATAGATGGCCAGTGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCAAAGAGTCTCCACCt           <  4:66688/146‑1 (MQ=255)                                                                                                                                         ccGATTGATCGATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCAAAGAGTCTCCACCTc          <  2:328002/148‑1 (MQ=255)                                                                                                                                                aTCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGCTGGGGCTAAAATCTAAGAGTCTCCACCTCAACCCgcg  >  6:258618/1‑147 (MQ=255)                                                                                                                                                  cAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCAAAGAGTCTCCACCTCAAACCGAg  >  5:146593/1‑147 (MQ=255)                                                                                                                                                  cAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCAAAGAGTCTCCACCTCAAACCGAg  <  6:146593/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTTTGAGCGACTGTGCCTGCGCTTGCCGATTGATCAATGCAGACAAAGCTGTTTCAACATCTTCCGAAGCCTCCAGCACAGACTGCCGATAGATGGCCAGGGCTTCGGCTTCACGTCCTTTGGCCTGGTCTATCTGCGCGTTGATGCGGCTAAAATCAAAGAGTCTCCACCTCAAACCGAG  >  NC_002947/3731884‑3732175

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