Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,350,381 G→C G457G (GGC→GGG cadC ← lysine decarboxylation/transport transcriptional activator CadC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,350,3810GC97.1% 109.4 / ‑3.4 34G457G (GGC→GGGcadClysine decarboxylation/transport transcriptional activator CadC
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (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

GTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTT  >  NZ_CP009273/4350298‑4350465
                                                                                   |                                                                                    
gtgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                <  2:37471/90‑1 (MQ=255)
gtgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                <  2:228533/90‑1 (MQ=255)
gtgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                <  1:4318/90‑1 (MQ=255)
 tgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCa                                                                               >  2:474012/1‑90 (MQ=255)
  gtTTGTCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCaa                                                                              >  2:172348/1‑90 (MQ=255)
  gtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCaa                                                                              <  2:619535/90‑1 (MQ=255)
  gtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCaa                                                                              >  2:522469/1‑90 (MQ=255)
  gtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCaa                                                                              <  2:261743/90‑1 (MQ=255)
   gttGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAac                                                                             >  2:644392/2‑90 (MQ=255)
        cccTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATa                                                                        >  2:163428/1‑90 (MQ=255)
          cTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAAt                                                                      >  2:622084/1‑90 (MQ=255)
          cTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAAt                                                                      >  2:604462/1‑90 (MQ=255)
          cTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAAt                                                                      >  2:594214/1‑90 (MQ=255)
                gTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGcc                                                                <  1:169979/90‑1 (MQ=255)
                gTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGcc                                                                >  2:44084/1‑90 (MQ=255)
                 tAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCa                                                               >  1:516469/1‑90 (MQ=255)
                             ggTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAAGACATTTCAAg                                                   >  1:521100/1‑90 (MQ=255)
                                 agaTATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCa                                               <  1:594214/90‑1 (MQ=255)
                                 agaTATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCa                                               <  1:189913/90‑1 (MQ=255)
                                       gCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCa                                         <  2:561355/90‑1 (MQ=255)
                                                  ttCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGc                              <  1:209829/90‑1 (MQ=255)
                                                   tCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGcc                             <  1:280945/90‑1 (MQ=255)
                                                        gTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTa                        <  1:712376/90‑1 (MQ=255)
                                                          tCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAg                      <  1:385187/90‑1 (MQ=255)
                                                           cATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGTCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTa                        >  1:516261/1‑87 (MQ=255)
                                                           cATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGTCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTa                        <  2:516261/87‑1 (MQ=255)
                                                             tCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACt                   >  2:635763/1‑90 (MQ=255)
                                                                 cTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATc               <  2:381143/90‑1 (MQ=255)
                                                                 cTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATc               <  2:27415/90‑1 (MQ=255)
                                                                  ttcatttcatAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCt              >  1:244901/1‑90 (MQ=255)
                                                                     atttcatAAACCTTCCCAAGTAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGtt           >  1:259302/1‑90 (MQ=255)
                                                                             aaCCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGtttt         >  1:233202/1‑84 (MQ=255)
                                                                             aaCCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGtttt         <  2:233202/84‑1 (MQ=255)
                                                                              aCCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGttttacctttt  <  1:474012/90‑1 (MQ=255)
                                                                                   |                                                                                    
GTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTT  >  NZ_CP009273/4350298‑4350465

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