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.8% 145.4 / ‑1.0 48E369G (GAA→GGA) gspLtype II secretion system protein GspL
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (22/24);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCC  >  NZ_CP009273/3458301‑3458470
                                                                                     |                                                                                    
aaGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGgaaaa                                                                                  <  1:776337/90‑1 (MQ=255)
   tCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACa                                                                               <  2:471414/90‑1 (MQ=255)
                  aTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                                <  2:191737/90‑1 (MQ=255)
                  aTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                                <  2:519857/90‑1 (MQ=255)
                    gAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAg                                                              >  2:373545/1‑90 (MQ=255)
                    gAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAg                                                              >  2:459941/1‑90 (MQ=255)
                    gAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAg                                                              >  2:530080/1‑90 (MQ=255)
                    gAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAg                                                              >  1:103902/1‑90 (MQ=255)
                      gTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                            <  1:399263/90‑1 (MQ=255)
                      gTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                            <  1:271726/90‑1 (MQ=255)
                      gTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                            <  2:120037/90‑1 (MQ=255)
                      gTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                            <  1:632158/90‑1 (MQ=255)
                      gTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAgtg                                                            <  1:641238/90‑1 (MQ=255)
                           aaGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTa                                                       <  2:15852/90‑1 (MQ=255)
                                tGAAGCTACTATCGATCGTTACTGTGAGTTAATACAGTCATGGCTTCCGATGGGAAAACCAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  2:582400/90‑1 (MQ=255)
                                tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  1:275644/90‑1 (MQ=255)
                                tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  >  1:432434/1‑90 (MQ=255)
                                tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  2:534830/90‑1 (MQ=255)
                                tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  2:601585/90‑1 (MQ=255)
                                tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  1:530080/90‑1 (MQ=255)
                                tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  1:680177/90‑1 (MQ=255)
                                  aaGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAg                                                >  1:777022/1‑90 (MQ=255)
                                         tatCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                            >  1:81923/1‑71 (MQ=255)
                                         tatCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAgtg                                                            <  2:81923/71‑1 (MQ=255)
                                         tatCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGaa                                         <  1:564064/90‑1 (MQ=255)
                                             gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCa                                     >  1:108774/1‑90 (MQ=255)
                                             gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCa                                     >  1:453205/1‑90 (MQ=255)
                                             gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCa                                     >  2:253615/1‑90 (MQ=255)
                                             gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAATGGTGTATGGACAGTAAGGAACTCa                                     >  2:200983/1‑90 (MQ=255)
                                              atcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAg                                    >  1:568441/1‑90 (MQ=255)
                                              atcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAg                                    >  1:666378/1‑90 (MQ=255)
                                              atcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAg                                    >  2:793823/1‑90 (MQ=255)
                                               tcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAgg                                   >  2:718722/1‑90 (MQ=255)
                                               tcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAgg                                   >  1:211520/1‑90 (MQ=255)
                                                   ttCTGTGAGTTAACACAATCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGaaa                               <  1:169122/90‑1 (MQ=255)
                                                       gtgAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGAt                           >  1:466241/1‑90 (MQ=255)
                                                          aGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTaa                        >  1:402547/1‑90 (MQ=255)
                                                            ttAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAAt                      <  1:459941/90‑1 (MQ=255)
                                                            ttAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAATGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAAt                      <  1:200983/90‑1 (MQ=255)
                                                               acacAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAt                   >  1:420508/1‑90 (MQ=255)
                                                                 acaGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAGATGATTAAATCAtgg                 <  2:453205/90‑1 (MQ=255)
                                                                 acaGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtgg                 <  1:765876/90‑1 (MQ=255)
                                                                   aGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtggtg               >  1:594161/1‑90 (MQ=255)
                                                                          gCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGaa        <  1:235189/90‑1 (MQ=255)
                                                                              ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  >  2:594832/1‑44 (MQ=255)
                                                                              ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  1:594832/44‑1 (MQ=255)
                                                                              ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt    >  2:691329/1‑90 (MQ=255)
                                                                                gATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATcc  >  2:344594/1‑90 (MQ=255)
                                                                                gATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATcc  <  1:58455/90‑1 (MQ=255)
                                                                                     |                                                                                    
AAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCC  >  NZ_CP009273/3458301‑3458470

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