Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,978,163 A→G 12.0% V161A (GTG→GCG)  flgH ← flagellar L‑ring protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,978,1630AG12.0% 126.8 / 7.4 50V161A (GTG→GCG) flgHflagellar L‑ring protein
Reads supporting (aligned to +/- strand):  ref base A (19/25);  new base G (3/3);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAGCGAAAGGTCACCGCCACCGAACGG  >  NC_002947/4978021‑4978311
                                                                                                                                               |                                                                                                                                                    
tgTACCGGGCCAGGTGAGGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTCTCTTCACGCTGCTGCCAATTTCCAACCATGCCGCTTGTGCGACCCAGCTCTTCGCCGGTGTTCAGCGTCAGCCAATTCTCGCCGCGCACGGCgat                                                                                                                                               <  1:876027/149‑1 (MQ=255)
ggTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCgat                                                                                                                                               <  1:497771/150‑1 (MQ=255)
    ccGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGatgatg                                                                                                                                           <  2:103067/150‑1 (MQ=255)
              tGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCa                                                                                                                                 <  1:794769/150‑1 (MQ=255)
              tGATGCGCGCGTCGGCCACACCGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCa                                                                                                                                 <  2:510180/150‑1 (MQ=255)
                            gCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCacg                                                                                                                   <  1:1048131/150‑1 (MQ=255)
                               acacGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCacggtg                                                                                                                >  2:655893/1‑150 (MQ=255)
                                     gTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCacggtgacggtg                                                                                                          >  2:1152687/1‑150 (MQ=255)
                                              aCGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCg                                                                                                 >  1:587911/1‑150 (MQ=255)
                                              aCGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCg                                                                                                 >  1:629544/1‑150 (MQ=255)
                                                ggTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTggg                                                                                                                                   <  2:1157931/114‑1 (MQ=255)
                                                ggTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTggg                                                                                                                                   >  1:1157931/1‑114 (MQ=255)
                                                 gtgtTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGatgatg                                                                                                                                           <  2:200474/105‑1 (MQ=255)
                                                 gtgtTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGatgatg                                                                                                                                           >  1:200474/1‑105 (MQ=255)
                                                               gATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGGCAGCCCCCTCTCTCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTg                                                                                >  2:704091/1‑150 (MQ=255)
                                                                     gTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCacggtga                                                                                                               >  1:694940/1‑113 (MQ=255)
                                                                     gTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCacggtga                                                                                                               <  2:694940/113‑1 (MQ=255)
                                                                          gcgcGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTg                                                                      <  1:1152687/150‑1 (MQ=255)
                                                                           cgcgAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCg                                                                                                                         >  1:198244/1‑98 (MQ=255)
                                                                           cgcgAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCg                                                                                                                         <  2:198244/98‑1 (MQ=255)
                                                                            gcgAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCt                                                                    <  2:587911/150‑1 (MQ=255)
                                                                            gcgAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCt                                                                    >  2:224442/1‑150 (MQ=255)
                                                                                 ccATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGc                                                               <  1:9055/150‑1 (MQ=255)
                                                                                   aTGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTc                                                                                                                          <  1:107303/89‑1 (MQ=255)
                                                                                   aTGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGccc                                                             <  1:655893/150‑1 (MQ=255)
                                                                                   aTGCCGGCAATGCGAACAAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTc                                                                                                                          >  2:107303/1‑89 (MQ=255)
                                                                                      ccGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTg                                                          <  2:650559/150‑1 (MQ=255)
                                                                                       cccccaTGCGCACCACCCCTTCGCCTGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACCGCGATGATGACGGTGGACAGTACCTCGGCCACGGTGACGGTGATCTATCCTGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTgg                                                         <  2:401863/145‑1 (MQ=255)
                                                                                          cAATGCGCACCAGCTCTTCGCCGGTTTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTcg                                                      >  2:214131/1‑150 (MQ=255)
                                                                                             tGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTcgcgc                                                   <  2:1134392/150‑1 (MQ=255)
                                                                                                     aGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGc                                                                     <  1:1040480/124‑1 (MQ=255)
                                                                                                     aGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGc                                                                     >  2:1040480/1‑124 (MQ=255)
                                                                                                          ttCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCTTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCtata                                      >  2:1149199/1‑150 (MQ=255)
                                                                                                              ccGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCAGGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCg                                  <  1:472843/150‑1 (MQ=255)
                                                                                                              ccGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTgg                                                         >  2:853576/1‑127 (MQ=255)
                                                                                                              ccGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTgg                                                         <  1:853576/127‑1 (MQ=255)
                                                                                                              ccGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCg                                  >  2:353007/1‑150 (MQ=255)
                                                                                                                     tCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCa                           <  2:680786/150‑1 (MQ=255)
                                                                                                                           tCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCt                                                            >  2:962761/1‑111 (MQ=255)
                                                                                                                           tCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCt                                                            <  1:962761/111‑1 (MQ=255)
                                                                                                                               ccACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTgg                                                         <  2:50978/110‑1 (MQ=255)
                                                                                                                               ccACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTgg                                                         >  1:50978/1‑110 (MQ=255)
                                                                                                                                  atcatcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAg                          >  1:762802/5‑138 (MQ=255)
                                                                                                                                  atcatcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAg                          <  2:762802/134‑1 (MQ=255)
                                                                                                                                    catcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGcc                                                              <  2:543713/98‑1 (MQ=255)
                                                                                                                                    catcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGcc                                                              >  1:543713/3‑100 (MQ=255)
                                                                                                                                    catcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCta                                        <  1:1106947/120‑1 (MQ=255)
                                                                                                                                    catcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCta                                        >  2:1106947/3‑122 (MQ=255)
                                                                                                                                        gccgcGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAGCGAAAGGTCACCGCCAcc        >  1:1012074/1‑150 (MQ=255)
                                                                                                                                              cACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAGCGAAAGGTCACCGCCACCGAACgg  <  2:1012074/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAGCGAAAGGTCACCGCCACCGAACGG  >  NC_002947/4978021‑4978311

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