Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,202,478 T→G Y345D (TAT→GAT)  yehP → VWA domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,202,4780TG94.1% 48.6 / ‑3.8 17Y345D (TAT→GAT) yehPVWA domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (9/7);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTT  >  NZ_CP009273/2202402‑2202561
                                                                            |                                                                                   
tCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATAc                                                                        <  1:421890/90‑1 (MQ=255)
      gAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCa                                                                  <  2:525316/90‑1 (MQ=255)
           aGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGgcgc                                                             <  1:255380/90‑1 (MQ=255)
           aGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGgcgc                                                             <  2:213410/90‑1 (MQ=255)
             tgtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGGCGCTg                                                           >  2:27337/1‑90 (MQ=255)
             tgtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGGCGCTg                                                           >  2:506331/1‑90 (MQ=255)
                     gagCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGGCGCTGGTTAATGt                                                   >  2:485682/1‑90 (MQ=255)
                      agCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGgcgc                                                             >  1:524655/1‑79 (MQ=255)
                      agCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGgcgc                                                             <  2:524655/79‑1 (MQ=255)
                          gCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGgcgc                                                             <  1:510647/75‑1 (MQ=255)
                          gCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGgcgc                                                             >  2:510647/1‑75 (MQ=255)
                                    gCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAgc                                    >  2:64354/1‑90 (MQ=255)
                                           cTGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATg                             <  2:343417/90‑1 (MQ=255)
                                            tGGCAGCGCTCGATAGCACCGCAACACCTTGCGATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGa                            >  1:294185/1‑90 (MQ=255)
                                                              ccGCAACACCTTGCGATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcat          >  1:308351/1‑90 (MQ=255)
                                                                gCAACACCTTGCTATGCCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatca        <  1:197853/90‑1 (MQ=255)
                                                                      ccTTGCGATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCtt  >  1:258498/1‑90 (MQ=255)
                                                                            |                                                                                   
TCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTT  >  NZ_CP009273/2202402‑2202561

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