Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	464,322	G→C	9.4%	H259D (CAC→GAC) ‡	davB ←	lysine 2‑monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	464,322	0	G	C	9.4%	158.8 / 6.7	53	H259D (CAC→GAC) ‡	davB	lysine 2‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (21/27);  new base C (2/3);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGC  >  NC_002947/464177‑464464                                                                                                                                                   |                                                                                                                                               gcgcGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCtg                                                                                                                                             <  8:226785/149‑1 (MQ=255)         aaTGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGgtcgt                                                                                                                                     >  1:55917/1‑149 (MQ=255)          aTGCGCTTGACCCCGGTACGCGTTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCGAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGgtcgtc                                                                                                                                    <  6:45200/149‑1 (MQ=255)                  gACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGtt                                                                                                                              <  7:1184/147‑1 (MQ=255)                  gACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTgg                                                                                                                            >  3:234899/1‑149 (MQ=255)                   aCCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTgg                                                                                                                            <  4:178934/148‑1 (MQ=255)                     cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtc                                                                                                                                       <  8:11528/135‑1 (MQ=255)                     cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtc                                                                                                                                       >  7:11528/1‑135 (MQ=255)                                    gccgccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCgg                                                                                                              >  3:279470/1‑145 (MQ=255)                                        ccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTc                                                                                                      <  6:122736/149‑1 (MQ=255)                                        ccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcggg                                                                                                                                    <  4:345182/119‑3 (MQ=255)                                             cAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCa                                                                                                 <  2:168344/149‑1 (MQ=255)                                                   gCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGt                                                                                           >  4:177842/1‑149 (MQ=255)                                                     tCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATc                                                                                               <  4:158378/143‑1 (MQ=255)                                                     tCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATc                                                                                               >  3:158378/1‑143 (MQ=255)                                                       aGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTggg                                                                                       >  5:13400/1‑149 (MQ=255)                                                         gCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTgggg                                                                                      >  3:38150/1‑148 (MQ=255)                                                                cccTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTGTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTcc                                                                                >  1:86438/1‑148 (MQ=255)                                                                cccTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGa                                                                              <  2:55917/149‑1 (MQ=255)                                                                     tGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTcc                                                                          >  1:224965/1‑148 (MQ=255)                                                                     tGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTc                                                                           <  7:5902/147‑1 (MQ=255)                                                                       gCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTcc                                                                          <  2:86438/146‑1 (MQ=255)                                                                       gCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCaa                                                                       >  4:342501/1‑149 (MQ=255)                                                                         cAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAcc                                                                     >  2:5872/1‑149 (MQ=255)                                                                           aTGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAccgc                                                                   <  1:26811/149‑1 (MQ=255)                                                                            tGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATc                                                                                               >  8:25965/1‑120 (MQ=255)                                                                            tGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATc                                                                                               <  7:25965/120‑1 (MQ=255)                                                                            tGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAccgcc                                                                  <  2:339327/149‑1 (MQ=255)                                                                               aCGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGt                                                               <  3:177842/149‑1 (MQ=255)                                                                                           ggTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGAc                                                   >  6:249857/1‑149 (MQ=255)                                                                                            gTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGAcc                                                  >  8:48929/1‑149 (MQ=255)                                                                                            gTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGAcc                                                  >  4:117745/1‑149 (MQ=255)                                                                                               cGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTgg                                               <  3:117745/149‑1 (MQ=255)                                                                                                    cgcCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGaa                                          <  4:279470/149‑1 (MQ=255)                                                                                                    cgcCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGCTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGaa                                          <  3:342501/149‑1 (MQ=255)                                                                                                            gCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTcgc                                  >  6:304289/1‑149 (MQ=255)                                                                                                             ccT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTcgc                                  <  2:224965/148‑1 (MQ=255)                                                                                                              cT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCgg                                <  4:172803/149‑1 (MQ=255)                                                                                                                   tgGCACTTGTTCCACGCCCCCGACCACCAGGTGCGGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCt                            <  8:102876/149‑1 (MQ=255)                                                                                                                     gCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTTGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGTAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCTGTGCTgg                          <  6:13400/149‑1 (MQ=255)                                                                                                                       aCTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTgg                          >  8:14776/1‑147 (MQ=255)                                                                                                                       aCTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTgg                          >  3:12322/1‑147 (MQ=255)                                                                                                                            ttCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAc                                                                      <  2:321614/98‑1 (MQ=255)                                                                                                                            ttCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAc                                                                      >  1:321614/1‑98 (MQ=255)                                                                                                                            ttCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCa                   <  5:249857/149‑1 (MQ=255)                                                                                                                                 cGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTg              <  4:12322/149‑1 (MQ=255)                                                                                                                                      ccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCGCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGGTGGAAGGTCGGCTTGGCaa          >  4:209548/1‑148 (MQ=255)                                                                                                                                         accaccAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAaga      <  7:48929/149‑1 (MQ=255)                                                                                                                                         accaccAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCt                            >  2:195715/1‑127 (MQ=255)                                                                                                                                         accaccAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCt                            <  1:195715/127‑1 (MQ=255)                                                                                                                                             ccaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAgcgc  <  5:304289/149‑1 (MQ=255)                                                                                                                                             ccaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAgcg   <  6:84147/148‑1 (MQ=255)                                                                                                                                             ccaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAgcg   >  5:84147/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                               GCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCT‑TGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGC  >  NC_002947/464177‑464464

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