Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,972,615 A→G D79G (GAT→GGT)  lysR → DNA‑binding transcriptional regulator LysR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,972,6150AG78.1% 57.3 / 13.3 32D79G (GAT→GGT) lysRDNA‑binding transcriptional regulator LysR
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base G (12/13);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

CGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCAATCTTT  >  NZ_CP009273/2972548‑2972702
                                                                   |                                                                                       
cGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGaa                                                                   <  2:119713/90‑1 (MQ=255)
  aTTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAg                                                                 >  1:1074899/1‑90 (MQ=255)
        tCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTgcg                                                           >  1:733598/1‑90 (MQ=255)
          cTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGcgcg                                                         <  1:497739/90‑1 (MQ=255)
            aCCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGAGCCGCAGAAAGTCTGCGCGAg                                                       >  2:14254/1‑90 (MQ=255)
            aCCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCGTTGTCAGGTACGCAGAAAGTCTGCGCGAg                                                       >  2:521236/1‑90 (MQ=255)
            aCCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAg                                                       >  1:87128/1‑90 (MQ=255)
               gTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGttt                                                    >  1:430120/1‑90 (MQ=255)
               gTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGttt                                                    >  1:537213/1‑90 (MQ=255)
                tGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTc                                                   >  2:472647/1‑90 (MQ=255)
                  cAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGCGGCGCAGAAAGTCTGCGCGAGTTTCGc                                                 <  1:74709/90‑1 (MQ=255)
                  cAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGc                                                 <  2:507861/90‑1 (MQ=255)
                       aCTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAggg                                            <  1:319171/90‑1 (MQ=255)
                             tCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACt                                      <  2:735632/90‑1 (MQ=255)
                                         agTGCAACGATCCTGGTACGGACTGGGTCGCATTGTCAGAGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATtgcctg                          <  1:14254/90‑1 (MQ=255)
                                              aaCGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAAc                                       <  2:791999/72‑1 (MQ=255)
                                              aaCGATCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAAc                                       >  1:791999/1‑72 (MQ=255)
                                                  aTCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACt                                      <  1:645437/69‑1 (MQ=255)
                                                  aTCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACt                                      >  2:645437/1‑69 (MQ=255)
                                                  aTCCTGGTACGGACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTCCCTGCCTTCCGGt                 >  2:112749/1‑90 (MQ=255)
                                                            ggACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGc         <  1:658239/88‑1 (MQ=255)
                                                            ggACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGc         >  2:658239/1‑88 (MQ=255)
                                                            ggACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCaa       <  1:682663/90‑1 (MQ=255)
                                                            ggACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCaa       <  2:173609/90‑1 (MQ=255)
                                                            ggACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCaa       <  1:481381/90‑1 (MQ=255)
                                                            ggACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCaa       <  1:472647/90‑1 (MQ=255)
                                                            ggACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCaa       >  1:173609/1‑90 (MQ=255)
                                                             gACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCAAt      >  1:78010/1‑90 (MQ=255)
                                                             gACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCAAt      >  1:588174/1‑90 (MQ=255)
                                                             gACTGGGTCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCAAt      >  1:131812/1‑90 (MQ=255)
                                                             gACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCAAt      >  2:355171/1‑90 (MQ=255)
                                                                 ggATCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCAATCttt  >  1:39650/1‑90 (MQ=255)
                                                                   |                                                                                       
CGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGCGCGAGTTTCGCCAGGGAGAACTGTCTATTGCCTGCCTGCCGGTCTTTTCGCAATCTTT  >  NZ_CP009273/2972548‑2972702

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