Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,682,977 T→C intergenic (‑112/+53) qseG ← / ← qseE two‑component system QseEF‑associated lipoprotein QseG/two component system sensor histidine kinase QseE/GlrK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,682,9770TC96.2% 169.9 / ‑0.5 52intergenic (‑112/+53)qseG/qseEtwo‑component system QseEF‑associated lipoprotein QseG/two component system sensor histidine kinase QseE/GlrK
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (27/23);  total (28/24)
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.71e-01

ATGTGTCGCATATTCACCAGACTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACTGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGT  >  NZ_CP009273/2682855‑2683105
                                                                                                                          |                                                                                                                                
aTGTGTCGCATATTCACCAGACTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTg                                                                                                                  <  1:624355/139‑1 (MQ=255)
         atatTCACCAGACTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACaaataaa                                                                                                         <  1:152609/139‑1 (MQ=255)
              cACCAGACTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAg                                                                                                     <  2:306094/138‑1 (MQ=255)
               aCCAGACTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGcc                                                                                                   >  2:357590/1‑139 (MQ=255)
                 cAGACTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTg                                                                                                 <  1:610563/139‑1 (MQ=255)
                    aCTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGaa                                                                                                                     >  2:400135/1‑116 (MQ=255)
                    aCTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGaa                                                                                                                     <  1:400135/116‑1 (MQ=255)
                    aCTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTga                                                                                                >  2:504628/1‑137 (MQ=255)
                    aCTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTga                                                                                                >  2:623681/1‑137 (MQ=255)
                        aaaGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTgatagga                                                                                           <  1:790824/138‑1 (MQ=255)
                         aaGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTgataggata                                                                                         <  2:276669/139‑1 (MQ=255)
                                tCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGt                                                                                   >  1:608855/1‑138 (MQ=255)
                                      tGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCat                                                                             >  2:254947/1‑138 (MQ=255)
                                          cGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGcc                                                                                                                               <  2:495864/84‑1 (MQ=255)
                                          cGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGcc                                                                                                                               >  1:495864/1‑84 (MQ=255)
                                          cGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTAttt                                                                        >  2:192863/1‑139 (MQ=255)
                                                       gACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTcgac                                                           <  1:224249/139‑1 (MQ=255)
                                                       gACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTcgac                                                           <  2:839319/139‑1 (MQ=255)
                                                            tCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACTGCCCGGGGCCGAACTGACAAATAAAATAGCCTg                                                                                                 <  1:271125/96‑1 (MQ=255)
                                                            tCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACTGCCCGGGGCCGAACTGACAAATAAAATAGCCTg                                                                                                 >  2:271125/1‑96 (MQ=255)
                                                            tCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACaaataaaat                                                                                                       <  1:137480/90‑1 (MQ=255)
                                                            tCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACaaataaaat                                                                                                       >  2:137480/1‑90 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGaa                                                                                                                     >  1:660303/1‑66 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGaa                                                                                                                     <  2:660303/66‑1 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACaa                                                                                                              >  1:809575/1‑73 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACaa                                                                                                              <  2:809575/73‑1 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTga                                                                                                <  1:602695/87‑1 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTga                                                                                                >  2:602695/1‑87 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCCTTATTTCGTGTTTTTCGACGACGGTAATTCAATg                                            >  2:343/1‑139 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATg                                            >  2:689707/1‑139 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATg                                            >  2:756927/1‑139 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATg                                            >  1:620934/1‑139 (MQ=255)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACAATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGaa                                                                                                                     >  2:691671/1‑65 (MQ=17)
                                                                      gacaGCGCGGAGAAACCGGAGCGTACAATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGaa                                                                                                                     <  1:691671/65‑1 (MQ=17)
                                                                          gcgcGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCg                                                                                                                       <  2:258874/60‑1 (MQ=255)
                                                                          gcgcGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCg                                                                                                                       >  1:258874/1‑60 (MQ=255)
                                                                           cgcgGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGcc                                                                                                                               <  2:766669/51‑1 (MQ=255)
                                                                           cgcgGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGcc                                                                                                                               >  1:766669/1‑51 (MQ=255)
                                                                           cgcgGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTcgac                                                           >  2:92257/1‑119 (MQ=255)
                                                                           cgcgGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTcgac                                                           >  2:237381/1‑119 (MQ=255)
                                                                           cgcgGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTcgac                                                           <  1:92257/119‑1 (MQ=255)
                                                                           cgcgGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTcgac                                                           <  1:237381/119‑1 (MQ=255)
                                                                              ggAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGaaaca                                    <  2:608855/139‑1 (MQ=255)
                                                                                          gCGTACATATCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGccc                        >  2:406580/1‑139 (MQ=255)
                                                                                                   tCGTACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCtcgt                <  1:406580/138‑1 (MQ=255)
                                                                                                      tACGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGAc            >  1:32698/1‑139 (MQ=255)
                                                                                                        cGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTgataggata                                                                                         >  1:764296/1‑60 (MQ=255)
                                                                                                        cGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTgataggata                                                                                         <  2:764296/60‑1 (MQ=255)
                                                                                                        cGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGAc            >  2:427981/1‑137 (MQ=255)
                                                                                                        cGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGAc            >  2:33136/1‑137 (MQ=255)
                                                                                                        cGTGAAGATTTCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGAc            >  1:557289/1‑137 (MQ=255)
                                                                                                                tttCGAGCACCGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGt  <  1:858568/139‑1 (MQ=255)
                                                                                                                          |                                                                                                                                
ATGTGTCGCATATTCACCAGACTTAAAGCCTATCCCAGTGGGCGTAATTGTTGCAGACAGTCTGGACATGGACAGCGCGGAGAAACCGGAGCGTACATATCGTACGTGAAGATTTCGAGCACTGCCCGGGGCCGAACTGACAAATAAAATAGCCTGATAGGATAGGCTGTTAGCATTATTTCGTGTTTTTCGACGACGGTAATTCAATGCGGAAACAAACGTCTTGCCCGCTCTCGTCGACCAGATACAGT  >  NZ_CP009273/2682855‑2683105

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: