Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,136,783 T→C V170V (GTT→GTC)  flgI → putative flagellar basal body protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,136,7830TC100.0% 204.2 / NA 54V170V (GTT→GTC) flgIputative flagellar basal body protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/27);  total (27/27)

CGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCACCGCG  >  NC_000913/1136647‑1136918
                                                                                                                                        |                                                                                                                                       
cGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCAt                                                                                                                                       <  2:606558/139‑1 (MQ=255)
cGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCAt                                                                                                                                       >  1:346039/1‑139 (MQ=255)
          cGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCCGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAAc                                                                                                                                  >  1:1275589/1‑134 (MQ=255)
          cGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAAc                                                                                                                                  <  2:1275589/134‑1 (MQ=255)
                            tGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGc                                                                                                           <  1:1476580/139‑1 (MQ=255)
                              cgcTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGt                                                                                                         <  1:1307839/139‑1 (MQ=255)
                               gcTGGCGTAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAAGTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTc                                                                                                        >  2:1093781/1‑139 (MQ=255)
                                    cgcAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   <  1:849063/139‑1 (MQ=255)
                                    cgcAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   <  1:1429453/139‑1 (MQ=255)
                                                      ttGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACt                                                                                 <  2:346039/139‑1 (MQ=255)
                                                        ggcggcGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCccagc                                                                                                                    >  2:536740/1‑102 (MQ=255)
                                                        ggcggcGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCccagc                                                                                                                    <  1:536740/102‑1 (MQ=255)
                                                          cggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGaa                                                                                                                                   >  2:708440/1‑85 (MQ=255)
                                                          cggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGaa                                                                                                                                   <  1:708440/85‑1 (MQ=255)
                                                          cggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCCGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAac                                                                             >  1:801033/1‑139 (MQ=255)
                                                          cggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAac                                                                             >  2:1533209/1‑139 (MQ=255)
                                                          cggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAac                                                                             >  2:756015/1‑139 (MQ=255)
                                                           ggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAacg                                                                            <  1:614010/139‑1 (MQ=255)
                                                                          gCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCa                                                             <  1:8438/139‑1 (MQ=255)
                                                                            tGGCAGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGaa                                                                                                                                   <  2:795068/67‑1 (MQ=255)
                                                                            tGGCAGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGaa                                                                                                                                   >  1:795068/1‑67 (MQ=255)
                                                                                            cAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTg                                           <  1:756015/139‑1 (MQ=255)
                                                                                             aGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGa                                          <  1:1093781/139‑1 (MQ=255)
                                                                                               gTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGaca                                        >  2:636030/1‑139 (MQ=255)
                                                                                                  aaCCAACTGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCa                                     >  2:1607036/1‑139 (MQ=255)
                                                                                                         tGAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACcg                              <  1:636030/139‑1 (MQ=255)
                                                                                                          gAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGccc                                                                                                                       >  1:1458960/1‑49 (MQ=255)
                                                                                                          gAACGGTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGccc                                                                                                                       <  2:1458960/49‑1 (MQ=255)
                                                                                                              ggTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                         <  1:1607036/139‑1 (MQ=255)
                                                                                                              ggTGGACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCgcgtgc                         >  2:1444171/1‑139 (MQ=255)
                                                                                                                  gACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATgg                                                          <  2:336404/102‑1 (MQ=255)
                                                                                                                  gACGGATCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATgg                                                          >  1:336404/1‑102 (MQ=255)
                                                                                                                       aTCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTg                                           <  2:291111/112‑1 (MQ=255)
                                                                                                                       aTCACCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTg                                           >  1:291111/1‑112 (MQ=255)
                                                                                                                          aCCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGt                                                                                                         >  2:1741561/1‑47 (MQ=255)
                                                                                                                          aCCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGt                                                                                                         <  1:1741561/47‑1 (MQ=255)
                                                                                                                          aCCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   >  1:612618/1‑53 (MQ=255)
                                                                                                                          aCCAATGGTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   <  2:612618/53‑1 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   >  2:310163/1‑46 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   <  1:310163/46‑1 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGaa                                                                       <  1:424513/74‑1 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGaa                                                                       >  2:424513/1‑74 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGa                                          <  2:799501/103‑1 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGa                                          >  1:799501/1‑103 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCAc      >  2:947917/1‑139 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCAc      >  2:632273/1‑139 (MQ=255)
                                                                                                                                 gTGCGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCAc      >  1:693461/1‑139 (MQ=255)
                                                                                                                                 gTGCGGACATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGTAAATCGCTg                                           >  1:1382545/1‑102 (MQ=255)
                                                                                                                                 gTGCGGACATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGTAAATCGCTg                                           <  2:1382545/102‑1 (MQ=255)
                                                                                                                                    cGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   <  2:14397/43‑1 (MQ=255)
                                                                                                                                    cGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   >  1:14397/1‑43 (MQ=255)
                                                                                                                                    cGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   <  1:914323/43‑1 (MQ=255)
                                                                                                                                    cGGTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGaa                                                                                                   >  2:914323/1‑43 (MQ=255)
                                                                                                                                     ggTCATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCACcgcg  >  1:801840/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGATTTCAGCATGGCGCAGCAAATCGCTGACACCATCAACCGCGTGCGTGGATATGGCAGCGCCACCGCG  >  NC_000913/1136647‑1136918

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑