Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,972,691 A→G intergenic (‑237/‑543) flhD ← / → uspC flagellar transcriptional regulator FlhD/universal stress protein UspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,972,6910AG95.5% 320.7 / 7.2 89intergenic (‑237/‑543)flhD/uspCflagellar transcriptional regulator FlhD/universal stress protein UspC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (45/40);  minor base . (2/2);  total (47/42)
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.35e-01

TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGT  >  NZ_CP009273/1972605‑1972777
                                                                                       |                                                                                      
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                     >  1:793202/1‑90 (MQ=255)
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                     >  1:168949/1‑90 (MQ=255)
      tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTg                                                                               >  2:172203/1‑90 (MQ=255)
       aGAAATGGGTTTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGc                                                                              <  2:2403781/90‑1 (MQ=255)
       aGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGc                                                                              >  1:2429817/1‑90 (MQ=255)
             gggTCTTTACACTTATTTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGCTGCGATTTTc                                                                       <  2:2829250/90‑1 (MQ=255)
               gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCa                                                                      >  1:2377735/1‑90 (MQ=255)
                 cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaat                                                                    >  2:95375/1‑90 (MQ=255)
                 cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaat                                                                    >  2:922307/1‑90 (MQ=255)
                  tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaata                                                                   >  2:2129948/1‑90 (MQ=255)
                  tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaata                                                                   >  1:1125812/1‑90 (MQ=255)
                  tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaata                                                                   >  1:2690286/1‑90 (MQ=255)
                   ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTCCTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                                  <  2:1140562/90‑1 (MQ=255)
                   ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                                  <  2:2429817/90‑1 (MQ=255)
                     acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATg                                                                <  2:1318232/90‑1 (MQ=255)
                        cTTATCTAAGATTTTTCCTAAATCGACGTAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                             <  1:868969/90‑1 (MQ=255)
                        cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                             <  1:1207238/90‑1 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                            >  1:2823943/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                            >  1:2577309/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                            >  2:1925466/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                            >  2:1989253/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                            >  2:2594514/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                            >  2:2644517/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGCTGCGATTTTCAATAATGCGTGa                                                           >  2:1523875/1‑90 (MQ=255)
                          tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                                  <  2:477190/83‑1 (MQ=255)
                          tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                                  >  1:477190/1‑83 (MQ=255)
                           aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGAt                                                          >  1:993924/1‑90 (MQ=255)
                            tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATg                                                         <  1:197115/90‑1 (MQ=255)
                                aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGa                                                     >  2:2840198/1‑90 (MQ=255)
                                 gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAt                                                    <  2:989242/90‑1 (MQ=255)
                                 gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAt                                                    <  2:1764339/90‑1 (MQ=255)
                                    ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                                  <  2:793202/90‑1 (MQ=255)
                                    ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcac                                                 >  1:2957494/1‑90 (MQ=255)
                                     tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcaca                                                >  2:2421084/1‑90 (MQ=255)
                                     tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcaca                                                >  1:620582/1‑90 (MQ=255)
                                     tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcaca                                                >  2:2565514/1‑90 (MQ=255)
                                     tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcacac                                               >  2:2234685/1‑90 (MQ=255)
                                        ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACaa                                             <  1:1716192/90‑1 (MQ=255)
                                         cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACaaa                                            <  2:391203/90‑1 (MQ=255)
                                           aaaTCGACGCAACTGTACTCGTCACTACACGCCCATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAac                                          <  1:2648740/90‑1 (MQ=255)
                                           aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAac                                          >  2:177805/1‑90 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAt                                   <  2:477620/90‑1 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAt                                   <  1:1527062/90‑1 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAt                                   <  2:2585604/90‑1 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAtt                                  <  1:1367822/90‑1 (MQ=255)
                                                   gCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAtt                                  <  1:1289605/90‑1 (MQ=255)
                                                   gCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAtt                                  <  1:2780820/90‑1 (MQ=255)
                                                           aCTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                          >  2:2023340‑M1/1‑88 (MQ=255)
                                                           aCTCGTCACTCACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                         <  2:552880‑M1/90‑4 (MQ=255)
                                                             tCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                        <  2:2947471‑M1/90‑5 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGa                                                                            <  2:2782634/37‑1 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGa                                                                            >  1:2782634/1‑37 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                                  >  2:1677226/1‑47 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                                  <  1:1677226/47‑1 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                                  >  1:1380418/1‑63 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                                  <  2:1380418/63‑1 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                       >  1:165824‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                       >  1:2258103‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                       >  1:669808‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                       >  2:2405100‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                       >  2:525083‑M1/1‑85 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                     >  1:136054‑M1/1‑83 (MQ=255)
                                                                 cACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAac                                          >  2:2559528/1‑68 (MQ=255)
                                                                 cACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAac                                          <  1:2559528/68‑1 (MQ=255)
                                                                   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                  >  2:1815862‑M1/1‑80 (MQ=255)
                                                                   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgc                   <  2:2436701‑M1/89‑10 (MQ=255)
                                                                   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgc                   >  1:2436701‑M1/1‑80 (MQ=255)
                                                                     acacGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaa                >  2:946029‑M1/1‑78 (MQ=255)
                                                                        cGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                                  <  2:241834/53‑1 (MQ=255)
                                                                        cGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                                  >  1:241834/1‑53 (MQ=255)
                                                                        cGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa             <  1:2840198‑M1/90‑16 (MQ=255)
                                                                          cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt           <  1:2565514‑M1/90‑18 (MQ=255)
                                                                          cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt           >  1:1208170‑M1/1‑73 (MQ=255)
                                                                          cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt           <  2:2823943‑M1/90‑18 (MQ=255)
                                                                          cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt           >  2:1594678‑M1/1‑73 (MQ=255)
                                                                          cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggggcgaacaagt           >  2:1763717‑M1/1‑73 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        <  2:993924‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        <  2:631902‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        >  1:2901193‑M1/1‑70 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        <  1:2071649‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        <  1:2023340‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        >  2:313589‑M1/1‑70 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        <  1:1925466‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        <  2:2300232‑M1/90‑21 (MQ=255)
                                                                              tACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccct       <  2:2625391‑M1/90‑22 (MQ=255)
                                                                               aCAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        >  1:366572‑M1/1‑68 (MQ=255)
                                                                               aCAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc        <  2:366572‑M1/88‑21 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata   >  1:1812223‑M1/1‑65 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata   >  1:1362449‑M1/1‑65 (MQ=255)
                                                                                   cGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatat  <  2:75624‑M1/90‑27 (MQ=255)
                                                                                   cGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatat  <  2:1362449‑M1/90‑27 (MQ=255)
                                                                                       |                                                                                      
TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGT  >  NZ_CP009273/1972605‑1972777

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