Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	77,117	G→A	R61R (CGC→CGT)	sgrR ←	DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	77,117	0	G	A	96.9%	98.0 / NA	32	R61R (CGC→CGT)	sgrR	DNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (9/22);  total (9/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 = 1.00e+00

ACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACCAGTTGATCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCGCGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGTCAAT  >  NC_000913/76877‑77356                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                aCGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCATCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACtt                                                                                                                                                                                                                                       <  1:92575/251‑1 (MQ=255)     cgcAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGc                                                                                                                                                                                                                                   >  2:149925/1‑251 (MQ=255)          tGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTccc                                                                                                                                                                                                                              >  2:116747/1‑251 (MQ=255)                  cAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACtt                                                                                                                                                                                                                                       >  2:114806/1‑234 (MQ=255)                  cAGGTAGCAGACTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACtt                                                                                                                                                                                                                                       <  1:114806/234‑1 (MQ=255)                              tACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTcc                                                                                                                                                                                                                                     >  1:112653/1‑224 (MQ=255)                              tACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTcc                                                                                                                                                                                                                                     <  2:112653/224‑1 (MQ=255)                                             agtagAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGc                                                                                                                                                                                           <  2:226833/251‑1 (MQ=255)                                                   gCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTGGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTg                                                                                                                                                                                     <  2:88115/251‑1 (MQ=255)                                                                         cTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGagcagc                                                                                                                                                               >  1:8123/1‑251 (MQ=255)                                                                                     gCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCt                                                                                                                                                    <  1:71594/250‑1 (MQ=255)                                                                                                               gCCGCACAGTCGCTTTGTCGCCAACCGACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGGCAGACGCGAGCGTTTGCCGCGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGGGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGc                                                                                                                         <  1:109110/251‑1 (MQ=255)                                                                                                                       gTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTAcc                                                                                                                 <  1:218911/251‑1 (MQ=255)                                                                                                                         cGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTAccc                                                                                                                <  2:59143/250‑1 (MQ=255)                                                                                                                                    caaccaacTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGc                                                                                                    <  1:180073/251‑1 (MQ=255)                                                                                                                                                   gttggtcgccaccctgcTGCTCCAGAAGGTCCCCCGCACGCTGTTGCTGCAGCGCCAGCCCGGTATAGCGGCATGCCAGACGCGCGCGTTTACCACGCCCGACTTCCGCTTCCCACCTCAGCCAGCCGCGATCCTGCATGGTGGTGAGCAGGGTGCGCATATGACGACGCGAGCAGATCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCATGCGATTTACCCTCGCGGCATTGACACCGGCGGCTGaa                                                                                      >  2:267188/17‑250 (MQ=255)                                                                                                                                                              cGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAg                                                                          <  1:166333/251‑1 (MQ=255)                                                                                                                                                                     gCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCGCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAg                                                                     <  2:40992/249‑1 (MQ=255)                                                                                                                                                                                                        gcCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACa                                <  1:102741/251‑1 (MQ=255)                                                                                                                                                                                                         ccAGCCCGGTATAGAGGAATGTCATACGCGAGCGTTGACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGAGGGCAGAAAAGGGGAACTCCTGTGCAAAAGACAg                               <  2:259096/251‑1 (MQ=255)                                                                                                                                                                                                          cAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGc                              <  1:144010/251‑1 (MQ=255)                                                                                                                                                                                                            gCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCaa                            <  1:43561/251‑1 (MQ=255)                                                                                                                                                                                                               cGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAattt                         <  1:215756/251‑1 (MQ=255)                                                                                                                                                                                                                   ataaagGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGGAATGATGTTGCAGACGAGCAGCTGGCATAAAAAGGGAACTCCTGTGCAAAAGACAGCAattttatt                     >  2:5168/5‑251 (MQ=255)                                                                                                                                                                                                                   ataGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGAt                                                                                                                      <  2:208944/154‑1 (MQ=255)                                                                                                                                                                                                                   ataGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGAt                                                                                                                      >  1:208944/1‑154 (MQ=255)                                                                                                                                                                                                                     agagGAATGTCAGTCGCGAGCGTTGCCCACGCCCGACTTCCGCCTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGTGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGTACAAGACAGCAattttatttt                   <  2:246261/251‑1 (MQ=255)                                                                                                                                                                                                                      gcggAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCGCGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCTGCGTGCGCATATGACGACGCGCGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTc                  <  2:100534/249‑1 (MQ=255)                                                                                                                                                                                                                         gAATGTCATACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCGGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCAGTGCCACAGGCGGAGGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAACGACAGCAATTTTATTTTCCCt               <  2:122284/251‑1 (MQ=255)                                                                                                                                                                                                                           aTGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCtat             <  2:8683/251‑1 (MQ=255)                                                                                                                                                                                                                                  aCGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGt      <  1:129259/251‑1 (MQ=255)                                                                                                                                                                                                                                      gAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGTCaat  >  1:257174/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                ACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACCAGTTGATCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCGCGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGAACTGTTGTTGCAGACGAGCAGATGGCATAAAAGGGGAACTCCTGTGCAAAAGACAGCAATTTTATTTTCCCTATATTAAGTCAAT  >  NC_000913/76877‑77356

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑