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,9720AG97.5% 122.5 / ‑4.2 40D45G (GAC→GGC) decRDNA‑binding transcriptional regulator DecR
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (23/16);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

ATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAAC  >  NZ_CP009273/463885‑464058
                                                                                       |                                                                                      
aTTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAGACGGCTGGAGGacggcg                                                                                      >  1:174431/1‑90 (MQ=255)
 ttGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGAcggcgg                                                                                     <  2:1550572/90‑1 (MQ=255)
 ttGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGAcggcgg                                                                                     <  1:1516411/90‑1 (MQ=255)
    cACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTAt                                                                                  <  2:746684/90‑1 (MQ=255)
        ctctctTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCtt                                                                              >  2:1157682/1‑90 (MQ=255)
        ctctctTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCtt                                                                              >  2:503902/1‑90 (MQ=255)
         tctctTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTa                                                                             >  2:634938/1‑90 (MQ=255)
          ctctTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTAt                                                                            <  1:132515/90‑1 (MQ=255)
                   ggCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTc                                                                  >  2:677485/1‑90 (MQ=255)
                    gCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTc                                                                  >  1:355262/1‑90 (MQ=255)
                    gCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTc                                                                  >  2:632057/1‑90 (MQ=255)
                    gCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTc                                                                  >  1:271827/1‑90 (MQ=255)
                      aCTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGc                                                                <  2:306551/90‑1 (MQ=255)
                        tGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGccc                                                              >  1:278592/1‑90 (MQ=255)
                        tGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGccc                                                              >  2:1555765/1‑90 (MQ=255)
                             gAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCctgctg                                                         >  1:350296/1‑90 (MQ=255)
                             gAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCctgctg                                                         >  2:1462914/1‑90 (MQ=255)
                                 ccGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATc                                                     <  1:516962/90‑1 (MQ=255)
                                 ccGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATc                                                     <  2:457936/90‑1 (MQ=255)
                                   gTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGc                                                                <  2:1325439/77‑1 (MQ=255)
                                   gTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGc                                                                >  1:1325439/1‑77 (MQ=255)
                                   gTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATc                                                     >  2:861990/1‑88 (MQ=255)
                                   gTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATc                                                     <  1:861990/88‑1 (MQ=255)
                                      aaTCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGaa                                                <  1:546433/90‑1 (MQ=255)
                                      aaTCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGaa                                                <  2:873626/90‑1 (MQ=255)
                                       aTCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGaaa                                               <  1:1572668/90‑1 (MQ=255)
                                           gACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAAt                                           >  2:986219/1‑90 (MQ=255)
                                           gACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAAt                                           >  1:70920/1‑90 (MQ=255)
                                                ccaccCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGcc                                      >  2:482488/1‑90 (MQ=255)
                                                  accCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTc                                    >  2:1295350/1‑90 (MQ=255)
                                                  accCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTc                                    >  2:194601/1‑90 (MQ=255)
                                                  accCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTc                                    >  2:1420468/1‑90 (MQ=255)
                                                      cTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGcc                                <  1:1025701/90‑1 (MQ=255)
                                                      cTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGcc                                <  1:302242/90‑1 (MQ=255)
                                                           tGGAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAcc                           <  1:503902/90‑1 (MQ=255)
                                                             gAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCg                                                                 >  1:63535/1‑50 (MQ=255)
                                                             gAAGCGCCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCg                                                                 <  2:63535/50‑1 (MQ=255)
                                                                  gcCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTtg                    >  2:753958/1‑90 (MQ=255)
                                                                  gcCTGAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTtg                    >  1:645482/1‑90 (MQ=255)
                                                                      gAAACGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCt                >  2:824896/1‑90 (MQ=255)
                                                                          cGGCTGGAGGACGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATa            >  2:490986/1‑90 (MQ=255)
                                                                               ggaggaCGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAAc       <  1:1451406/90‑1 (MQ=255)
                                                                                  ggaCGGCGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCa    <  1:634938/90‑1 (MQ=255)
                                                                                    acggcgGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAAc  <  2:645482/90‑1 (MQ=255)
                                                                                    acggcgGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAAc  <  1:227190/90‑1 (MQ=255)
                                                                                    acggcgGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAAc  <  1:1462914/90‑1 (MQ=255)
                                                                                       |                                                                                      
ATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAAC  >  NZ_CP009273/463885‑464058

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