Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,014,103 A→C Y179S (TAC→TCC)  fliM → flagellar motor switch protein FliM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,014,1030AC88.5% 72.9 / 4.1 26Y179S (TAC→TCC) fliMflagellar motor switch protein FliM
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (14/9);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

GGTCGCGAGTTTACCCATACCGAACAGCGCGTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTACGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATT  >  NZ_CP009273/2013994‑2014229
                                                                                                             |                                                                                                                              
ggTCGCGAGTTTACCCATACCGAACAGCGCGTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTACGTGCGTTCGGAAATGCAGGTGAAATTTa                                                                                                   >  2:207854/1‑139 (MQ=255)
  tCGCGAGTTTACCCATACCGAACAGCGCGTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTAcc                                                                                                 <  2:157444/139‑1 (MQ=255)
     cgAGTTTACCCATACCGAACAGCGCGTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCaa                                                                                               >  1:592648/1‑138 (MQ=255)
                   ccGAACAGCGCGTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTACGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGcc                                                                                <  1:413933/139‑1 (MQ=255)
                      aaCAGCGCGTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGaa                                                                             >  1:183739/1‑139 (MQ=255)
                      aaCAGCGCGTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGaa                                                                             >  2:161774/1‑139 (MQ=255)
                           cgcgTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACa                                                                        >  2:732672/1‑139 (MQ=255)
                           cgcgTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACa                                                                        >  2:538349/1‑139 (MQ=255)
                           cgcgTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACa                                                                        >  1:588821/1‑139 (MQ=255)
                             cgTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTACGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACAtt                                                                      <  1:207854/139‑1 (MQ=255)
                                catcaACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATtgtg                                                                   >  1:45706/1‑139 (MQ=255)
                                                aaaCTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTc                                                    >  2:683196/1‑138 (MQ=255)
                                                      gcgcTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCAtgtg                                              <  2:225132/138‑1 (MQ=255)
                                                      gcgcTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCAtgtg                                              >  2:384180/1‑138 (MQ=255)
                                                      gcgcTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCAtgtg                                              >  1:513327/1‑138 (MQ=255)
                                                      gcgcTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCAtgtg                                              >  1:724166/1‑138 (MQ=255)
                                                       cgcTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGcc                                                        <  1:562507/127‑1 (MQ=255)
                                                       cgcTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGcc                                                        >  2:562507/1‑127 (MQ=255)
                                                                            ccTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACa                                                                        <  2:292944/90‑1 (MQ=255)
                                                                            ccTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACa                                                                        >  1:292944/1‑90 (MQ=255)
                                                                            ccTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCg                        >  1:479324/1‑138 (MQ=255)
                                                                            ccTGGAAGGCGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAAGGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCg                        <  2:479324/138‑1 (MQ=255)
                                                                                     cGATTAATCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAatat              <  2:183739/139‑1 (MQ=255)
                                                                                            tCCGCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATctgcctg       <  2:202028/139‑1 (MQ=255)
                                                                                               gCTGGAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCa    <  1:550305/139‑1 (MQ=255)
                                                                                                   gAAGTTGAGTCCGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCAtt  <  2:513327/137‑1 (MQ=255)
                                                                                                             |                                                                                                                              
GGTCGCGAGTTTACCCATACCGAACAGCGCGTCATCAACCGCATGTTGAAACTGGCGCTTGAAGGCTATAGCGACGCCTGGAAGGCGATTAATCCGCTGGAAGTTGAGTACGTGCGTTCGGAAATGCAGGTGAAATTTACCAATATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATT  >  NZ_CP009273/2013994‑2014229

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