Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,211,170 (T)8→7 coding (1246/2187 nt) espL4 ← type III secretion system effector EspL4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,211,1630T.97.1% 135.8 / ‑1.1 35coding (1253/2187 nt)espL4type III secretion system effector EspL4
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (11/23);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTT  >  NZ_CP009273/4211078‑4211246
                                                                                     |                                                                                   
tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                 >  1:1046214/1‑90 (MQ=255)
    acaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcg                                                                             <  1:1445499/90‑3 (MQ=255)
    acaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcg                                                                             <  2:722021/90‑3 (MQ=255)
      aTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcgag                                                                           <  1:934916/90‑5 (MQ=255)
      aTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcgag                                                                           >  2:195905/1‑86 (MQ=255)
           tATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                                     <  1:1294825/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                                   <  2:865398/90‑1 (MQ=255)
               cTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCt                                                                 >  1:549017/1‑90 (MQ=255)
                 tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCtt                                                                >  2:1216123/1‑90 (MQ=255)
                 tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTc                                                               >  2:828936/1‑90 (MQ=255)
                      gATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGA‑‑TTTTTTCGAGAGTAATCCTTCATATcc                                                         <  1:387816/90‑1 (MQ=255)
                          tCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGt                                                      >  2:108728/1‑90 (MQ=255)
                           cGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGtt                                                     >  2:1484146/1‑90 (MQ=255)
                               aaaaaTAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTg                                                 >  2:1302112/1‑90 (MQ=255)
                                   ataatCCAGAAGAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACaa                                             <  1:1356479/90‑1 (MQ=255)
                                   ataatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACaa                                             <  2:890052/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                           <  2:288370/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                           >  1:1341045/1‑90 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                           <  2:868005/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCAAGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                           <  1:45906/90‑1 (MQ=255)
                                      atCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGa                                          >  2:215043/1‑90 (MQ=255)
                                      atCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGa                                          >  2:120098/1‑90 (MQ=255)
                                       tCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGaa                                         >  1:1161632/1‑90 (MQ=255)
                                           gAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGa                                     >  1:1321413/1‑90 (MQ=255)
                                               aaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAAt                                 <  1:637077/90‑1 (MQ=255)
                                                 ccATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAg                               <  1:1454686/90‑1 (MQ=255)
                                                    tttttATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGtct                            <  2:1572825/90‑1 (MQ=255)
                                                    tttttATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGtct                            <  1:1506299/90‑1 (MQ=255)
                                                       ttATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCaa                         <  2:1341045/90‑1 (MQ=255)
                                                        tATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAc                        <  1:448029/90‑1 (MQ=255)
                                                        tATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAc                        <  2:486893/90‑1 (MQ=255)
                                                              ttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATc                  <  1:1458431/90‑1 (MQ=255)
                                                               tGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATcc                 <  2:618832/90‑1 (MQ=255)
                                                                 ccTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGc               <  1:622287/90‑1 (MQ=255)
                                                                   tCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAt             <  1:1302112/90‑1 (MQ=255)
                                                                    ccAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAtt            <  2:1345559/90‑1 (MQ=255)
                                                                    ccAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAtt            <  1:231698/90‑1 (MQ=255)
                                                                         atatGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCAt       <  2:767330/90‑1 (MQ=255)
                                                                          tatGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATa      <  2:1435944/90‑1 (MQ=255)
                                                                              cATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAAttt  >  2:95294/1‑90 (MQ=255)
                                                                                     |                                                                                   
TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTT  >  NZ_CP009273/4211078‑4211246

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