Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	544,918	G→A	14.3%	R369H (CGC→CAC)	ylbE →	NO‑induced DUF1116 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	544,918	0	G	A	14.3%	63.6 / 4.6	28	R369H (CGC→CAC)	ylbE	NO‑induced DUF1116 protein
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base A (2/2);  total (13/15)
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.51e-01

CCTGTAAACACTCCGCAAGGTCTGTTTTTCACCGGCTTCTCGCAGGAGCAGGCGAACCCGGATATGGGCGATAGCGCGATTACCGAAACCTTTGGTATCGGAGGTGCGGCAATGATCGCAGCGCCTGGCGTAACGCGCTTTGTCGGTGCGGGTGGCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGG  >  CP009273/544644‑545192                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                   ccTGTAAACACTCCGCAAGGTCTGTTTTTCACCGGCTTCTCGCAGGAGCAGGCGAACCCGGATATGGGCGATAGCGCGATTACCGAAACCTTTGGTATCGGAGGTGCGGCAATGATCGCAGCGCCTGGCGTAACGCGCTTTGTCGGTGCGGGTGGCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATtcgtcg                                                                                                                                                                                                                                                                                    >  1:420781/1‑275 (MQ=255)                cgcgttCTGTTTTTCACCGGCTTCGCGCCGTCGCAGGCTCACCCGTATAGGGGCGATATCGAGAGTCCCGTAACCTTTGGTATCGGAGGTGCGGCAATGATCGCAGCGCCTGGCGTAACGCGTGTTGTCGGTGCGGGTGGCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGGTGGACATTCGTCGCGTGGTAGAAACCgg                                                                                                                                                                                                                                                                     <  2:650170/270‑1 (MQ=255)                                                        aCCCGGATATGGGCGATAGCGCGATTACCGAAACCTTTGGTATCGGAGGTGCGGCAATGATCGCAGCGCCTGGCGTAACGCGCTTTGTCGGTGCGGGTGGCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTCTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAAACGGACTGAAGCCCCGCATCAATACCGGTATCGCCCATgaacta                                                                                                                                                                                                                             >  2:220849/1‑272 (MQ=255)                                                                                                                 tagtCGCGGGGCCTGGCGTAACGCGGTTTGGCGGTGCGGGGGGCATGGAAGCGGCAATGGCGGTATCTGAAGATATGGCGGAAATTTACCTTGAACGCAATATGCAGGTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCACGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACAGTGCGGGCACCGCTGGCGTGCTTTGaa                                                                                                                                                                    <  2:18037/272‑1 (MQ=255)                                                                                                                          cgcCTGGCGTAACGCGCTTTGTCGGTGCGGGTGGCATGGAAGCGGCAATAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCACGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACAGTGCGGGCACCGCTGGCGTGCTTTGAACAGgcgcg                                                                                                                                                            >  1:18037/1‑273 (MQ=255)                                                                                                                                 ggaaggggggtttGTTGGGGAGGTTGGCATGGGGCCGGCAAGAGGNGTATCTGAAGGGTTGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCa                                                                                                                                                     <  2:420776/264‑1 (MQ=255)                                                                                                                                          cTTTGTCGGTGCGGGTGGCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAg                                                                                                                                                                                                                                                                            >  2:419923/1‑146 (MQ=255)                                                                                                                                              gTCGGTGCGGGTGGCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAAc                                                                                                                                                                                                                                                                        <  1:419928/146‑1 (MQ=255)                                                                                                                                                  gggggggTGGCATGGAAGCGGCAAGATCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCATGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGt                                                                                                                                   <  1:220853/271‑1 (MQ=255)                                                                                                                                                       ggTGGCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGaaa                                                                                                                                                                                                                                                                         >  2:305450/1‑136 (MQ=255)                                                                                                                                                           gCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCgg                                                                                                                                                                                                                                                                     <  1:305455/136‑1 (MQ=255)                                                                                                                                                                                          gaTGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTAc                                                                                                                                                                                                                                                               >  1:80663/1‑111 (MQ=255)                                                                                                                                                                                            tGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCggg                                                                                                                                                                                        >  1:247937/1‑180 (MQ=255)                                                                                                                                                                                             ggCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCttt                                                                                                                                                                       >  2:151415/1‑196 (MQ=255)                                                                                                                                                                                              gCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCa                                                                                                                                                                                                                                                           <  2:80663/111‑1 (MQ=255)                                                                                                                                                                                                ggAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCAcc                                                                                                                                                                                    <  2:247933/180‑1 (MQ=255)                                                                                                                                                                                                 gAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAAc                                                                                                                                                                   <  1:151418/196‑1 (MQ=255)                                                                                                                                                                                                                                     ccAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGt                                                <  1:397955/275‑1 (MQ=255)                                                                                                                                                                                                                                      cAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGAt                                                                                                                                                                                                                     >  1:28843/1‑109 (MQ=255)                                                                                                                                                                                                                                          cTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCggg                                                                                                                                                                                                                 <  2:28843/109‑1 (MQ=255)                                                                                                                                                                                                                                              gcttttCAGGCCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCACTACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCGGAAAGCATGGGGATTGGTTGAGGat                                                                                                                      >  1:642737/3‑195 (MQ=255)                                                                                                                                                                                                                                              gATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTAt                                                                                                                                                                                                                                          >  1:13642/1‑80 (MQ=255)                                                                                                                                                                                                                                              gATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCa                                                                                                                                                                                                               >  2:308275/1‑107 (MQ=255)                                                                                                                                                                                                                                                  ttCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGAtt                                                                                                                                                                                                           <  1:308280/107‑1 (MQ=255)                                                                                                                                                                                                                                                  ttCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGGAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGcc                                                                                                                                                                                                                                      <  2:13642/80‑1 (MQ=255)                                                                                                                                                                                                                                                             tGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCa                                                                                                                                                                                                                                    >  2:562884/1‑71 (MQ=255)                                                                                                                                                                                                                                                                 tGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATaaa                                                                                                                                                                                                                                <  1:562895/71‑1 (MQ=255)                                                                                                                                                                                                                                                                         gACATTCGTCACGTGGTAGAAACCGGCATTACGCCACTCATCAATACCg                                                                                                                                                                                                                                              >  1:52779/1‑49 (MQ=255)                                                                                                                                                                                                                                                                             ttCGTCACGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTAt                                                                                                                                                                                                                                          <  2:52779/49‑1 (MQ=255)                                                                                                                                                                                                                                                                                   gcgTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTgg                                                                                                                                                                               >  1:653233/1‑102 (MQ=255)                                                                                                                                                                                                                                                                                   gcgTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAAACTGATGACGGAAAgcgg  >  1:102347/1‑275 (MQ=255)                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                   CCTGTAAACACTCCGCAAGGTCTGTTTTTCACCGGCTTCTCGCAGGAGCAGGCGAACCCGGATATGGGCGATAGCGCGATTACCGAAACCTTTGGTATCGGAGGTGCGGCAATGATCGCAGCGCCTGGCGTAACGCGCTTTGTCGGTGCGGGTGGCATGGAAGCGGCAAGAGCGGTATCTGAAGAGATGGCGGAAATTTACCTTGAACGCAATATGCAGTTGCAGATCCCAAGCTGGGATTTTCAGGGCGCGTGCCTGGGGCTGGACATTCGTCGCGTGGTAGAAACCGGCATTACGCCACTCATCAATACCGGTATCGCCCATAAAGAGGCGGGGATCGGGCAGATTGGCGCAGGCACCGTGCGGGCACCGCTGGCGTGCTTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCACCGAATTATCGCCAGTCCTGGCGGTTAGCGGGAGTGTGGCGGCGGGCGATTAACCTGATGACGGAAAGCGG  >  CP009273/544644‑545192

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑