Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	147,479	A→G	14.3%	L204P (CTG→CCG)	PP_0139 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	147,479	0	A	G	14.3%	98.5 / 7.5	42	L204P (CTG→CCG)	PP_0139	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base G (3/3);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAACCACCTTCTACGTACTGGCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTGGCGGATTTCGAAGCGCTTGCTG  >  NC_002947/147332‑147616                                                                                                                                                      |                                                                                                                                          gAAACCACCTTCTACGTACTGGCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAg                                                                                                                                          <  1:81300/149‑1 (MQ=255) gAAACCACCTTCTACGTACTGGCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAg                                                                                                                                          <  3:125708/149‑1 (MQ=255)      caccTTCTACGTACTGGCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCGGTtgtg                                                                                                                                     >  7:91034/1‑149 (MQ=255)      caccTTCTACGTACTGGCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCGGTtgtg                                                                                                                                     >  7:95887/1‑149 (MQ=255)             tACGTACTGGCCGTGTACTTTGGTGCCGTGGGCAGTCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGAt                                                                                                                               >  1:47752/1‑149 (MQ=255)               cGTACTGGCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCTCTCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGa                                                                                                                                <  1:72604/146‑1 (MQ=255)                    tGGCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTgg                                                                                                                       >  4:272180/1‑149 (MQ=255)                      gCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTgg                                                                                                                       <  1:299691/147‑1 (MQ=255)                      gCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGt                                                                                                                     <  4:9688/149‑1 (MQ=255)                        cGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTCCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGa                                                                                                                   >  4:67990/1‑149 (MQ=255)                         gtgtACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTGTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGAt                                                                                                                  >  5:307467/1‑149 (MQ=255)                          tgtACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCGGTTGTGCCTGATTGGCCTGGGTGATg                                                                                                                 <  8:95887/149‑1 (MQ=255)                          tgtACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCGGTTGTGCCTGATTGGCCTGGGTGATg                                                                                                                 <  8:91034/149‑1 (MQ=255)                          tgtACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATg                                                                                                                 <  2:47752/149‑1 (MQ=255)                             aCTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGtt                                                                                                                                        <  5:249341/123‑1 (MQ=255)                             aCTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGtt                                                                                                                                        >  6:249341/1‑123 (MQ=255)                              cTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTc                                                                                                              >  6:251033/1‑148 (MQ=255)                                   gTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGtt                                                                                                         <  3:272180/148‑1 (MQ=255)                                                 ccAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCGGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCcca                                                                                           >  5:309018/1‑149 (MQ=255)                                                        gTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTg                                                                                    <  6:307467/149‑1 (MQ=255)                                                                        gtTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTc                                                                                                              >  2:82054/1‑107 (MQ=255)                                                                        gtTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTc                                                                                                              <  1:82054/107‑1 (MQ=255)                                                                               aCTGCTGGCCGAGTTTTCAGGTGTAGTCGAGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCgg                                                             <  6:296718/149‑1 (MQ=255)                                                                               aCTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCA‑TCGTCTGCTACTGGTTCTTTGACTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCgg                                                             >  5:296718/1‑149 (MQ=255)                                                                                   ctgGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCCGGCTTCGGCCAGGCCAGCGAc                                                         <  3:67990/149‑1 (MQ=255)                                                                                                 aGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGctggcctg                                                                        >  3:136483/1‑120 (MQ=255)                                                                                                 aGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGctggcctg                                                                        <  4:136483/120‑1 (MQ=255)                                                                                                 aGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCActgct                                                  >  5:315434/1‑142 (MQ=255)                                                                                                 aGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCActgct                                                  <  6:315434/142‑1 (MQ=255)                                                                                                            cggcgATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTg                                      >  8:297793/1‑143 (MQ=255)                                                                                                            cggcgATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTg                                      <  7:297793/143‑1 (MQ=255)                                                                                                             ggcgATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCa                               <  5:251033/149‑1 (MQ=255)                                                                                                                gATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCg                            >  3:54360/1‑149 (MQ=255)                                                                                                                             cTACTGGTTCTTTGCCTGAAGCTTCGGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTGGCGGAttt                <  6:309018/149‑1 (MQ=255)                                                                                                                                cTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCGGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTGGCGGATTTCg              >  5:300117/1‑148 (MQ=255)                                                                                                                                  ggTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTgtcgtc                                 <  6:330966/127‑1 (MQ=255)                                                                                                                                  ggTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTgtcgtc                                 >  5:330966/1‑127 (MQ=255)                                                                                                                                       tttGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTg                        >  7:290035/1‑131 (MQ=255)                                                                                                                                       tttGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTg                        <  8:290035/131‑1 (MQ=255)                                                                                                                                           ccTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTGGCGGATTTCGAAGCGCTTGCTg  >  1:24778/1‑149 (MQ=255)                                                                                                                                           ccTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGATGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTGGCGGATTTCGAAGCGCTTGCTg  >  3:15498/1‑149 (MQ=255)                                                                                                                                           ccTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGAGATGGTCCAGTTCACGACTTGCCCCACCAGCTGATCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTGGCGGATTTCGAAGCGCTTGCTg  >  1:109078/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                          GAAACCACCTTCTACGTACTGGCCGTGTACTTTGGTGCCGTGGGCA‑TCCAGCGCGTGCGCCATGCGGTCGGTTGCGCACTGCTGGCCGAGTTTTCAGGGGTAGTCGCGGCGATCT‑TCGTCTGCTACTGGTTCTTTGCCTGAAGCTTCAGTTGTGCCTGATTGGCCTGGGTGATGGTCCAGTTCACGACTTGCCCCACCAGCTGGTCGCTGGCCTGGCCGAAGCCGGCGACCACTGCTGATACCTGTTTGTCGTCCAGCGGTTGGCGGATTTCGAAGCGCTTGCTG  >  NC_002947/147332‑147616

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