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,5080TG97.1% 103.4 / ‑3.7 34T250P (ACC→CCC) hcaRDNA‑binding transcriptional regulator HcaR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (15/18);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT  >  NZ_CP009273/2661423‑2661590
                                                                                     |                                                                                  
ggCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTa                                                                                >  2:822016/1‑90 (MQ=255)
 gCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTaa                                                                               <  1:455619/90‑1 (MQ=255)
   ttCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTaaaa                                                                             <  1:925026/90‑1 (MQ=255)
     cATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAatt                                                                           >  1:948775/1‑90 (MQ=255)
      aTCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAatta                                                                          <  2:987006/90‑1 (MQ=255)
       tctcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAattat                                                                         <  1:1097578/90‑1 (MQ=255)
              tttCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCAtata                                                                  <  2:483946/90‑1 (MQ=255)
               ttcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATaa                                                                 >  1:542806/1‑90 (MQ=255)
                tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAAc                                                                <  2:884117/90‑1 (MQ=255)
                tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAAc                                                                <  1:640412/90‑1 (MQ=255)
                  ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                              <  1:623295/90‑1 (MQ=255)
                  ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                              <  2:846036/90‑1 (MQ=255)
                  ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                              <  1:110462/90‑1 (MQ=255)
                  ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCg                                                              <  2:345244/90‑1 (MQ=255)
                    ccACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCggg                                                            >  1:1093101/1‑90 (MQ=255)
                         ccATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCa                                                       >  2:521741/1‑90 (MQ=255)
                           aTCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCaaa                                                     >  1:1578168/1‑90 (MQ=255)
                               gTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGa                                                 <  2:674011/90‑1 (MQ=255)
                               gTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGa                                                 <  2:1130019/90‑1 (MQ=255)
                                 aaaGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACg                                               >  1:1503148/1‑90 (MQ=255)
                                 aaaGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACg                                               >  2:17161/1‑90 (MQ=255)
                                    gCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACACCCTGTCCGGGATTAAAATTATTCATATAACTGGGTATCAAAGTGACGccc                                            <  1:1479447/90‑1 (MQ=255)
                                     cAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCa                                           >  1:1460233/1‑90 (MQ=255)
                                      aaTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAg                                          >  2:1044315/1‑90 (MQ=255)
                                       aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc                                         >  2:1552280/1‑90 (MQ=255)
                                       aTGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc                                         >  2:1485974/1‑90 (MQ=255)
                                            aGGCACATTTCCGGCAATGGGACGAAAAACAACCTTTCCGGGATTAAAATTATTCATAGAACCGGGTATCAAAGTGACGCCCACCCCCAt                                    >  2:1479447/1‑90 (MQ=255)
                                                           aaTGGGACGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCAt                     <  2:962133/90‑1 (MQ=255)
                                                                  cGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCAcc              <  2:1460233/90‑1 (MQ=255)
                                                                  cGAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCAcc              <  2:839980/90‑1 (MQ=255)
                                                                   gAAAAACAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCa             <  2:246119/90‑1 (MQ=255)
                                                                     aaaaCAACCTGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGa           >  1:935478/1‑90 (MQ=255)
                                                                              tGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGtt  >  1:1202669/1‑90 (MQ=255)
                                                                              tGTCCGGGATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGtt  >  1:1502261/1‑90 (MQ=255)
                                                                                     |                                                                                  
GGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCATGGTCACCAGAATATTCGTT  >  NZ_CP009273/2661423‑2661590

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