Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,978,160 C→T 12.4% R162H (CGC→CAC)  flgH ← flagellar L‑ring protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,978,1600CT12.4% 120.1 / 7.1 48R162H (CGC→CAC) flgHflagellar L‑ring protein
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base T (3/3);  total (21/27)
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.76e-01

AACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAGCGAAAGGTCACCGCCACC  >  NC_002947/4978013‑4978305
                                                                                                                                                    |                                                                                                                                                 
aacgaacCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCgcgc                                                                                                                                                 >  1:680786/1‑150 (MQ=255)
        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:629544/1‑150 (MQ=255)
                                                      aCGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCg                                                                                           >  1:587911/1‑150 (MQ=255)
                                                        ggTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTggg                                                                                                                             <  2:1157931/114‑1 (MQ=255)
                                                        ggTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTggg                                                                                                                             >  1:1157931/1‑114 (MQ=255)
                                                         gtgtTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGatgatg                                                                                                                                     >  1:200474/1‑105 (MQ=255)
                                                         gtgtTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGatgatg                                                                                                                                     <  2:200474/105‑1 (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                                                   <  1:853576/127‑1 (MQ=255)
                                                                                                                      ccGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTgg                                                   >  2:853576/1‑127 (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                                                        >  1:543713/3‑100 (MQ=255)
                                                                                                                                            catcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGcc                                                        <  2:543713/98‑1 (MQ=255)
                                                                                                                                            catcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCta                                  <  1:1106947/120‑1 (MQ=255)
                                                                                                                                            catcGCCGTGCGCGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCta                                  >  2:1106947/3‑122 (MQ=255)
                                                                                                                                                gccgcGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAGCGAAAGGTCACCGCCAcc  >  1:1012074/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGGTCAGCCACTTCTCGCCGCGCACGGCGATGATGCCGTTGGGCAGTACCTCGGCCACGGTGACGGTGATCGAACCGGTCAGCGTGTTGCCCTGGGTGGCCTTGCTGTCGCCCTTGGTCGCGCGATCGCCGCTATAACCGGCCTCCAGCGAAAGGTCACCGCCACC  >  NC_002947/4978013‑4978305

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: