Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,139,406 Δ1 bp coding (1029/1644 nt) flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,139,4060T.100.0% 215.4 / NA 53coding (1029/1644 nt)flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (22/31);  total (22/31)

GGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGT  >  NC_000913/1139272‑1139528
                                                                                                                                      |                                                                                                                          
ggACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGt                                                                                                                       >  1:901728/1‑139 (MQ=255)
      cTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                 >  1:1315233/1‑139 (MQ=255)
        ggCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCg                                                                                                               >  2:34699/1‑139 (MQ=255)
          cgcTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGt                                                                                                             >  2:230353/1‑139 (MQ=255)
          cgcTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGt                                                                                                             <  1:547417/139‑1 (MQ=255)
           gcTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTg                                                                                                            <  1:507853/139‑1 (MQ=255)
                cATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACg                                                                                                       <  1:1411223/139‑1 (MQ=255)
                      ccGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACt                                                                                                 <  1:1643290/139‑1 (MQ=255)
                         aGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGAt                                                                                              >  1:1265870/1‑139 (MQ=255)
                           gCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGc                                                                                            <  2:1362831/139‑1 (MQ=255)
                                cAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcg                                                                                       <  2:901728/139‑1 (MQ=255)
                                    acCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGt                                                                                   <  1:1197665/139‑1 (MQ=255)
                                     cccAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTa                                                                                  >  1:686457/1‑139 (MQ=255)
                                            caAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGa                                                                           <  1:1329811/139‑1 (MQ=255)
                                            caAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGa                                                                           <  1:338660/139‑1 (MQ=255)
                                                    gATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                       >  1:1475213/1‑135 (MQ=255)
                                                    gATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                       <  2:1475213/135‑1 (MQ=255)
                                                               aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                        <  2:1523666/139‑1 (MQ=255)
                                                               aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                        <  2:1265870/139‑1 (MQ=255)
                                                               aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                        <  2:161725/139‑1 (MQ=255)
                                                               aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                        <  1:510788/139‑1 (MQ=255)
                                                               aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                        <  2:446613/139‑1 (MQ=255)
                                                                 cGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCg                                                      >  2:608561/1‑139 (MQ=255)
                                                                    cGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGata                                                   >  1:582586/1‑139 (MQ=255)
                                                                      aTGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGataat                                                 >  2:1165191/1‑139 (MQ=255)
                                                                        gCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGc                                                                                            <  1:1140639/94‑1 (MQ=255)
                                                                        gCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGc                                                                                            >  2:1140639/1‑94 (MQ=255)
                                                                         ccGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATaataat                                              <  1:608561/139‑1 (MQ=255)
                                                                                aGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGc                                       >  2:348131/1‑139 (MQ=255)
                                                                                 gATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCa                                      <  2:21565/139‑1 (MQ=255)
                                                                                         tGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCAccc                              >  1:1030446/1‑139 (MQ=255)
                                                                                           cTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGc                            <  1:34699/139‑1 (MQ=255)
                                                                                              tCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCtt                         <  2:1030446/139‑1 (MQ=255)
                                                                                                     gCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                 <  1:393293/44‑1 (MQ=255)
                                                                                                     gCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                 >  2:393293/1‑44 (MQ=255)
                                                                                                     gCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGataata                                                >  2:1474630/1‑109 (MQ=255)
                                                                                                     gCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGataata                                                <  1:1474630/109‑1 (MQ=255)
                                                                                                       ccGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAAt                <  2:1635540/139‑1 (MQ=255)
                                                                                                        cgcgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCg                                                          >  1:43374/1‑96 (MQ=255)
                                                                                                        cgcgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCg                                                          <  2:43374/96‑1 (MQ=255)
                                                                                                         gcgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATAc              <  1:230353/139‑1 (MQ=255)
                                                                                                          cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGc                                                                                                                    >  1:1090026/1‑36 (MQ=255)
                                                                                                          cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGc                                                                                                                    <  2:1090026/36‑1 (MQ=255)
                                                                                                          cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGa                                                                                               <  1:454097/57‑1 (MQ=255)
                                                                                                          cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGa                                                                                               >  2:454097/1‑57 (MQ=255)
                                                                                                          cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                       <  1:6768/81‑1 (MQ=255)
                                                                                                          cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                       >  2:6768/1‑81 (MQ=255)
                                                                                                           ggTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTg                        >  2:387990/1‑127 (MQ=255)
                                                                                                           ggTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTg                        <  1:387990/127‑1 (MQ=255)
                                                                                                                  cAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTAc     >  1:390926/1‑139 (MQ=255)
                                                                                                                     aaCACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGt  <  1:348131/139‑1 (MQ=255)
                                                                                                                            aaaaCAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGa                                                                                               <  1:678608/39‑1 (MQ=255)
                                                                                                                            aaaaCAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGa                                                                                               >  2:678608/1‑39 (MQ=255)
                                                                                                                                      |                                                                                                                          
GGACAACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGT  >  NC_000913/1139272‑1139528

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