Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,134,272 A→G T147A (ACC→GCC)  flgK → flagellar hook‑associated protein FlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,134,2720AG91.7% 103.8 / 2.0 36T147A (ACC→GCC) flgKflagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (14/19);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

TACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGT  >  NZ_CP009273/1134139‑1134399
                                                                                                                                     |                                                                                                                               
tACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGAAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGAt                                                                                                                            <  2:618284/139‑1 (MQ=255)
tACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGCCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGAt                                                                                                                            <  2:793676/139‑1 (MQ=255)
    tctTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCaa                                                                                                                         >  2:590193/1‑138 (MQ=255)
        cGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCt                                                                                                                    >  1:1276460/1‑139 (MQ=255)
               aaCACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGGGAGTAACGCGGAAGACCGGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGAc                                                                                                             <  1:1261034/139‑1 (MQ=255)
                acacAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGAc                                                                                                             >  1:590064/1‑138 (MQ=255)
                         cAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACaaa                                                                                                    <  1:1131183/138‑1 (MQ=255)
                           ggATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAg                                                                                                 >  1:710718/1‑139 (MQ=255)
                                tcttcACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCaa                                                                                            >  1:146662/1‑139 (MQ=255)
                                     aCCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATcg                                                                                       <  2:343058/139‑1 (MQ=255)
                                        aGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACaaa                                                                                                    <  2:177687/123‑1 (MQ=255)
                                        aGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACaaa                                                                                                    >  1:177687/1‑123 (MQ=255)
                                         gCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGAt                                                                                   <  1:500521/139‑1 (MQ=255)
                                          ccTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGa                                                                                    >  2:500521/1‑137 (MQ=255)
                                               aaaCGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGc                                                                             <  2:146662/139‑1 (MQ=255)
                                                    cTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGt                                                                                               >  1:884633/1‑116 (MQ=255)
                                                    cTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGt                                                                                               <  2:884633/116‑1 (MQ=255)
                                                      ggTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGAc                                                                                                             <  1:1087175/100‑1 (MQ=255)
                                                      ggTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGAc                                                                                                             >  2:1087175/1‑100 (MQ=255)
                                                          aGTAACGCGGAAGCCCCGGCACCGCGCCAGGCGCTAATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgat                                                                   <  1:369848/138‑1 (MQ=255)
                                                          aGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgatc                                                                  <  2:349813/139‑1 (MQ=255)
                                                          aGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGAAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgatc                                                                  <  2:737584/139‑1 (MQ=255)
                                                               cgcgGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgatcagat                                                              >  2:830958/1‑138 (MQ=255)
                                                                 cgGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATcaa                                                           >  2:298752/1‑139 (MQ=255)
                                                                    aaGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATcaacaa                                                        <  2:991878/139‑1 (MQ=255)
                                                                    aaGACCCGGCAGCGCGCCAGGCGCTGATGGGGAAATCAGAAGGTTTGGTGAATCAGTTAAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATcaacaa                                                        <  2:729004/139‑1 (MQ=255)
                                                                      gACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGAt                                                                                                                            <  1:513702/69‑1 (MQ=255)
                                                                      gACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGAt                                                                                                                            >  2:513702/1‑69 (MQ=255)
                                                                                  cgcCAGGCGCTGATTGGTAAATCAGAAGGATTGGTGAACCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACaaa                                                                                                    <  1:336603/81‑1 (MQ=255)
                                                                                  cgcCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACaaa                                                                                                    >  2:336603/1‑81 (MQ=255)
                                                                                             gATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCt                                >  1:1030284/1‑138 (MQ=255)
                                                                                                          gAAGGATTGGTGAATCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAAtct                  <  2:710718/139‑1 (MQ=255)
                                                                                                                   gTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGa         <  1:788524/139‑1 (MQ=255)
                                                                                                                        tCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATa                                                                                  >  1:428825/1‑61 (MQ=255)
                                                                                                                        tCAGTTTAAAACCGCCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATa                                                                                  <  2:428825/61‑1 (MQ=255)
                                                                                                                          aGTTAAAAACCGCCGATCAATAGCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGGGCCAGCGTTGAGCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGAGt  <  2:275780/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
TACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGT  >  NZ_CP009273/1134139‑1134399

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