Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,797,772	G→C	R7G (CGC→GGC)	phoR ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,797,772	0	G	C	93.7%	212.1 / 6.5	63	R7G (CGC→GGC)	phoR	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base C (31/28);  total (33/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.00e-01

AAATATATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATCAACAATTCT  >  CP000730/1797637‑1797909                                                                                                                                        |                                                                                                                                          aaataTATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCcgtggt                                                                                                                                      <  2:26584/141‑1 (MQ=255)  aataTATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCgtggtg                                                                                                                                     >  2:125218/1‑141 (MQ=255)      tatGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAAc                                                                                                                                 <  1:424773/141‑1 (MQ=255)      tatGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAAc                                                                                                                                 <  2:403966/141‑1 (MQ=255)       atGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACt                                                                                                                                >  2:168262/1‑141 (MQ=255)        tGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCCGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACtt                                                                                                                               >  1:307263/1‑141 (MQ=255)            cTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTtcatc                                                                                                                           >  2:416848/1‑141 (MQ=255)               gCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCAtt                                                                                                                        <  1:293962/141‑1 (MQ=255)                  tCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGtt                                                                                                                     >  1:541237/1‑141 (MQ=255)                    tGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCt                                                                                                                   >  2:403348/1‑141 (MQ=255)                     gTTGCTAGAGCTTGTTCGGGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCtt                                                                                                                  <  2:592528/141‑1 (MQ=255)                      ttGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                                                                                 >  2:50066/1‑141 (MQ=255)                      ttGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                                                                                 >  1:50500/1‑141 (MQ=255)                      ttGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                                                                                 <  2:178117/141‑1 (MQ=255)                       tGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTACTTTa                                                                                                                <  2:529338/141‑1 (MQ=255)                            gagCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCt                                                                                                           >  1:544955/1‑141 (MQ=255)                              gCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCttt                                                                                                         >  1:415035/1‑141 (MQ=255)                               cTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTc                                                                                                        >  2:269601/1‑141 (MQ=255)                               cTTGATCGTGTTTTTCTGATAAAACATATGTGATGGAAGCCTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCTATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTc                                                                                                        <  2:701278/141‑1 (MQ=255)                                 tGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCa                                                                                                      <  1:20728/141‑1 (MQ=255)                                 tGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCa                                                                                                      <  1:272218/141‑1 (MQ=255)                                   ttCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAt                                                                                                    >  2:490046/1‑141 (MQ=255)                                              cTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTaa                                                                                         >  1:409360/1‑141 (MQ=255)                                               tGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAg                                                                                        >  1:632504/1‑141 (MQ=255)                                                     aaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGc                                                                                  >  1:349633/1‑141 (MQ=255)                                                         tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAg                                                                              <  1:204053/141‑1 (MQ=255)                                                         tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAg                                                                              <  1:30266/141‑1 (MQ=255)                                                         tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACCTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAg                                                                              <  1:318677/141‑1 (MQ=255)                                                          atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGt                                                                             <  2:575033/141‑1 (MQ=255)                                                           tgtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGtt                                                                            >  2:79349/1‑141 (MQ=255)                                                            gtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTGCTCTATTGCTTATCAATAACATTAAGCCGTGGTGAAAATTCATAATTGTTATTTAGGTCTTTTAAATTTATAGCATAAGCAACGAACAGttt                                                                           >  1:339413/1‑141 (MQ=255)                                                                 ggAAGCGTGAATAATAATGCCTAATGTAACAAAACTTATAATTAATATACTTCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTtaata                                                                      <  2:339413/141‑1 (MQ=255)                                                                    aGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGt                                                                   >  2:50922/1‑141 (MQ=255)                                                                       gTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtgt                                                                <  1:167369/141‑1 (MQ=255)                                                                            taataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGctt                                                           <  1:279110/141‑1 (MQ=255)                                                                              ataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGcttct                                                         <  1:623112/141‑1 (MQ=255)                                                                                aatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTa                                                       >  1:656974/1‑141 (MQ=255)                                                                                aatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTa                                                       >  2:147616/1‑141 (MQ=255)                                                                                  tGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAAAATTAAGCCGTGGTGAGACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg                                                     >  2:477023/1‑141 (MQ=255)                                                                                     cTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGAtt                                                  <  2:94829/141‑1 (MQ=255)                                                                                      tAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTa                                                 <  1:514462/141‑1 (MQ=255)                                                                                        aTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTAtc                                               >  1:696620/1‑141 (MQ=255)                                                                                        aTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTAtc                                               >  1:696653/1‑141 (MQ=255)                                                                                          gTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTAtctt                                             <  1:256364/141‑1 (MQ=255)                                                                                          gTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTAtctt                                             <  2:544955/141‑1 (MQ=255)                                                                                                     gATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCttt                                                                                                         <  1:349324/70‑1 (MQ=255)                                                                                                     gATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCttt                                                                                                         >  2:349324/1‑70 (MQ=255)                                                                                                      aTAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTCGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTGAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTAtctc                                 <  1:608411/141‑1 (MQ=255)                                                                                                       taattaatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATctct                                <  1:269601/141‑1 (MQ=255)                                                                                                         attaatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTa                              <  1:79349/141‑1 (MQ=255)                                                                                                           taatATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTaaa                            >  2:458080/1‑141 (MQ=255)                                                                                                                  cTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGa                     >  2:413848/1‑141 (MQ=255)                                                                                                                    gCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGata                   <  1:458080/141‑1 (MQ=255)                                                                                                                          cAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACa             >  2:485517/1‑141 (MQ=255)                                                                                                                          cAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACa             >  1:22493/1‑141 (MQ=255)                                                                                                                          cAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATATTCTAATTTATCTCTTAAATGACTGATATGTACa             >  1:554990/1‑141 (MQ=255)                                                                                                                            aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATctct                                <  1:402217/120‑1 (MQ=255)                                                                                                                            aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATctct                                >  2:402217/1‑120 (MQ=255)                                                                                                                            aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGAATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATc           >  2:672103/1‑141 (MQ=255)                                                                                                                               aCATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATcaac        <  1:469289/141‑1 (MQ=255)                                                                                                                               aCATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATcaac        <  1:403348/141‑1 (MQ=255)                                                                                                                                cATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATcaaca       <  1:605553/141‑1 (MQ=255)                                                                                                                                    aaGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATCAACAATTc   >  1:222207/1‑141 (MQ=255)                                                                                                                                     aGCCGTGGTGAAACTTCATCATTGCTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACGGTTTTAATAAGTTGTGGCTTCTTAGGTTTATCTTCTAATTTATCTCTTAAATGCCTGATATGAACATCAACAATTCt  >  1:518739/1‑141 (MQ=255)                                                                                                                                        |                                                                                                                                          AAATATATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATCAACAATTCT  >  CP000730/1797637‑1797909

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