Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,389,951	C→G	11.4%	P25A (CCG→GCG)	PP_1212 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,389,951	0	C	G	11.4%	136.2 / 9.2	53	P25A (CCG→GCG)	PP_1212	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/26);  new base G (3/3);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGACT  >  NC_002947/1389811‑1390098                                                                                                                                             |                                                                                                                                                    ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                             <  5:277054/149‑1 (MQ=255) ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATc                                                                                                                                             >  6:107419/1‑149 (MQ=255) ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATc                                                                                                                                             >  5:275360/1‑149 (MQ=255)    aCCATATATAATCCCTCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTTCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAt                                                                                                                                          <  2:95345/149‑1 (MQ=255)    aCCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCg                                                                                                                                            <  5:51164/147‑1 (MQ=255)    aCCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGa                                                                                                                                           <  7:102527/148‑1 (MQ=255)           ataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGAt                                                                                                                                              >  1:153264/1‑138 (MQ=255)           ataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGAt                                                                                                                                              <  2:153264/138‑1 (MQ=255)               tCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTg                                                                                                                               >  3:278830/1‑149 (MQ=255)                 ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                              >  8:209687/1‑148 (MQ=255)                         tgctgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGcc                                                                                                                      <  6:286266/148‑1 (MQ=255)                            tgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCa                                                                                                                  >  2:327374/1‑149 (MQ=255)                                  ccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGaaa                                                                                                            <  7:318756/149‑1 (MQ=255)                                       agccagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgct                                                                                                       <  6:275360/149‑1 (MQ=255)                                              gTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTc                                                                                                 <  7:209687/148‑1 (MQ=255)                                                  tGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                                            <  1:327374/149‑1 (MQ=255)                                                  tGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                                            <  4:278830/149‑1 (MQ=255)                                                    gTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                                            >  6:128647/1‑147 (MQ=255)                                                     ttGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAg                                                                                         >  5:60678/1‑149 (MQ=255)                                                      tGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                                        >  4:285140/1‑149 (MQ=255)                                                        aCGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAg                                                                                         <  6:276730/146‑1 (MQ=255)                                                        aCGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAg                                                                                         >  5:276730/1‑146 (MQ=255)                                                           agagaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                                   <  5:107419/149‑1 (MQ=255)                                                                  cccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAg                                                                            <  6:67116/149‑1 (MQ=255)                                                                   ccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAtt                                                                                                                                         <  4:181591/87‑1 (MQ=255)                                                                   ccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAtt                                                                                                                                         >  3:181591/1‑87 (MQ=255)                                                                         cccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGATGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAAc                                                                     >  8:22969/1‑149 (MQ=255)                                                                          ccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCa                                                                                                                            >  5:142132/1‑93 (MQ=255)                                                                          ccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCa                                                                                                                            <  6:142132/93‑1 (MQ=255)                                                                               aTGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTg                                                                <  7:121508/148‑1 (MQ=255)                                                                               aTGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGAAGAAGATCAGCGCCGAGACGCAGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAACCTGCAGTCGGTGCCAGCCTTTCGCCTCAGCGGTAACGGTTgg                                                               >  8:143003/1‑149 (MQ=255)                                                                                tGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTgg                                                               <  3:8636/148‑1 (MQ=255)                                                                                   cTATCAATGTGCGTCCTGCGAGAACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGACAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                            >  2:231585/1‑148 (MQ=255)                                                                                        aaTGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAcc                                                      >  5:178843/1‑149 (MQ=255)                                                                                         aTGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCAGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAcc                                                      >  1:265885/1‑148 (MQ=255)                                                                                            tgCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGAc                                                   <  7:22969/148‑1 (MQ=255)                                                                                              cGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                                <  6:60678/149‑1 (MQ=255)                                                                                                tCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCaa                                              <  6:178843/149‑1 (MQ=255)                                                                                                tCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCaa                                              <  3:285140/149‑1 (MQ=255)                                                                                                           ccGCATGGAAGTGCTGCAGAAGATCAGAGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGTGCg                                    >  5:339270/1‑148 (MQ=255)                                                                                                                  gAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGa                             <  7:143003/148‑1 (MQ=255)                                                                                                                     gTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCt                         <  2:265885/149‑1 (MQ=255)                                                                                                                      tgctgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTg                        <  6:280965/149‑1 (MQ=255)                                                                                                                           cAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGACAGGCTTTCGCCTGAGCTGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCa                     <  1:231585/147‑1 (MQ=255)                                                                                                                             gaagaTCAGCGCCGCGCGGCTGATCGATTGCCCGGCCTGTCAGGCGCGGGCACTGAAAAAGCTGCTGTCGGTGCCAGGGTTTCGCGTGAGCGGTAGCGGTGGGTACGAATCCGAGTTCTAGGCGGGGGGGGAAAAGAATCTGGGAggcg                 >  2:299950/1‑149 (MQ=255)                                                                                                                             gaagaTCAGAGCCGCGCCGCTGATCGATTGCCCGGCCAGTCACGCGCCGGCACTGAAAAAGCAGCTGTCGGTGACAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAACGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAgccg                 >  2:284625/1‑149 (MQ=255)                                                                                                                                gaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                                       >  7:300685/1‑124 (MQ=255)                                                                                                                                gaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                                       <  8:300685/124‑1 (MQ=255)                                                                                                                                gaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc               >  8:262900/1‑148 (MQ=255)                                                                                                                                gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                                            >  4:171005/1‑71 (MQ=255)                                                                                                                                gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                                            <  3:171005/71‑1 (MQ=255)                                                                                                                                  tCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                               >  6:235674/1‑130 (MQ=255)                                                                                                                                  tCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                               <  5:235674/130‑1 (MQ=255)                                                                                                                                         cgcgccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGcc      >  2:290270/1‑148 (MQ=255)                                                                                                                                             ccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGACt  >  5:58337/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                                    CTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGACT  >  NC_002947/1389811‑1390098

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