Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,382,535	G→T	7.9%	L256I (CTC→ATC)	PP_1204 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,382,535	0	G	T	7.9%	150.5 / 4.2	51	L256I (CTC→ATC)	PP_1204	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base T (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

CGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCT  >  NC_002947/1382389‑1382675                                                                                                                                                    |                                                                                                                                             cGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGt                                                                                                                                            >  8:84621/1‑149 (MQ=255)   aCTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCa                                                                                                                                          >  1:131404/1‑149 (MQ=255)       cAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATc                                                                                                                                        <  4:243592/147‑1 (MQ=255)       cAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                       >  7:208336/1‑148 (MQ=255)        aGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCTGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                       >  1:189319/1‑147 (MQ=255)          gAACTCCTGCGGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                       >  2:25630/1‑145 (MQ=255)            aCTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGATGACCTCGAAGAtg                                                                                                                                 <  2:32048/149‑2 (MQ=255)            aCTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATc                                                                                                                                 <  3:129131/149‑1 (MQ=255)            aCTCCTGATGCAGCTCCGGATCGTTGCGGGTCAGCTCGATGAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGATGACCTCGAAGAtg                                                                                                                                 <  8:162923/149‑2 (MQ=255)                   ctgcAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCAt                                                                                                                          <  8:208336/149‑1 (MQ=255)                         cTCCGGATCGTTGCGGGTCAGCGCGATCAGGCTCTGTTCCATTTCACTGGCTTACTCTTCAAGGCCCAGTTCGGA‑CAGGCGCGTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCtt                                                                                                                    <  7:250468/149‑1 (MQ=255)                             ggATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGag                                                                                                                <  2:131404/149‑1 (MQ=255)                                  gTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGagcagc                                                                                                            <  1:45281/148‑1 (MQ=255)                                    tGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTa                                                                                                         >  4:207675/1‑149 (MQ=255)                                        gggcAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTgc                                                                                                     <  8:222121/146‑1 (MQ=255)                                               tCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGt                                                                                              >  5:97226/1‑149 (MQ=255)                                                             ttCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTACACCCACTGGCAGACATCATCATCTTCGAGGTCATAGAAGATCAGCTCATGGGATTAGAGCAGCTTATTGCGCAGGGTCGCTCTTACCaa                                                                                  >  6:220355/1‑147 (MQ=255)                                                                 aTTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCt                                                                            >  4:224139/1‑149 (MQ=255)                                                                 aTTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCt                                                                            <  1:25630/149‑1 (MQ=255)                                                                  tttCACTGGCTTCCTCTTCAAGGCCCAGTTCGGAGCAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCt                                                                            >  1:36519/1‑149 (MQ=255)                                                                  tttCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                           <  7:84621/149‑1 (MQ=255)                                                                   ttCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTgg                                                                          >  4:146653/1‑149 (MQ=255)                                                                      aCTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGagcag                                                                                                             <  2:236991/111‑1 (MQ=255)                                                                      aCTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGagcag                                                                                                             >  1:236991/1‑111 (MQ=255)                                                                      aCTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCAt                                                                       >  2:268773/1‑149 (MQ=255)                                                                            ttCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGGTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                           >  1:219834/1‑139 (MQ=255)                                                                            ttCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                           <  2:219834/139‑1 (MQ=255)                                                                               ctctTCAAGGCCCAGTTCGGA‑CATGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGt                                                               <  7:191718/148‑1 (MQ=255)                                                                                    cAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATAAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCa                                                         >  5:260591/1‑149 (MQ=255)                                                                                             gTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCg                                                >  5:115483/1‑149 (MQ=255)                                                                                                 ggA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCaa                                            <  6:115483/149‑1 (MQ=255)                                                                                                  gA‑CAGGCGCTTGACCCGGCGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGagca                                                                                                              <  6:167454/82‑1 (MQ=255)                                                                                                  gA‑CAGGCGCTTGACCCGGCGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGagca                                                                                                              >  5:167454/1‑82 (MQ=255)                                                                                                   a‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGc                                                                    >  7:92691/1‑123 (MQ=255)                                                                                                   a‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGc                                                                    <  8:92691/123‑1 (MQ=255)                                                                                                     cAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTa                                         <  2:36519/149‑1 (MQ=255)                                                                                                        gcgcTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACtt                                      >  5:114669/1‑149 (MQ=255)                                                                                                                 ccGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtgca                             <  2:189319/149‑1 (MQ=255)                                                                                                                  ctgTGAACCCACTGGCCGACATCATAAACTACGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGCTGTCGCTGACAAACAGGCTGGCATCGCCTCGTACCGCATTGACCTCGTCAAGTAATtcaaggtcaac                            >  6:153916/3‑145 (MQ=255)                                                                                                                       cACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGc                       <  3:207675/149‑1 (MQ=255)                                                                                                                          ccACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                           >  2:103891/1‑94 (MQ=255)                                                                                                                          ccACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                           <  1:103891/94‑1 (MQ=255)                                                                                                                                gCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTgg               <  1:268773/148‑1 (MQ=255)                                                                                                                                   gACATCTTCATCTTCGAGGTCAACGAAGATCAGCTCATGTGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGa           <  5:220355/149‑1 (MQ=255)                                                                                                                                     catcatcatCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAgg         <  6:114669/149‑1 (MQ=255)                                                                                                                                         atcatcTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGtt             <  4:74233/141‑1 (MQ=255)                                                                                                                                         atcatcTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGtt             >  3:74233/1‑141 (MQ=255)                                                                                                                                         atcatcTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCt     >  7:173838/1‑149 (MQ=255)                                                                                                                                            atcTTCGATGACCTCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCt  >  4:151569/1‑149 (MQ=255)                                                                                                                                            atcTTCGATGACCTCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCt  >  2:53793/1‑149 (MQ=255)                                                                                                                                            atcTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCt  >  4:226957/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                             CGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGA‑CAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCT  >  NC_002947/1382389‑1382675

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