Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 463,972 A→G D45G (GAC→GGC)  decR → DNA‑binding transcriptional regulator DecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273463,9720AG79.0% 29.7 / 5.0 19D45G (GAC→GGC) decRDNA‑binding transcriptional regulator DecR
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base G (9/6);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAAC  >  NZ_CP009273/463886‑464058
                                                                                      |                                                                                      
ttGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCCCCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGAcggcgg                                                                                     <  1:13036/90‑1 (MQ=255)
ttGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGAcggcgg                                                                                     <  2:87531/90‑1 (MQ=255)
  gCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTa                                                                                   >  1:431266/1‑90 (MQ=255)
       ctctctTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCtt                                                                              >  1:620017/1‑90 (MQ=255)
        tctctTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTa                                                                             >  2:820820/1‑90 (MQ=255)
                   gCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTc                                                                  >  2:97861/1‑90 (MQ=255)
                   gCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTc                                                                  >  1:552196/1‑90 (MQ=255)
                   gCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTc                                                                  >  1:598490/1‑90 (MQ=255)
                       tGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGccc                                                              >  2:806909/1‑90 (MQ=255)
                                  gTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCg                                                   >  1:720750/1‑90 (MQ=255)
                                     aaTCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGACCCGGaa                                                <  1:806909/90‑1 (MQ=255)
                                           aCAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATa                                          <  2:254972/90‑1 (MQ=255)
                                           aCAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATa                                          <  1:585241/90‑1 (MQ=255)
                                               ccaccCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGcc                                      >  2:614677/1‑90 (MQ=255)
                                                 accCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTc                                    >  2:95652/1‑90 (MQ=255)
                                                      ttGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCt                               >  1:346812/1‑90 (MQ=255)
                                                          tGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAcc                           <  2:552196/90‑1 (MQ=255)
                                                          tGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAcc                           <  1:482650/90‑1 (MQ=255)
                                                             aaGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAcc                           >  2:533122/1‑87 (MQ=255)
                                                             aaGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAcc                           <  1:533122/87‑1 (MQ=255)
                                                                 gcCTGAAACGGCTGGAGGACGACGGGATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTtg                    >  2:743803/1‑90 (MQ=255)
                                                                                   acggcgGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAAc  <  2:370481/90‑1 (MQ=255)
                                                                                      |                                                                                      
TTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAAC  >  NZ_CP009273/463886‑464058

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