Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,339,900	C→G	10.7%	R376R (CGC→CGG)	pcnB →	poly(A) polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,339,900	0	C	G	10.7%	157.7 / 10.6	56	R376R (CGC→CGG)	pcnB	poly(A) polymerase
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base G (3/3);  total (29/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 = 8.57e-01

GCAGAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGCGC  >  NC_002947/5339753‑5340047                                                                                                                                                    |                                                                                                                                                    gCAGAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCa                                                                                                                                                    >  6:343289/1‑149 (MQ=255)   aGAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGc                                                                                                                                                  <  6:104412/149‑1 (MQ=255)    gAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCg                                                                                                                                                 <  5:148430/149‑1 (MQ=255)        cAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCTTGCCGCTGCGGCGCGGCag                                                                                                                                             >  6:97273/1‑148 (MQ=255)          gggCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCaatc                                                                                                                                           <  2:316643/149‑3 (MQ=255)              gTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCTTGCCGCTGCGGCGCGGCAGgcgctc                                                                                                                                       <  5:97273/149‑3 (MQ=255)                   gccgATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACa                                                                                                                                  <  7:80712/149‑1 (MQ=255)                          ccGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATgctgct                                                                                                                           <  7:263490/149‑1 (MQ=255)                           cGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGctgctg                                                                                                                          >  1:312726/1‑149 (MQ=255)                                aTGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACaa                                                                                                                     >  2:128304/1‑149 (MQ=255)                                     cGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCcgc                                                                                                                <  5:364901/149‑1 (MQ=255)                                      ggTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCcgcg                                                                                                               >  5:121423/1‑149 (MQ=255)                                           ggCCCACGACCTGATCGCCGAGCAGTGCGCGCTCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGGGAGATCTGGGACATGCATGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTcc                                                                                                          <  3:283929/149‑1 (MQ=255)                                           ggCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTcc                                                                                                          <  1:105622/149‑1 (MQ=255)                                           ggCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTcc                                                                                                          >  3:182216/1‑149 (MQ=255)                                                     cTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTa                                                                                                <  8:174853/149‑1 (MQ=255)                                                       gATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGCGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACg                                                                                              >  2:148374/1‑149 (MQ=255)                                                       gATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGctgctg                                                                                                                          <  6:136287/121‑1 (MQ=255)                                                       gATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGctgctg                                                                                                                          >  5:136287/1‑121 (MQ=255)                                                            ccGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCTGCGGCGCGGCag                                                                                                                                             <  3:161353/97‑2 (MQ=255)                                                            ccGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCTGCGGCGCGGCag                                                                                                                                             >  4:161353/1‑96 (MQ=255)                                                            ccGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTc                                                                                         <  8:171906/149‑1 (MQ=255)                                                            ccGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTc                                                                                         <  5:343289/149‑1 (MQ=255)                                                              gAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCct                                                                                       <  2:312726/149‑1 (MQ=255)                                                               accAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTcc                                                                                        <  1:128304/145‑1 (MQ=255)                                                                    tGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCctgctgct                                                                                 >  5:274827/1‑149 (MQ=255)                                                                         gcgcATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTCGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGt                                                                             >  1:202156/1‑148 (MQ=255)                                                                            cATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGaaa                                                                         <  2:202156/149‑1 (MQ=255)                                                                                   aTCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCgg                                                                  >  6:112274/1‑149 (MQ=255)                                                                                                 tACCCTGCCGATCCGCGAGATCTGGTGCATGCATGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATTCTTCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCctgct                                                                                    <  7:290753/117‑1 (MQ=255)                                                                                                 tACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACcgac                                                     <  6:121423/148‑1 (MQ=255)                                                                                                 tACCCTGCCGATCCGCGAGATCTGCGGCATGCAGGAGCGACTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTcctcct                                                                                    >  8:290753/1‑117 (MQ=255)                                                                                                    ccTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTc                                                                                                           >  8:167667/1‑91 (MQ=255)                                                                                                    ccTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTc                                                                                                           <  7:167667/91‑1 (MQ=255)                                                                                                       gCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCAGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGATTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTgg                                              >  5:144632/1‑149 (MQ=255)                                                                                                          gATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGc                                            >  8:203171/1‑148 (MQ=255)                                                                                                              cgcgAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtg                                       >  4:108835/1‑149 (MQ=255)                                                                                                               gcgAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtgg                                      >  6:282324/1‑149 (MQ=255)                                                                                                                cgAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtggt                                     >  5:344164/1‑149 (MQ=255)                                                                                                                cgAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtggt                                     >  4:340562/1‑149 (MQ=255)                                                                                                                 gagaTCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtggtg                                    >  7:240891/1‑149 (MQ=255)                                                                                                                  agaTCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtggtg                                    >  3:329592/1‑148 (MQ=255)                                                                                                                     tCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGAcc                                >  2:79164/1‑149 (MQ=255)                                                                                                                       tGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGa                              >  4:108583/1‑149 (MQ=255)                                                                                                                          gACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTa                           <  6:274827/149‑1 (MQ=255)                                                                                                                                 aGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTgg                                              <  4:265854/123‑1 (MQ=255)                                                                                                                                 aGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTgg                                              >  3:265854/1‑123 (MQ=255)                                                                                                                                 aGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGAc                    <  5:112274/149‑1 (MQ=255)                                                                                                                                   gcgcgcTTGCCGCTGCGGCGCGGCAGGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGga                                                               <  3:207014/102‑1 (MQ=37)                                                                                                                                   gcgcgcTTGCCGCTGCGGCGCGGCAGGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGga                                                               >  4:207014/3‑104 (MQ=37)                                                                                                                                        ccTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCgg                                                                  <  7:113605/96‑1 (MQ=255)                                                                                                                                        ccTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCgg                                                                  >  8:113605/1‑96 (MQ=255)                                                                                                                                        ccTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAAc              >  2:230051/1‑148 (MQ=255)                                                                                                                                        ccTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACg             <  4:329592/149‑1 (MQ=255)                                                                                                                                         cTGCCGCGCCGCAGCGGCAAGCGCGCCTACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTTCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGa            <  4:182216/149‑1 (MQ=255)                                                                                                                                         cTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGa            >  7:45582/1‑149 (MQ=255)                                                                                                                                            ccgcgccgcAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACa          <  7:203171/148‑1 (MQ=255)                                                                                                                                                   gcAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAgcgc  <  1:79164/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCAGAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGATCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGCGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGCGC  >  NC_002947/5339753‑5340047

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