Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,452,746 G→T G350V (GGC→GTC)  gspE → type II secretion system ATPase GspE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,452,7460GT92.0% 169.8 / 8.5 50G350V (GGC→GTC) gspEtype II secretion system ATPase GspE
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (23/23);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

CGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCC  >  NZ_CP009273/3452610‑3452876
                                                                                                                                        |                                                                                                                                  
cgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGtcg                                                                                                                                  >  1:888853/1‑139 (MQ=255)
cgcgcCATACTTCGCAAGGACCCGGATGTCGTCATGGTGGGGGAAATGCGTGATACAGAAGCCGCACAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTCCATCAgg                                                                                                                                    <  2:513199/137‑1 (MQ=255)
  cgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCa                                                                                                                                <  1:736244/139‑1 (MQ=255)
               cAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTcc                                                                                                                   <  2:677557/139‑1 (MQ=255)
                  gACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCcgc                                                                                                                 >  1:1054011/1‑138 (MQ=255)
                        gATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATgg                                                                                                          >  1:1116597/1‑139 (MQ=255)
                            tcgtcATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGt                                                                                                      >  2:249689/1‑139 (MQ=255)
                            tcgtcATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGt                                                                                                      <  2:344402/139‑1 (MQ=255)
                            tcgtcATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGt                                                                                                      >  1:344402/1‑139 (MQ=255)
                                      gggggAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTc                                                                                            >  1:270005/1‑139 (MQ=255)
                                       ggggAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTc                                                                                            <  1:249689/138‑1 (MQ=255)
                                         ggAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTg                                                                                         >  2:1181901/1‑139 (MQ=255)
                                             aTTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCttt                                                                                     >  1:1160276/1‑139 (MQ=255)
                                                          aaaCCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGc                                                                        <  1:237804/139‑1 (MQ=255)
                                                               gcgcAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGa                                                                   >  1:627918/1‑139 (MQ=255)
                                                                 gcAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGAtt                                                                 <  2:888853/139‑1 (MQ=255)
                                                                  cAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTa                                                                >  1:638299/1‑139 (MQ=255)
                                                                        gCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcgcg                                                           <  2:48789/138‑1 (MQ=255)
                                                                        gCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcgcg                                                           <  1:1011893/138‑1 (MQ=255)
                                                                         cGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCa                                                         <  1:125908/139‑1 (MQ=255)
                                                                               aGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTc                                                    <  2:1116597/138‑1 (MQ=255)
                                                                               aGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCt                                                   >  1:570439/1‑139 (MQ=255)
                                                                                      gCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGa                                                                   >  2:401862/1‑116 (MQ=255)
                                                                                      gCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGa                                                                   <  1:401862/116‑1 (MQ=255)
                                                                                      gCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTtcgt                                            >  1:232657/1‑139 (MQ=255)
                                                                                          aCCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTtcgtcg                                          <  2:627918/137‑1 (MQ=255)
                                                                                           ccGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcgcg                                                           >  1:87810/1‑119 (MQ=255)
                                                                                           ccGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcgcg                                                           <  2:1267943/119‑1 (MQ=255)
                                                                                           ccGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcgcg                                                           <  2:87810/119‑1 (MQ=255)
                                                                                           ccGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcgcg                                                           >  1:1267943/1‑119 (MQ=255)
                                                                                                cATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTc                                  <  1:1080530/139‑1 (MQ=255)
                                                                                                  tCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGc                                                                        <  1:686622/99‑1 (MQ=255)
                                                                                                  tCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGc                                                                        >  2:686622/1‑99 (MQ=255)
                                                                                                            tcAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTgg                                                                         <  1:883663/88‑1 (MQ=255)
                                                                                                            tcAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTgg                                                                         >  2:883663/1‑88 (MQ=255)
                                                                                                            tcAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGcc                                        <  2:138049/121‑1 (MQ=255)
                                                                                                            tcAACACTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGcc                                        >  1:138049/1‑121 (MQ=255)
                                                                                                                cacTCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAAtt                  <  2:570439/139‑1 (MQ=255)
                                                                                                                   tCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCAc               >  1:39305/1‑139 (MQ=255)
                                                                                                                   tCCACACTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCAc               >  1:258832/1‑139 (MQ=255)
                                                                                                                     cacaCTAACAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTACCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGc             >  1:528340/1‑139 (MQ=255)
                                                                                                                            aCAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTAt       <  2:648640/138‑1 (MQ=255)
                                                                                                                            aCAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATc      <  2:270005/139‑1 (MQ=255)
                                                                                                                             cAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTcc                                                                             >  2:854500/1‑67 (MQ=255)
                                                                                                                             cAGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTcc                                                                             <  1:854500/67‑1 (MQ=255)
                                                                                                                             cAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCa     >  2:249890/1‑139 (MQ=255)
                                                                                                                              aGTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCAc    <  1:191113/139‑1 (MQ=255)
                                                                                                                               gTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGc             <  2:1018335/129‑1 (MQ=255)
                                                                                                                               gTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGc             >  1:1018335/1‑129 (MQ=255)
                                                                                                                               gTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCAcc   >  1:398364/1‑139 (MQ=255)
                                                                                                                               gTGCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCAcc   >  2:147619/1‑139 (MQ=255)
                                                                                                                                 gCATCAGTCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCAccc  >  2:368121/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                  
CGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCC  >  NZ_CP009273/3452610‑3452876

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