Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 354,025 T→C E126G (GAA→GGA)  cynR ← transcriptional regulator CynR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273354,0250TC87.5% 32.5 / ‑1.2 16E126G (GAA→GGA) cynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (7/7);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAA  >  NZ_CP009273/353943‑354101
                                                                                  |                                                                            
cGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAg                                                                       >  2:459269/1‑90 (MQ=255)
    gAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGCCATTCCCTGTAGCTgg                                                                   <  1:277513/90‑1 (MQ=255)
         cacaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGt                                                              <  2:110111/90‑1 (MQ=255)
          acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTg                                                             <  1:375795/90‑1 (MQ=255)
          acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                             <  1:385973/90‑1 (MQ=255)
          acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                             <  2:64725/90‑1 (MQ=255)
           caGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCt                                                                           >  1:292684/1‑75 (MQ=255)
           caGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCt                                                                           <  2:292684/75‑1 (MQ=255)
                  gAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTggg                                                     >  1:220291/1‑90 (MQ=255)
                  gAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTggg                                                     >  2:262390/1‑90 (MQ=255)
                                   ccAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAAt                                    >  2:273469/1‑90 (MQ=255)
                                     aaCTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATcc                                  <  1:499508/90‑1 (MQ=255)
                                                 cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg                      <  2:143552/90‑1 (MQ=255)
                                                   aaGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCACAGAAATCCGCCCTTAAGGGGct                    >  1:553454/1‑89 (MQ=255)
                                                             tCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAg          >  2:185372/1‑90 (MQ=255)
                                                                     ctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaaa  <  2:186083/90‑1 (MQ=255)
                                                                                  |                                                                            
CGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAA  >  NZ_CP009273/353943‑354101

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