Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,039,724 A→G T436A (ACA→GCA)  yeeJ → inverse autotransporter adhesin YeeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,039,7240AG95.7% 76.5 / ‑4.3 23T436A (ACA→GCA) yeeJinverse autotransporter adhesin YeeJ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (13/9);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGCAATGCAGAAACAGCTTGACCCGAATGAAGTCGCTGCACGGCGTAGCCTTGCAGGCAGCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGACAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGT  >  NZ_CP009273/2039588‑2039848
                                                                                                                                        |                                                                                                                            
cgcAATGCAGAAACAGCTTGACCCGAATGAAGTCGCTGCACGGCGTAGCCTTGCAGGCAGCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCa                                                                                                                            <  2:423577/139‑1 (MQ=255)
                cTTGACCCGAATGAAGTCGCTGCACGGCGTAGCCTTGCAGGCAGCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAgg                                                                                                              >  1:117540/1‑137 (MQ=255)
                cTTGACCCGAATGAAGTCGCTGCACGGCGTAGCCTTGCAGGCAGCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAgg                                                                                                              <  2:117540/137‑1 (MQ=255)
                        gAATGAAGTCGCTGCACGGCGTAGCCTTGCAGGTAGCCGTTATGATCTGGTGGATCGCACCAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGga                                                                                                    <  1:197866/139‑1 (MQ=255)
                                                   tGCAGGCAGCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTa                                                                         >  2:440228/1‑139 (MQ=255)
                                                     caggcagCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTAc                                                                        >  1:349355/1‑138 (MQ=255)
                                                     caggcagCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTAc                                                                        >  2:434074/1‑138 (MQ=255)
                                                          agCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAAcc                                                                   >  1:421875/1‑138 (MQ=255)
                                                                         gtggATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTg                                                       <  1:187751/135‑1 (MQ=255)
                                                                         gtggATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTg                                                       <  1:174178/135‑1 (MQ=255)
                                                                         gtggATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAg                                                   >  1:17948/1‑139 (MQ=255)
                                                                                gCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCtata                                             <  2:349355/138‑1 (MQ=255)
                                                                                                  tGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCa                                                                                                                            <  2:438979/41‑1 (MQ=255)
                                                                                                  tGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCa                                                                                                                            >  1:438979/1‑41 (MQ=255)
                                                                                                    gAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTg                         >  2:226666/1‑138 (MQ=255)
                                                                                                         tCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCtata                                             <  1:183418/113‑1 (MQ=255)
                                                                                                         tCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCtata                                             >  2:183418/1‑113 (MQ=255)
                                                                                                         tCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCt                   >  1:178780/1‑139 (MQ=255)
                                                                                                         tCGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCt                   >  1:272625/1‑139 (MQ=255)
                                                                                                          cGCAAAAAAGAACTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTg                  >  2:226664/1‑139 (MQ=255)
                                                                                                                     aCTGGTTCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCaaa       <  1:440228/139‑1 (MQ=255)
                                                                                                                     aCTGGTTCGCCTGACCCTGACAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCaaa       <  2:74493/139‑1 (MQ=255)
                                                                                                                           tCGCCTGACCCTGGCAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGt  >  1:157516/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                            
CGCAATGCAGAAACAGCTTGACCCGAATGAAGTCGCTGCACGGCGTAGCCTTGCAGGCAGCCGTTATGATCTGGTGGATCGCAACAACAATATCGTTCTGGAATATCGCAAAAAAGAACTGGTTCGCCTGACCCTGACAGACCCCGTGACAGGGAAGTCAGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGT  >  NZ_CP009273/2039588‑2039848

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