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,3860AG100.0% 66.4 / NA 21E369G (GAA→GGA) gspLtype II secretion system protein GspL
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/8);  total (13/8)

CGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAT  >  NZ_CP009273/3458305‑3458468
                                                                                 |                                                                                  
cGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAg                                                                            >  1:202838/1‑90 (MQ=255)
    gACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGTGAAAAACAGaaaa                                                                        <  2:258424/90‑1 (MQ=255)
         tAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATc                                                                   <  1:230254/90‑1 (MQ=255)
         tAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATc                                                                   >  2:202463/1‑90 (MQ=255)
              aTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                              <  1:55158/90‑1 (MQ=255)
                          agTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATgg                                                  <  2:316930/90‑1 (MQ=255)
                            tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                          <  1:232064/80‑1 (MQ=255)
                            tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                          >  2:232064/1‑80 (MQ=255)
                            tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                >  1:208782/1‑90 (MQ=255)
                              aaGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAg                                              <  2:222973/90‑1 (MQ=255)
                                         gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCa                                   >  1:146082/1‑90 (MQ=255)
                                             gTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                          >  2:100752/1‑63 (MQ=255)
                                             gTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                          <  1:100752/63‑1 (MQ=255)
                                               ttCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGaaa                             >  2:201943/1‑90 (MQ=255)
                                                               aGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtggtg             >  1:289196/1‑90 (MQ=255)
                                                                   cTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCa         >  1:94960/1‑90 (MQ=255)
                                                                     ggCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGa       >  2:270068/1‑90 (MQ=255)
                                                                          ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtgg               >  2:153544/1‑77 (MQ=255)
                                                                          ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtgg               <  1:153544/77‑1 (MQ=255)
                                                                          ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt  >  2:132308/1‑90 (MQ=255)
                                                                          ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt  >  1:237480/1‑90 (MQ=255)
                                                                                 |                                                                                  
CGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAT  >  NZ_CP009273/3458305‑3458468

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