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,0250TC97.4% 120.7 / ‑4.5 39E126G (GAA→GGA) cynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (16/22);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGG  >  NZ_CP009273/353939‑354106
                                                                                      |                                                                                 
gCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCt                                                                                <  1:912766/90‑1 (MQ=255)
gCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCt                                                                                <  1:1223535/90‑1 (MQ=255)
gCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCt                                                                                <  2:148380/90‑1 (MQ=255)
  tCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGt                                                                              <  2:56153/90‑1 (MQ=255)
   ccGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCTCTGTa                                                                             >  1:1121979/1‑90 (MQ=255)
   ccGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTa                                                                             >  1:56773/1‑90 (MQ=255)
   ccGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTa                                                                             >  1:188942/1‑90 (MQ=255)
     ggCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGc                                                                           >  1:1105288/1‑90 (MQ=255)
             cacaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGt                                                                   <  1:325799/90‑1 (MQ=255)
             cacaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGACAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGt                                                                   <  1:293582/90‑1 (MQ=255)
              acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                                  <  2:900684/90‑1 (MQ=255)
              acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                                  <  1:365834/90‑1 (MQ=255)
              acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                                  <  1:107716/90‑1 (MQ=255)
                  gcgcGAAGGCAATCCCAGCGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGc                                                              <  2:613089/90‑1 (MQ=255)
                  gcgcGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAg                                                                      <  1:1558440/82‑1 (MQ=255)
                  gcgcGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAg                                                                      >  2:1558440/1‑82 (MQ=255)
                            aaTCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAgcg                                                    <  1:820138/90‑1 (MQ=255)
                              tCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAgcgcg                                                  <  2:1105288/90‑1 (MQ=255)
                                       ccAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAAt                                         >  2:927714/1‑90 (MQ=255)
                                       ccAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAAt                                         >  2:448436/1‑90 (MQ=255)
                                         aaCTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATcc                                       <  1:598313/90‑1 (MQ=255)
                                           cTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGc                                     <  2:1123331/90‑1 (MQ=255)
                                             cgtcgCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTgg                                                                        >  1:148692/1‑53 (MQ=255)
                                             cgtcgCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTgg                                                                        <  2:148692/53‑1 (MQ=255)
                                              gtcgCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCAt                                  >  1:1487336/1‑90 (MQ=255)
                                                 gcgGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTaa                               <  1:927714/90‑1 (MQ=255)
                                                     cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg                           <  2:1305907/90‑1 (MQ=255)
                                                            tatCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATaa                    >  1:213768/1‑90 (MQ=255)
                                                                 tagATTTTCACCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAg               <  1:1315876/88‑1 (MQ=255)
                                                                 tCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAg               <  2:213768/90‑1 (MQ=255)
                                                                 tCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAg               >  2:275572/1‑90 (MQ=255)
                                                                 tCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAg               >  1:552717/1‑90 (MQ=255)
                                                                   gATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGcc                         >  1:896454/1‑78 (MQ=255)
                                                                   gATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGcc                         <  2:896454/78‑1 (MQ=255)
                                                                   gATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCt             >  1:984059/1‑90 (MQ=255)
                                                                       ttCTCCTGCGACATTCTCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTa         <  2:1121979/90‑1 (MQ=255)
                                                                       ttCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTa         >  1:207216/1‑90 (MQ=255)
                                                                        tctcCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaa        <  1:1132039/90‑1 (MQ=255)
                                                                              gCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTgg  >  1:77088/1‑90 (MQ=255)
                                                                                      |                                                                                 
GCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGG  >  NZ_CP009273/353939‑354106

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