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,0250TC77.3% 32.5 / 8.8 22E126G (GAA→GGA) cynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base C (5/12);  total (5/17)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

TCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAA  >  NZ_CP009273/353941‑354100
                                                                                    |                                                                           
tCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGt                                                                        <  2:202025/90‑1 (MQ=255)
      gAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTgg                                                                  >  1:680999/1‑90 (MQ=255)
           cacaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGt                                                             <  2:510876/90‑1 (MQ=255)
            acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTg                                                            <  2:471502/90‑1 (MQ=255)
            acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                            <  2:636746/90‑1 (MQ=255)
            acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                            <  1:57534/90‑1 (MQ=255)
            acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                            <  1:702142/90‑1 (MQ=255)
                gcgcGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGc                                                        <  1:370805/90‑1 (MQ=255)
                                      cAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATc                                  <  1:809092/90‑1 (MQ=255)
                                          tcgtcgCGGCAAAGCATATCCTGGATTTTCCCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGcc                              <  2:668259/90‑1 (MQ=255)
                                          tcgtcgCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGcc                              <  1:383976/90‑1 (MQ=255)
                                                   cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg                     <  2:297407/90‑1 (MQ=255)
                                                   cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg                     <  2:418418/90‑1 (MQ=255)
                                                   cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg                     <  1:173485/90‑1 (MQ=255)
                                                            tCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAg            >  2:291368/1‑90 (MQ=255)
                                                                      tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg                     >  1:487236/1‑71 (MQ=255)
                                                                      tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg                     <  2:487236/71‑1 (MQ=255)
                                                                      tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaa  >  1:639768/1‑90 (MQ=255)
                                                                      tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaa  <  1:469055/90‑1 (MQ=255)
                                                                      tctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaa  <  1:290942/90‑1 (MQ=255)
                                                                       ctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAg                        >  2:226559/1‑67 (MQ=255)
                                                                       ctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAg                        <  1:226559/67‑1 (MQ=255)
                                                                                    |                                                                           
TCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAA  >  NZ_CP009273/353941‑354100

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