Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,458,386 A→G E369G (GAA→GGA)  gspL → type II secretion system protein GspL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,458,3860AG95.7% 68.5 / ‑3.3 23E369G (GAA→GGA) gspLtype II secretion system protein GspL
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (14/8);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAT  >  NZ_CP009273/3458303‑3458468
                                                                                   |                                                                                  
gTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAAc                                                                              <  1:190656/90‑1 (MQ=255)
      gACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGaaaa                                                                        <  1:249929/90‑1 (MQ=255)
                aTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                              <  1:206631/90‑1 (MQ=255)
                aTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCCTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                              <  2:455479/90‑1 (MQ=255)
                  gAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAg                                                            >  2:218375/1‑90 (MQ=255)
                  gAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAg                                                            >  1:467599/1‑90 (MQ=255)
                  gAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAg                                                            >  2:198076/1‑90 (MQ=255)
                    gTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                          <  2:238653/90‑1 (MQ=255)
                         aaGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTa                                                     <  1:151749/90‑1 (MQ=255)
                              tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                <  1:218375/90‑1 (MQ=255)
                              tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                >  2:177921/1‑90 (MQ=255)
                                           gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCa                                   >  1:47511/1‑90 (MQ=255)
                                                 ttCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                >  2:51917/1‑71 (MQ=255)
                                                 ttCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                <  1:51917/71‑1 (MQ=255)
                                                                 aGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtggtg             >  2:181030/1‑90 (MQ=255)
                                                                 aGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtggtg             >  1:375729/1‑90 (MQ=255)
                                                                  gTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                              <  1:459820/40‑1 (MQ=39)
                                                                  gTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                              >  2:459820/1‑40 (MQ=39)
                                                                  gTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtggtgg            >  2:211381/1‑90 (MQ=255)
                                                                            ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt  >  2:147874/1‑90 (MQ=255)
                                                                            ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt  >  1:429927/1‑90 (MQ=255)
                                                                            ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt  >  1:256857/1‑90 (MQ=255)
                                                                            ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGGATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt  >  2:423375/1‑90 (MQ=255)
                                                                                   |                                                                                  
GTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAT  >  NZ_CP009273/3458303‑3458468

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