Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,389,951	C→G	7.8%	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	7.8%	160.9 / 3.7	51	P25A (CCG→GCG)	PP_1212	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base G (2/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGA  >  NC_002947/1389807‑1390096                                                                                                                                                  |                                                                                                                                                  cttgctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCgccgcgccgc                                                                                                                                                <  4:118445/148‑1 (MQ=255)     ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATc                                                                                                                                           >  2:291402/1‑149 (MQ=255)         ccATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAtt                                                                                                                                       >  4:357942/1‑149 (MQ=255)          cATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTg                                                                                                                                      >  7:47071/1‑149 (MQ=255)              tataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCg                                                                                                                                  >  4:182169/1‑149 (MQ=255)               ataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATcttctg                                                                                                                                      >  1:145217/1‑139 (MQ=255)               ataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATcttctg                                                                                                                                      <  2:145217/144‑6 (MQ=255)                  aTCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCt                                                                                                                              <  8:209476/149‑1 (MQ=255)                     ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCTCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGa                                                                                                                                         <  6:10537/135‑1 (MQ=255)                     ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGa                                                                                                                                         >  5:10537/1‑135 (MQ=255)                     ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                            <  1:291402/148‑1 (MQ=255)                         tttttGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAt                                                                                                                                        <  7:342554/132‑1 (MQ=255)                         tttttGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAt                                                                                                                                        >  8:342554/1‑132 (MQ=255)                               ctgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCgg                                                                                                                  <  5:286992/148‑1 (MQ=255)                                    ccccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTg                                                                                                             >  7:178416/1‑148 (MQ=255)                                    ccccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTg                                                                                                             >  1:160930/1‑148 (MQ=255)                                                 cGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTc                                                                                               <  4:175528/149‑1 (MQ=255)                                                   ttATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCgg                                                                                             >  8:243699/1‑149 (MQ=255)                                                   ttATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCAAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCgg                                                                                             >  8:270019/1‑149 (MQ=255)                                                                gagaCACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTg                                                                                                             >  3:267408/1‑120 (MQ=255)                                                                gagaCACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTg                                                                                                             <  4:267408/120‑1 (MQ=255)                                                                 agaCACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGc                                                                               <  8:47071/149‑1 (MQ=255)                                                                       ccATGCCCCTTTATGACTATCAATGTGTCGTCCTGCGAGCCCCGCATGGAAGTGCTGCAGAAAAACAGCGCCGCGCCGCTGAACAATTGCCCGGCCTGTCAGGCGCCGGCAATGAAAAAGCTGATGTCGTTGCCAGGCTTTCGCCTGAg                                                                          >  1:15894/1‑149 (MQ=255)                                                                                tttATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                                 >  3:59720/1‑132 (MQ=255)                                                                                tttATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                                 <  4:59720/132‑1 (MQ=255)                                                                                       cTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                           >  1:65506/1‑147 (MQ=255)                                                                                       cTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                           <  2:65506/147‑1 (MQ=255)                                                                                        tATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCa                                                                                                                          <  6:1280/83‑1 (MQ=255)                                                                                        tATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCa                                                                                                                          >  5:1280/1‑83 (MQ=255)                                                                                          tCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTg                                                                                                             >  5:61135/1‑94 (MQ=255)                                                                                          tCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTg                                                                                                             <  6:61135/94‑1 (MQ=255)                                                                                             aTGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAcc                                                    >  5:214036/1‑148 (MQ=255)                                                                                                   cGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                              <  8:178416/149‑1 (MQ=255)                                                                                                       cTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGcc                                                                              <  4:125850/113‑1 (MQ=255)                                                                                                       cTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGcc                                                                              >  3:125850/1‑113 (MQ=255)                                                                                                          cGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                                          <  4:88767/98‑1 (MQ=255)                                                                                                           gAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCgg                                      <  6:214036/149‑1 (MQ=255)                                                                                                            aGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                                     <  2:160930/149‑1 (MQ=255)                                                                                                            aGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                                     <  8:342740/149‑1 (MQ=255)                                                                                                             gCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCgggg                                    <  3:182169/149‑1 (MQ=255)                                                                                                                        aaGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAg                                           >  5:24890/1‑131 (MQ=255)                                                                                                                        aaGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAg                                           <  6:24890/131‑1 (MQ=255)                                                                                                                         aGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATc                        <  7:243699/149‑1 (MQ=255)                                                                                                                              tgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCg                                       >  6:7249/1‑129 (MQ=255)                                                                                                                               gcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAg                  >  3:117101/1‑149 (MQ=255)                                                                                                                                cAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAgg                 >  5:127852/1‑149 (MQ=255)                                                                                                                                 agaagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggc                >  6:231065/1‑149 (MQ=255)                                                                                                                                 agaagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggc                >  6:22018/1‑149 (MQ=255)                                                                                                                                  gaagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCTCTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAGAAGAATCTGGCAggcg               >  7:267997/1‑149 (MQ=255)                                                                                                                                   aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcgg              >  8:336850/1‑149 (MQ=255)                                                                                                                                      aTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                                     <  6:112875/123‑1 (MQ=255)                                                                                                                                      aTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                                     >  5:112875/1‑123 (MQ=255)                                                                                                                                          gcgccgcgccGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAg       >  7:136375/1‑149 (MQ=255)                                                                                                                                             ccgcgccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGcc    <  7:336850/149‑1 (MQ=255)                                                                                                                                               gcgccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGa  <  3:357942/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTG‑CGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGA  >  NC_002947/1389807‑1390096

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