Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,661,508 T→G T250P (ACC→CCC)  hcaR ← DNA‑binding transcriptional regulator HcaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,661,5080TG92.9% 39.3 / ‑5.1 14T250P (ACC→CCC) hcaRDNA‑binding transcriptional regulator HcaR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (5/8);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGC  >  NZ_CP009273/2661427‑2661592
                                                                                 |                                                                                    
tCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAat                                                                              <  1:230126/90‑1 (MQ=255)
      ccccTTTCTTCCACGCCATCAGGAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATGATTCa                                                                        >  1:402802/1‑90 (MQ=255)
            tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAAc                                                                  <  1:444636/90‑1 (MQ=255)
            tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAAc                                                                  <  2:162864/90‑1 (MQ=255)
                     ccATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCa                                                         >  1:122013/1‑90 (MQ=255)
                           gTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGa                                                   <  1:283673/90‑1 (MQ=255)
                                gCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGccc                                              <  1:468759/90‑1 (MQ=255)
                                  aaTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAg                                            >  2:245689/1‑90 (MQ=255)
                                        aGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCAt                                      >  1:422539/1‑90 (MQ=255)
                                        aGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCAt                                      >  1:76888/1‑90 (MQ=255)
                                          gCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGc                                    <  2:42772/90‑1 (MQ=255)
                                               tttCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCccacc                               <  2:404341/90‑1 (MQ=255)
                                                              cGAAAAACAACCTGTCCGGGATTAAAATTACTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCAcc                <  2:290101/90‑1 (MQ=255)
                                                                            tCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGc  >  2:395064/1‑90 (MQ=255)
                                                                                 |                                                                                    
TCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTTGC  >  NZ_CP009273/2661427‑2661592

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