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.92.9% 44.9 / ‑1.0 14coding (1253/2187 nt)espL4type III secretion system effector EspL4
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (5/8);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTA  >  NZ_CP009273/4211078‑4211237
                                                                                     |                                                                          
tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                        >  1:149963/1‑90 (MQ=255)
    acaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcg                                                                    <  1:265092/90‑3 (MQ=255)
    acaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcg                                                                    <  2:328475/90‑3 (MQ=255)
           tATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCCCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                            <  1:323352/90‑1 (MQ=255)
            aTCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAAt                                                           >  1:228150/1‑90 (MQ=255)
            aTCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAAt                                                           >  2:65744/1‑90 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                          <  1:171699/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                          <  2:50070/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                          <  1:63717/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                          <  1:288839/90‑1 (MQ=255)
               cTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATcc                                                         >  1:243345/1‑90 (MQ=255)
               cTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCt                                                        >  2:254885/1‑90 (MQ=255)
                                      atCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGa                                 >  2:363847/1‑90 (MQ=255)
                                                            ttttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata           <  1:254885/90‑1 (MQ=255)
                                                                   tCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAt    <  2:354061/90‑1 (MQ=255)
                                                                    ccAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAtt   <  2:82442/90‑1 (MQ=255)
                                                                     cAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTa  >  2:45522/1‑90 (MQ=255)
                                                                                     |                                                                          
TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTA  >  NZ_CP009273/4211078‑4211237

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